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INPUT FILTER COMPENSATION FOR 
SWITCHING REGULATORS 

by 

Shriram S. Kelkar 
(ABSTRACT) 

An input filter is often required between a switching regu- 
lator and its power source due to the need of preventing the 
regulator switching current from being reflected bade into 
the source. The presence of the input filter often results 
in various performance difficulties such as loop instabili- 
ty, degradation of transient response, audiosusceptibility 
and output impedance characteristics. These problems are 
cause.- mainly by the interaction between, the peaking of the 
output impedance of the input filter and the regulator con- 
trol loop. Conventional single-stage and two-stage input 
filters can be designed to minimize the peaking effect, how- 
ever this often resuu in a penalty of weight or loss in- 
crease in the input filter. 

A novel input filter compensation scheme for a buck re- 
gulator that eliminates the interaction between the input 
filter output impedance and the regulator control loop is 
presented. The scheme is implemented using a feedforward 


loop that senses the input filter state variables and uses 
this information to modulate the duty cycle signal. The 
feedforward design process presented is seen to be straight- 
forward and the feedforward easy to implement. Extensive ex- 
perimental data supported by analytical results show that 
significant performance improvement is achieved with the use 
of feedforward in the following performance categories: loop 
stability, audiosusceptibility, output impedance and tran- 
sient response. 

The use of feedforward results in isolating the switch- 
ing regulator from its power source thus eliminating all in- 
teraction between the regulator and equipment upstream. In 
addition the use of feedforward removes some of the input 
filter design constraints and makes the input filter design 
process simpler thus making it possible to optimize the in- 
put filter. The concept of feedforward compensation can also 
be extended to other types of switching regulators. 
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Chapter I 
INTRODUCTION 


Switching-mode dc-dc regulators are coming into increasing 
use as power supplies because of the significant reduction 
of weight/ size and increase in equipment efficiency that can 
be attained. The dissipative regulators used earlier had 
the advantage of design simplicity/ but suffer from low ef- 
ficiency. and higher weight due to the large transformer and 
filter needed to achieve the required voltage/current regu- 
lation. 

Since the late sixties, the rapid expansion of the com- 
puter and communication industries and increasing complexity 
and sophistication of various systems necessitated the de- 
velopment of higher performance switched mode power sup- 
plies. The incentive for performance improvement prompted 
the initial development of multiple loop control schemes, 
such as the standardized control module (SCM) for dc-dc con- 
verters [2,7,8] and the current-injected control scheme [9]. 
The ever continuing search for performance improvement forms 
the underlying theme of this dissertation. 

The work presented in this dissertation is mainly con- 
cerned with developing a control scheme to alleviate the 
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problem brought about due to the use of an input filter in 
switching regulators. The switching regulator input current 
has a substantial ripple component at the switching frequen- 
cy and this necessitates the use of an input filter to 
smooth the pulsating current drawn from the supply. Furth- 
ermore the input filter also serves to attenuate noise pre- 
sent in the supply voltage from being propagated through the 
regulator to the payload downstream. The presence of the 
input filter, however, often results in various performance 
difficulties such as loop instability, degradation of tran- 
sient response, audiosusceptibility (closed loop input-to- 
output gain) and output impedance characteristics 
[3,4:, 5, 6, 10] . These problems are caused mainly by the inter- 
action between the resonant peaking of the output impedance 
of the input filter and the regulator control loop. Conven- 
tional single-stage and two-stage input filters can be de- 
signed such that the peaking effect is minimized, however 
such a design is often accompanied with a penalty of weight 
and loss in the input filter [3,4,5,6,10]. This disserta- 
tion presents a different approach via a feedforward control 
scheme to mitigate the undesirable interaction between the 
input filter and the regulator control loop. 

The concept of pole-zero cancellation is used, in this 
dissertation, to develop a novel feedforward control loop 
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that senses the input filter state variables and processes 
this information in a manner designed to cancel the detri- 
mental effect of peaking of the output impedance of the in- 
put filter. The feedforward loop working in conjunction 
with the feedback loops developed earlier [2,7,81 constitute 
a total state control scheme that eliminates the interaction 
between the input filter and the regulator control loop. Em- 
ploying the novel feedforward compensation scheme presented 
in this dissertation, a high performance converter together 
with an effective input filter design (minimum weight and 
loss) can be accomplished concurrently. A buck regulator 
employing a feedforward control loop working in conjunction 
with the feedback loops was used to obtain measurement - j that' 
showed significant improvement in the following performance 
categories: 

1. Eoop stability (open loop gain and phase margins); 

2 . Audiosusceptibi lity; 

3. Output impedance; and 

4. Transient response. 

In this dissertation the problem caused by input filter 
interaction and conventional input filter design techniques 
are discussed in Chapter 2 followed in Chapter 3 by the con- 
cept of pole-zero cancellation developed earlier [2,7,8]. 
Chapters 4 and 5 present the modeling of the power stage 


with input filter and the implementation of the feedforward 
for a buck regulator respectively. Measurements of open loop 
gain and phase that confirm the analytical prediction of 
performance improvement using the feedforward scheme devel- 
oped, are discussed in Chapter 6 along with other measure- 
ments of audiosusceptibility, output impedance and transient 
response. Chapter 7 presents experimental and analytical re- 
sults pertaining to transient response while Chapter 8 dis- 
cusses the use of the feedforward loop in stabilizing a re- 
gulator system made unstable due to input filter 
interaction. Chapter 9 is concerned with extending the con- 
cept of feedforward compensation to other types of control 
and to other types of regulators i.e. the buck-boost regula- 
tors, Finally, Chapter 10 presents the conclusions and sug- 
gestions for future work. 


Chapter II 


INPUT FILTER RELATED PROBLEMS AND CONVENTIONAL 
DESIGN TECHNIQUES 


2.1 INPUT FILTER RELATED PROBLEMS 

An input filter is often required between a switching regu- 
lator and its power source. A buck type switching regulator 
with a single-stage input filter is shown in Figure 1 . The 
regulator input current has a substantial pulsating current 
component at the switching frequency as a result of the 
opening and closing of the switch and this component should 
be prevented from being reflected back into the source ; an 
input filter is required to provide high attenuation at 
switching frequency and thus smooth the current drawn from 
the source. The input filter also serves to isolate source 
voltage disturbances from being propagated to the switching 
regulator payload downstream. 

A presumably well-designed input filter, satisfying the 
above mentioned requirements, when used with a switching re- 
gulator can often cause significant performance degradations 
[3,4,5,6,10]. This is due primarily to the complex interac- 
tion between the switching regulator control loop, the input 
filter and the regulator output filter [3,4,5,6,10]. 
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Figure 1: Suck Converter with an Input Filter 


5 V 0 ~5v R L =0.86 ohm L=200 microH 

R x =20mohin C=1540microF R c =7mohm T p =50microsecs 
Lj-77 microH Rj=39.6 mohm C, =412 microF 


Cj =412 irdcroF 



The interaction between the control loop and the input 

filter is illustrated in Figure 2 . The switching regulator 

has been shown to have a nonlinear negative resistance, as 

illustrated in the figure, [3]. The input current i r to the 

switching regulator is related nonlinear ly to the input vol- 

diy> i 

tage e r and- the input resistance --- « -- . Under certain 

afij* 

conditions the input filter- switching regulator combination 
can become a negative resistance oscillator, producing large 
amplitude voltage excursions across capacitor C. When this 
happens serious degradation of regulator performance could 
occur, [3J, including loss of stability. 

The effect of the input filter is more clearly seen us- 
ing a small signal model. 

The averaging technique [1] is used to relate the low 
frequency modulation component of the source voltage and 
control signal to the corresponding frequency components of 
the converter output voltage. Using the continuous inductor 
current buck regulator with input filter of Figure 1 , as an 
example, a small signal model using the dual-input describ- 
ing function can be developed, as shown in Figure 3 , [6]. 
In this model the effect of the input filter is character- 
ised by the following two parameters: the forward transfer 
characteristic of the input filter H(s) and the output impe- 
dance of the input filter Z(s). 
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In Figure 1 the output filter is made up of R^, L, R c and C, 
R l is the load resistance, D is the steady state duty cycle 
ratio D = T oa /T, V E and Ij are the steady state regulator 
input voltage and current respectively, and the lower case 
letters with a caret above them denote modulation signals. 

The small signal model of Figure 3 is used to illus- 
trate briefly the complex interaction between the input fil- 
ter, output filter and the control loop and the problems 
caused by the interaction. For detailed analysis please re- 
fer to [ 6 ] . 

2.1.1 Input Filter Interaction — Loop Stability and 
Transient Response 

The stability of a switching regulator can be examined by 
the open loop gain G T ( s ) : 

G t (s) = F c (s)F p (s)F E (s)F M (s) (2-1) 

where F c (s)F p (s) is the duty cycle- to -output describing 
function v Q /d , and F E (s), F M (s) are the transfer functions 
of the error processor and the pulse modulator respectively. 
The peaking of the output impedance of the input filter 2(s) 
has the following effects: 

(1) The duty-cycle power stage gain F c (s) includes the 
output impedance Z(s) — 

F c (s) = Vj-Sfs)^ or 


(2-2) 
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F C <S) 



Z (s) 


(2-3) 


The first term in the brackets is the negative input 

impedance of the regulator. At the input filter resonant 
frequency, Z(s) reaches a peak value and if this value is 
large enough the result could be a reduction in loop gain or 
even worse a negative duty cycle power stage gain F c (s). Re- 
duction in loop gain could lead to loop instability, whereas 
a negative F c (s) together with the negative feedback loop 
will result in a positive feedback unstable system. 

(2) The power stage transfer function F p (s) includes the 
output impedance 2{s) — 

[R_ + l/sC]//R- 

F (s) = — ~ — — where (2-4) 

p D^2(s) + z i (s) 


2^(s) » R £ + sL + rR c + 1 /sC]//R l (2-5) 

= input impedance of the regulator. 

Excessive Z(s) at the input filter resonant frequency can 
significantly reduce Fp(s), and thus the loop gain. 

Figures 4 (61, illustrate the effect of peaking of Z(s) 
on the duty cycle- to-output transfer function F c (s)Fp{s) if 
an improperly designed input filter is employed. 
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At the input filter resonant frequency the peaking of the 
output impedance 2(s) causes a sharp change in the gain and 
phase of the duty cycle- to-output transfer function. This 
could result in loop instability and degradation of tran- 
sient response from a presumably well damped system to an 
oscillatory one; control of the peaking effect of the output 
impedance 3{s) is necessary to avoid these problems. 

2.1.2 Input Filter Interaction— Audiosusceptibility and 
Output Impedance 

The audiosusceptibility refers to the switching regula- 
tor's ability to attenuate small signal sinusoidal distur- 
bances present at the input so as not to affect the regulat- 
ed output voltage. The audiosusceptibility performance is of 
considerable importance, as the regulator generally shares 
the input bus with other on-line equipment. The operation of 
this equipment generates noise voltages on the input line 
which must be attenuated by the closed- loop regulator so 
that operation of the various paylcads at the regulator out- 
put will not be adversely affected. The audiosusceptibility 
is expressed in terms of the closed loop input- to- output 
transfer function G„(s): 

„ , , _ V s) P I< S) V S > _ F I ts) V s) (2-6) 

A lS) ~ Vjts) 1+F C (S)F (s)F e (S)F m Cs) 1 + G t (s> 
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where Fj(s) = DH(s) = input voltage gain of the power stage. 
G a (s) and thus the audio susceptibility are affected by the 
resonant peaking of the output impedance Z(s) and of the 
forward transfer function of the input filter with the regu- 
lator disconnected H(s), because Fj(s) is a function of H(s) 
whereas F c (s) and F p (s) are functions of Z(s). The reduc- 
tion of loop gain at the resonant frequency can thus severe- 
ly degrade the audiosusceptibility. Figure 5 , [6], illus- 
trates the audiosusceptibility of the buck regulator with 
and without an input filter. 

The output impedance of the regulator should be small 
so that the regulator behaves like an ideal voltage source, 
however the output impedance is increased by the peaking of 
the output impedance of the input filter. 


V s) 


V s) 


1 + G T (s> 


(2-7) 


where Z p £s) is the output impedance of the power stage with 
the control loop open. 

At the resonant frequency the output impedance of the regu- 
lator Z Q (s) is increased. This is a consequence of the loss 
of loop gain G T (s) as a result of peaking. 
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The peaking of H<s) and Z(s) of the input filter thus 
results in a reduction in the loop gain, this in turn af- 
fects stability, transient response, audiosusceptibility and 
the output impedance of the regulator. 

A buck type switching regulator with a two stage input 
filter is shown an Figure 6 . Figure 7 , [6], shows the mea- 
sured values of open loop gain and phase as a function of 
the frequency. ( The input filter damping resistance R D is 
not employed to purposely illustrate the effect of input 
filter interaction with the regulator control loop. ) Signi- 
ficant changes in the open loop gain and phase characteris- 
tics at the resonant frequencies of both the first stage and 
the second stage of the input filter are observed, £6J. Thus 
it is seen that the peaking of the output impedance at reso- 
nant frequency of the two stage input filter can also cause 
serious performance degradation. 

2,2 INPUT FILTER DESIGN CONSIDERATIONS 

The design of the input filter is made more complicated by 
the necessity of satisfying the following constraints, which 
result from the interaction between the input filter and the 
regulator control loop discussed in section 2.1 : 

1. The amount of regulator switching current reflected 
back into the source should be limited (conducted in- 
terference requi rement ) . 
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Figure 6: f^^^^ituhing regulator with a two stage 


Vj =40 v V Q =20v Rr=7 ohm Re =0. 05 ohm 

C=45SmicroF R c =0.068ohm T p =50microsecs 
^-1.1 mH Rt =0 . 05 mohm L a =0.27 tnH 
R*=0.02 ohm fc - 225 microF =20 micro F 
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Figure 7: Open. loop gain and phase characteristics without 

the input filter damping resistance R^. 

(parameter values used are as shown in Fig. 6) 
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2. The peaking of the output impedance o£ the input fil- 
ter 2(s) should be limited to a safe value to avoid 
significant loop gain reduction. 

3. The peaking of the transfer function of the input 
filter H(s) should be limited to achieve a satisfac- 
tory rejection rate of audio signals propagating from 
input to output. 

4. The input filter weight and energy loss should be 
limited to low values. 

5. The Nyquist stability criterion has to be satisfied; 
thus the closed loop poles should be in the left half 
plane for stable operation - 

U + F c (s)F p (s)F B (s)F M (s)l > 0 

6. The closed-loop input- to-output transfer characteris- 
tic (audiosusceptibility) and transient response due 
to a sudden line/ load change should not be degraded 
by a noticeable amount. 

An input filter design that satisfies one constraint 
may often result in violating some other constraint. For ex- 
ample , an input filter design that limits performance degra- 
dation (degradation of stability, transient response and au- 
diosusceptibility) often results in higher weight and 
increased losses in the input filter. A satisactory input 
filter design trades off one or more of the performance deg- 
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radations for size, weight and loss. A near- to- optimal de- 
sign thus requires many trial and error design attempts. 
Some of the conventional input filter design techniques are 
presented next. 

2.2.1 Conventional Input Filter Design Techniques 
A single stage input filter as shown in Figure 8(a) can be 
designed to avoid performance degradation — but this would 
result in larger filter Lj and Cj thus resulting in weight 
and size increase. The filter is simple and commonly used 
but it cannot often satisfy the stringent requirement on au- 
diosusceptibility without size/weight penalty. Resonant 
peaking of the filter of Figure 8(b)/ [3]/ is lowered by ad- 
ding resistance R, but this lowers efficiency because the 
pulse current flowing through Cj increases losses. Another 
design uses a resistance R in parallel across C 1 , [4], but 

this results in a large C 1 . 

The optimal design of a single stage input filter thus 
is rather difficult without tradeoff between performance 
degradations and the weight and loss limitations. 

The degradation of the power stage transfer function 
F p (s) due to peaking of 2{s) can be avoided if there is suf- 
ficient separation of the input filter resonance frequency 
61 1 - and the output filter resonant frequency 

^Tor 
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= -== * {4,5,10]. Figure 9 , [4,5,10], shows three pos- 

i/lTc 

sible combinations of c<) 0 and Wj . F p (s) is related to both 
2(s) and the input impedance of the regulator Z£(s) thus 


[R + 1/SCJ//R. 

p (S) = ~ (2-8) 

p uZ (s) + 2 i (s) 


where 3^(3) » + sL + [Rq * 1/sC]//R l 

= input impedance of the regulator. 

3(s) and 2j(s) peak at the frequencies to 0 and tdj respective- 
ly. If the two resonance frequencies are the same as in Fig- 
ure 9 then both Z(s) and 2j(s) peak at the same frequency 
and thus at that frequency the transfer function F p (s) would 
be affected, i.e. reduced, to the maximum possible extent. 
Shifting the two frequencies to 0 and i*)j apart as shown in 
Figure 9 will result in reducing the effect of peaking on 
F p (s). Reducing toj would result in increasing the size and 
weight of the input filter. A high value of tdj is desirable 
from the point of view of weight and size reduction but this 
can result in severe performance degradation -- from Figure 
4 it is clear that the gain of the duty cycle-to-output 
transfer function F c (s)F p (s) decreases with increasing fre- 
quency and thus the effect of peaking of Z(s) on the gain of 
E c (s)F p (s) would be more pronounced if 3{s) peaks at a high- 
er «i. 
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The reduction of the loop gain at higher input filter 
resonant frequency ojj often results in poor audio suscepti- 
bility, oscillatory transient response or even an unstable 
system. The choice of Uj thus involves a trade off between 
meeting performance specifications and size/weight. 


2.2.2 An Optimal Configuration 


A two- stage input filter configuration has been described, 
[3,6] and is shown in Eigure 10 . The first stage consisting 
of Lj,Cj,R 3 and Rj controls the resonant peaking of the fil- 
ter. The second stage consisting of L a , C 2 supplies most of 
the pulse current required by the regulator. As shown in the 
literature, [61, the two-stage input filter is capable of 
reducing H(s) and 2(s) at resonant frequency without signi- 
ficantly increasing weight and loss, unlike the single-stage 
input filter. Computer optimization techniques have .been 
utilized to optimally design the two-stage £ilter[6j. It has 
been shown that the two- stage filter is much lighter than 
its single-stage counterpart under identical design const- 
raints. Also it has been shown that for the same filter 
weight the single stage filter has a significantly higher 
peaking of H(s) and 2{s). Figure 11 shows the gain and phase 
of the duty cycle- to -output describing function of a power 
stage with a two stage input filter, [6j. 
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Figure 11: Duty cycle-to-output describing function of 

power stage with two-stage input filter. 

(parameter values used are as shown in Fig.l with the fol 
lowing two- stage input filter parameters -- 
Lj- 232 microH R, =0.0276 ohm L«j,=77 microH 

R z =0.0119 ohm C =1G0 microF R,=1.73 ohm 

C 2 =30 microF ) 1 * 
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Figure 4 shows the gain and phase of the duty 
cycle- to-output describing function of a power stage with a 
single-stage input f liter , and the two- stage filter of Fig- 
ure 11 was designed to have the same weight as the single- 
stage input filter of Figure 4 . Comparing the two figures 
the improvement in performance regarding the duty cycle-to- 
output transfer function is dramatic. 

It can therefore be concluded that the two- stage filter 
provides the best compromise among the conflicting require- 
ments of an input filter. 

2.2.3 Input Filter Compensation Via a Feedforward Loop 
Limiting interaction between the input filter and the regu- 
lator control loop is possible with the addition of a feed- 
forward control loop. At this point it is important to em- 
phasise that such a scheme is designed to eliminate the 
effect of peaking of the output impedance of the input fil- 
ter 2(s) , since the peaking of Z(s) interacts with the con- 
trol loop. The forward transfer function H(s) does not in- 
teract with the regulator control loop, as is evident from 
Figure 3 and therefore the peaking of H(s) cannot be cont- 
rolled in any way by adding a feedforward loop. The peaking 
of H(s) can only be controlled by proper filter design. In 
this dissertation a feedforward loop is implemented for a 


y 
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switching buck regulator employing multiple control loops. 
The feedforward loop will be designed to cancel the detri- 
mental effect of input filter interaction and thus improve 
the regulator performance and stability. The use of the 
feedforward control loop thus removes some of the conflict- 
ing design constraints mentioned above and makes an optimal 
input filter design more easily attainable. 

2.2.4 Objectives of Feedforward Loop Design 

The proposed feedforward loop will be designed such that : 

1. It will eliminate input filter interaction with the 
regulator loop. This will result in improvement in 
stability margins, audio susceptibility, output impe- 
dance and transient response of the regulator. 

2. It will allow the input filter to be optimized. Some 
of the design constraints that make an optimal filter 
design difficult to attain are removed with the addi- 
tion of feedforward. 

3. It will eliminate equipment interaction. Figure 12 
shows a switching regulator and its preregulator 
which may be a rectifier and a filter. The dynamic 
output impedance of the preregulator will interact 
with the switching regulator and may cause problems 
like loop instability, degradation of audiosuscepti- 


bility, output impedance and transient response. The 
addition of a feedforward loop will eliminate such 
interaction, thus isolating the switching regulator 
from equipment upstream. 






Chapter III 


CONCEPT OF POLE-ZERO CANCELLATION USED IN TOTAL 

STATE CONTROL 


3.1 CONCEPT OF POLS- ZERO CANCELLATION . 

The concept o£ pole-zero cancellation that was developed 
earlier [2,7,8] was implemented by feedback loops that sense 
the regulator output filter state variables and process this 
information to achieve better performance, and a control 
adaptive to filter parameter and load changes. 

Figure 13 [2 ] , shows a two- loop controlled switching 
buck regulator. The dc loop senses the converter output vol- 
tage and compares it with the referance voltage to generate 
a dc error signal for voltage regulation. The ac loop senses 
the ac voltage across the output filter inductor to generate 
an ac signal. Both ac and dc signals are processed through 
an operational amplifier summing junction to provide a total 
error signal at the output of the operational amplifier 
integrator. It is apparent that the error signal at the out- 
put of the integrator contains information regarding the 
output filter state variables — the inductor current and 
the capacitor voltage. 
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Figure 13: 


Two loop controlled switching buck regulator 
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It was shown, [2], that the feedback control loops when 
properly designed can provide complex zeros to cancel com- 
pletely the compless poles presented by the low-pass output 
filter of the power stage; It was also shown that the feed- 
back control loop has the ability to sense filter parameter 
changes and automatically provide pole-zero cancellation. To 
examine the adaptive nature of the control loops the open 
loop regulator transfer function G*£(s) is used 


G t (s) 


KZ (ioi) 
SP(j£U) 


(3-1) 


where K is a constant determined by the power stage and con- 
trol loop parameters and 

2 2 

Z C jti>) 5 1 + n i (3-2) 

P(joi) = 1 + “/ tll n2 ~ (3-3) 

P(j w) has complex poles corresponding to the output filter 
and 2(j(o) has complex zeros produced by the two loop feed- 
back control. 


hi 

/ -SL 
/ LC 

(3-4) 

n2 = 

1 

/LC 

(3-5) 


“nl 

(3-6) 


5 1 “ 2 T z 


- !S2(-L. + 
2 u 




R C C + R^C) 


(3-7) 


(3-9) 


r z = («, ♦ R S )C 2 + 5 ^ 


Ry - Jj^R^/R^) + R^j 


(R^ + R2 ) 


R, 


(3-10) 


L, C, and R c form the output filter as in Figure 13 . The 
control parameters can be chosen such that 

(3-11) 


u nl " w n 2 


C 1 ” ? 2 


(3-12) 


thus resulting in 


P(ju) = Z(jw) 


(3-13) 


ana 


s< s > = I 


(3-14) 

The open loop transfer function is of first order and is 
completely independent of output filter parameters. The 
adaptive nature of the control loop is apparent from the 
fact that the complex zeros imitate the change in the com- 
plex poles due to component tolerance, aging or temperature 
variations , thus preserving the pole-zero cancellation. 
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3.2 TOTAL STATE CONTROL 

The concept of pole-aero cancellation of the output filter 
characteristics led to the idea of using similar means to 
control the effect of peaking of the output impedance of the 
input filter. The objective of the work reported in this 
dissertation is thus to develop a feedforward loop that 
senses the input filter state variables and uses the infor- 
mation contained therein to eliminate the interaction bet- 
ween the input filter and the regulator control loop. Such a 
feedforward loop working in conjunction with existing feed- 
back loops forms a total state control scheme, which is il- 
lustrated in Figure 14 . The feedforward loop senses the in- 
put filter state variables and feeds this information to the 
error processor. 

The other inputs to the error processor are the ac vol- 
tage across the output filter inductor and the output vol- 
tage - as shown in section 3.1 these contain information 
regarding the output filter state variables. 

The pulse modulator thus has as its input information re- 
garding the state variables of the output filter and also 
the input filter. The duty cycle signal d(t), which controls 
the switch in the power stage, is thus also affected by the 
input filter state variables. It is shown later in this dis- 
sertation, in Chapters 4 and 5, that the feedforward loop 
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Figure 14: Total State Control. 
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can be designed such that interaction between the input fil- 
ter and the regulator is eliminated. 


Chapter IV 

MODELING OF THE POWER STAGE WITH INPUT FILTER 


The first step in the design and analysis of the feed* 
forward loop is to develop the small signal model of the 
power stage with input filter for the buck-boost, buck and 
boost type of switching regulators, using the averaging 
technique, [II. The modeling is carried out in the continu- 
ous conduction operating mode, in which the inductor current 
is always nonzero. This mode is the prevalent operating mode 
for most de-dc converters. The discontinuous conduction op- 
erating mode in which the inductor current is zero for some 
time during the cycle occurs at light loads and is seldom 
used as the intended design at full load. 

The modeling is carried out in the following steps — 

1. State space equation formulation during T ow - and T QFF . 

2 . State space averaging and perturbation. 

3 . Linearization and derivation of the small signal 
equations and the small signal equivalent circuit and 
state space model. 
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4,1 BUCK-BOOST CONVERTER SMALL SIGNAL MODEL DERIVATION 
The buck-boost converter is shown in Figure 15 . In the 

buck-boost converter shown the input filter is composed of 
R lj/ LI, R C1 and Cl. The load Is represented by R ^ while Vj 
and V 0 are the input and output voltages respectively. Dur- 
ing T ow , the switch S is on and the circuit as shown in Fig- 
ure 16 (a) . 

The equations describing the circuit are 


dt 


dt 


T — § = flux in core 

L P 


Ta ( -«u ■■ 


R C1 J 


R C1 A Ici ^1 
LI L p <p “ LI + LI 


Np N p hi 1^ 


^ R C1 + 


(4-1) 

(4-2) 

(4-3) 


dv 

Cl 

. ^ 

dt 

Cl " c 


“ v c 

dt 


v o 

Ve ■ 


Clip 


(4-4) 


(4-5) 


(4-6) 
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During T OF p the switch S is off and the circuit is as shown 
in Figure 16 (b). 

The equations describing the circuit are 


"s 

(4-7) 

TT ■ ¥ f-«ti - R ci> - TT * ET 

(4-8) 

dtfi ~^ R S R C + R S®L + R cV* \ v C 

dt L S (E 0 + Sj) H s CR c * 

(4-9) 

dv ci _ h.i 

Cl Cl 

(4-10) 

dv c NgR.^ v Q 

d t a s (R c + Rj_) C (R c + Ej) 

(4-11) 

... Vc • . ¥l¥ 

0 R C + R L V R C + V 

(4-12) 


The following vectors are defined 


u = j> r ] 
Z “ fv Q ] 


(4-13) 


resulting in the following state space equations 
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t on 

T 

1 0 FF 

x = x + B, u 

x a Aj X + B 

2 * c l 5 

Z= c 2 £ 


(4-14) 


where 


A 1 


_ R Cl N p ~ 


LI 


-Cl 

®P 

Cl 


L 1 L t 


-*p 

CILp 


*< R C1 +R P> 1 


K. 


P 

0 


0 

-1 


c <*l + V 


(4-15) 


(R L1 +R C1 ) 


B 1^II 


LI 

0 

1 

Cl 


B 2 “ B 1 


0 n 


-R- 


Tj S 


0 

% 


CLsCR^V 


0 ]■ 


0 

-*L 


N g (R c +V 

0 

-1 

C(R C +R L ) 


(4-16) 


(4-17) 

(4-18) 


c 1 « 


*L 


H C +B L 


(4-19) 
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C 


2 


Ws n _2 k_ " 
V R C + V R c+Rl 


(4-20) 


% 


Vc+Vl+Vl 

r c +k l 


(4-21) 


The state space averaged model over the entire period T is 


where 


i: = [dA^ + d'A23 J£ + [dB^ + d'Bj] u 
X - [dC 1 +d , C 2 ] x 

d -- duty cycle ratio = T on / ( T on +T qF5 ) 


(4-22) 


d' - 1-d (4-23) 

T => T +T 

ON * X 0FF 

The state space averaged model is perturbed thus - 
d a D + d 
d' = D* - d 

u = U + u (4-24) 

y = y + y 

x “ X_ + 

Assuming that the perturbation is small 

A A 

d x 

— << 1, — << 1 etc. leads to the following small 
D X 


signal linearized model 
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0 = [DA X + D r A 2 ] X + IDBj^ + D'B^ Vj 
V Q » [DC^ + D ! C 2 ] X 

x * [DA 1 + D»A 2 ] x + [DB 1 + D’B 2 ] Vj (4-25) 

+ [A x - A 2 ) X + (B x - B 2 Vj] d 

V 0 = tC l * C 2 3 - * + fDC l * D?C 2 3 i 

Defining 

A ** DA^ *f* D'A 2 

B = DB X + D r B 2 (4-26) 

C = BC 1 + D'C 2 

results in 
♦ 

x 5 A K + B Vj 

+ CA 1 - A 2 ) X + (B x - B 2 ) Vj] d ' (4"27) 

Vq * [C x - C 2 ] X d + C x 

Using Laplace transforms results in 

x(s) = [SI - ApS Vj(s) 

+ [SI - A3 -1 [(A 1 - A 2 ) X + (B 1 - B 2 ) V x 3 d(s) 
v Q (s) = [C x - C 2 ] X d(s) + C x(s) 


(4-28) 
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The state space model is derived from the two equations 
above, and is shown in Figure 17 . 

To derive the small signal equivalent circuit, it is 
first necessary to use the equation for v 0 to substitute for 
v c in terms of v Q , in the. small signal equations described 
above. Simplifying, the four equations are obtained 


R--NJD A di_ , 

V I " (R C1 * L^ <P + LI 

Np A 


(4-29) 


DN p . N p A 

Lp Lp 


C1 ~~sr “ - it * - if * d 


(4-30) 


c ^ T <P - g ♦ a 


*L 


(4-31) 


V C1 N S 

N_ 


VI fciW.- . BR ci N s 2 

dt N p *1,1 * Lp * 


(4-32) 


. R s N s T . DD 'Ws r t * 

“ ^ VHc + *l) Vc 


^l^^S 5 . R ci d( * N s 


HE 

r 


d + 


L-. * d V 0 


W 5(, - d ' ) V C1 N S 
l s (r c +e l> ' N P 








5 C1 I L1 M S 


R C1 IL S 


W CD “ D ) V C1 N S 




< R C + V * 


N, 


+ V, 


Eigure IS: Buck-boost converter power stage small signal 

equivalent circuit 
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Using a fictitious current i described by 

L s i = N s * (4-33) 

an equivalent circuit can be made up that is described by 
the four equations given above. This circuit will use the 
current i flowing through L g and it is thus the small signal 
equivalent circuit for the buck-boost converter/ as shown in 
Figure IS. 


4.2 BUCK CONVERTER SMALL SIGNAL MODEL DERIVATION 
The procedure used in deriving the small signal model for 
the buck converter is exactly similar to that used for th© 
buck-boost converter. The buck converter is shown in Figure 
19 . 

During T QN the switch S is on and the circuit is as 
shown in Figure 20 (a). 

The equations describing the circuit are 


dl Ll = ~ <R L1 * . fci V C1 V I 

dt LI TJ. LI \ " LI LI 

dl L R C1 . ^1 ^C1 R L V G 

dt =! L 1 L1~L 1 L L'' LO^+I^) 

dv ci _ *L1 h 

dt ~ Cl Cl 


Vl V C 

dt ~ C(R C + Rj.) ” C(R C + RJ 


(4-34) 


(4-35) 


(4-36) 


(4-37) 


50 


„ „ Vc 1 ! . Vc 

0 *L * R C \ + R C 


(4-38) 


During Tq ff , the switch S is off, and the circuit is shown 
in Figure 20 (b) . The equations describing the circuit are 


di Ll = ~ (R L1 * R C1 ) . V 
dt LI Hi " LI 



(4-39) 


JVt Vc 

dt L LCRj. + R c *) 

dv Cl JjLl 
dt * Cl 


(4-40) 


(4-41) 


dv c Vl v c 

dt ° + R c ) CO^ + R c ) 


VA Vc 
\ * R c \ * R c 


Ro = R„ + 


Wl 


2 ^ R c + R l 


A “ R C1 + \ * R c + R l 


r c*l 


The following vectors are defined 

JU = [Vj] X 

Z " E V 0 I 


Hi 

H 

V C1 


(4-42) 


(4-43) 


(4-44) 


(4-45) 


(4-46) 
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resulting in the following state space equation 


T 

ION 

x a X + B^u 

i" c li 


T 

^OFF 

i " *2 x + I^u 

Z = c 2 ^ 


where 


" CE L1 +S C1 ) 

*C1 

-1 

A 

LI 

Ll 

Ll 

u 

R C1 

- R 1 
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L 

L 

L 

«“c + V 

1 

-1 

0 

0 

Cl 

Cl 



rt 

r l 

n 

-1 

U 

C(R l + E c ) 

u 



(4-47) 


(4 -48) 


-<fu+fca l 
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x 

Cl 


^2 

L 
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zl 

LI 


0 
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C <VV 


0 
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LCR^y 

0 

-1 

c ( r l + V 


( 4 - 49 ) 


54 


B, 



0 
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The state space averaged model over the entire period T is 


g - [dAj + d'AgJs + [dB 1 + d' B 2 ]u (4-54) 

£ - [dCj + d'C s ]x 

where d = duty cycle ratio = T OH /(T OW +T opp ) 


d' - 1 - d 


T 


= T + T 

a on a opf 


( 4 - 55 ) 


The state space averaged model is perturbed and linearized 


in exactly the same way as for the buck-boost converter. The 
resulting small signal linearized model is 


0 * [DA 1 + D’A 2 ]X + [DB^ + D , B 2 JV i 
g = [DA 1 + D'A^x + [DB 1 + D'B 2 ]v 
+ E(A X - A 2 )X + <B 1 - B 2 )V T ]d 
V 0 = [DC X + D'C 2 ]X 
v Q = [C x - C 2 ]Xd + [DC 1 + D f C 2 ]x 


( 4 - 56 ) 
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Defining 

A « DAj + D f A 2 

B — DBj + D’B 2 (4-57) 

C = DCj + D'C 3 

results in 

x = Ax + Bv^ 

+ C <A 2 - A 2 )X + (B 1 - B^Vjjd (4-5S) 

V 0 = tC l - C 2 J ^ + 4 

Using Laplace transforms gives 
x(s) = [SI - A] "S VjCs) 

+ [SI - A)” 1 [(A 1 - A 2 )X + (B x - B 2 )V I ]d(3) (4-59) 

v Q (s) = [C x - CjIX 3(s) + C ji(s) 

The state space model is derived front the above two equa- 
tions and is shown in Figure 21 . 

The procedure for deriving the small signal equivalent 
circuit is exactly similar to the one used for the buck- 
boost converter. The four equations that result are 

A 

V I “ ^L1 * ^Cl^l ~ DR C1 ^ + L1 “dF" (4-60) 

+ V C1 " R C1 V 

Dv C i = (DR C1 - DR cl i^ 1 + v Q 

A 

+ L ~dT " ~ + v ci )a 


(4-61) 
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A 


Cl- 


dv ci 

dt 


*U “ Bi L - dI L 


(4-62) 


dv„ 


4 

T. R. 


dt ~ ~ (4-63) 

The small signal equivalent circuit is described by the four 
equations above, as in the buck-boost converter, and is 
shown in Figure 22 , 


4,3 BOOST CONVERTER SMALL SIGNAL MODEL DERIVATION 
The procedure used in deriving the small signal model for 
the boost converter is exactly similar to that used for the 
buck-boost converter. The boost converter is shown in Figure 
23 . 

During T 0JJ , t * le switcl1 S is on, and the resulting cir- 
cuit is shown in Figure 24 (a). 

The equations describing the circuit are 


dl Ll ~ <R L1 * * 0 * , ,fci, V C1 V I 
dt ~ U 11 LI T, ' LI LI 


dv 


Cl 

dt 


ci _ ci 


fL 

dt 


R, 


r hi- 


( R d + r £ ) 


dv. 


C 
dt 


-v. 


CCR C + y 


v„ = 


v c*l 


r c + *l 



(4-64) 

(4-65) 

(4-66) 

(4-67) 

(4-68) 


\ 


< > 












59 


During T 0 pp the switch S is off and the resulting 
shown in Figure 24 (b) . 

The equations describing the circuit are 

dl Ll ~ (R L1 * . ^C1 . V C1 V I 

dt LI 11 LI V U LI 

dv Cl *L1 h 

dt Cl Cl 

dL L ^ R ci _ h +R+ R A } 

dt L *L1 L (R C1 + R Z + K. c + 

*L V C . V C1 
“ L(R C + Bj) L 

dv c R^i^ v c 

"dT = c(r c h* y " c(R c + rT) 


Vi . , Vc 

■ S C + h h R c + h 


The following vectors are defined 


x » 


*11 

V C1 


u = [v x ] 

z ■ Iv 0 l 


resulting in the following state space equations. 


T 

ON 

x = A^x + Bj_u 
2 = C 1 £ 


T 

OFF 

x = A 2 x + B 2 u 
Z * C 2— 


circuit is 

(4-69) 

(4-70) 

(4-71) 

(4-72) 

(4-73) 

(4-74) 

(4-75) 





where 
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(4-82) 


The state space averaged model over the entire period T is 
x * [dA^+d'A^x + [dB 1 +d , B 2 ]u 

Z “ [dC^+d'C^: 


(4-83) 


IX 


where 


d 

d* 

T 


duty cycle ratio = T ON /(T OH +T OPP ) 

1 - d 

T + T 
ON OFF 


(4-84) 


The state space averaged model is perturbed and linearised 
in exactly the same way as for the buck-boost converter. The 
resulting small signal linearized model is 


0 = AX + BVj 

x = A x + B Vj 

+ [(A x -A 2 ) X H* (B 1 -B 2 )V I ]d 


V 0= cx 

v 0 - [C x -C 2 ] X d + C x 

where A = DAj + d'A 2 

B. = DB X + B’B 2 
C - DC X 4- D’C 2 

Using Laplace transforms gives 


<4-85) 


(4-86) 

(4-87) 


(4-88) 
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x(s) - [SI - A]"* 1 B v x (s) 

+ [SI-A]’ 1 [(A 1 -A 2 >X + (B 1 -B 2 )V I ja(s) 

v 0 (3) = [^-C^X d(s) -5* Cx(s) 


The state space model is derived from the above two equa- 
tions and is shown, in Figure 25 . 

The procedure for deriving the small signal equivalent 
circuit is exactly similar to the one used for the buck- 
boost converter. The four equations that result are — 


d± Li A „ 

L1 “dt" “ -(R L1 *5’Bci)i L i + " V C1 + V I 


(4-90) 


d±. R EL 

L ~dt 3 r ci^li ” (R ci V" dd 'r c +r l *L 


+ v ci - DV 0 + 


R^^dCD-D’) 


R„ + 




+ V 


dv 


Cl 


i dt “ “ h 


(4-91) 


(4-92) 


dt 


. '0 
D ^ V 


(4-93) 


The small signal equivalent circuit is described by the four 
equations above as in the buck- boost converter,- and is shown 
in Figure 26 . 


65 


The power stage small signal models developed include 
the input filter state variables, whereas earlier models 
[3,4/5,61 had treated the input filter only in terms of its 
output impedance and transfer function. The models developed 
in this chapter are used to analyze and design a feedforward 
loop that includes the input filter state variables. 










Chapter V 


IMPLEMENTATION OF THE FEEDFORWARD LOOP FOR A 

BUCK REGULATOR 

This chapter first presents an analysis that leads to a de- 
sign of the feedforward loops for a buck regulator. A small 
signal model that includes a general form of the feedforward 
loops is developed first. Analysis of the small signal model 
leads to a design of the feedforward loop. Implementation of 
the design is next discussed and two feedforward circuits 
are presented. The buck regulator alone is treated in this 
chapter, however the analysis and design procedure would be 
similar for the boost and the buck-boost regulators. 

S.l STATS SPACE MODEL OF BUCK REGULATOR 

The state space model of the buck regulator is shown in Fig- 
ure 27 . In the figure the feedforward loop har ae its in- 

A A 

put the input filter state variables i £1 and v ci (' input 
filter inductor current and capacitor voltage respectively) . 
These two inputs are multiplied by the transfer functions 
c 2 (s) and c 3 (s) whose properties are yet to be determined. 
The feedback loop has as its inputs the output voltage and 
the output filter inductor current. The feedback control in 
this case is the two loop control {the standardized control 
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module or SCM) developed earlier 12/7,8} and discussed in 
Chapter 3. The error processor in Figure 27 is thus com- 
posed of the blocks labelled c % , c 3 and the feedback, and 
has as its input information regarding the output filter and 
input filter state variables. The pulse modulator is repre- 
sented by its transfer function F M [2,8]. The rest of Figure 
27 is the state space model of the buck power stage devel- 
oped in Chapter 4. 

The feedforward and feedback signals are added and fed 
to the pulse modulator. In physical terms this means sens- 
ing the small signal variations in input filter inductor 
current and capacitor voltage, processing these variations 
(as represented by the blocks c 2 (s) and c 3 (s) in Figure 27) 
and adding the processed variations to the feedback signal. 

The total state feedforward/feedback error signal is then 

\ 

used to modulate the duty cycle of the switch for loop gain 
correction. 

The transfer function v Q ( s)/v 2 (s) , Figure 27, is used 
to design the feedforward because it expresses clearly what 
the feedforward does; also the resulting design is indepen- 
dent of the feedback loop parameters. The generalized small 
signal model for the buck regulator is developed next and 
used to write the transfer function v 0 (s)/v x (s) . 
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5.2 GENERALIZED SMALL SIGNAL MODEL OF BUCK REGULATOR , 
The generalized small signal model of the buck regulator is 
shown in Figure 23 . The regulator is modelled according to 
the three basic functional' blocks: power stage, error pro- 
cessor and duty cycle pulse modulator. The power stage model 
consists of two inputs: disturbances from the line $ and 

A 

the duty cycle control d, and four outputs: the output vol- 

A A 

tage v 0 , the output filter inductor current i L , the input 
filter capacitor voltage v ci and the input filter inductor 

A 

current i L1 . The error processor has as its input informs- 
tion regarding the output filter and the input filter state 
variables. The transfer functions F 3 , F AC and F DC constitute 
the two loop standardized control module (SCM) developed 
earlier [2,7,8], whereas the feedforward loop gains c s and 
c 3 are as yet unknown. The error processor processes infor- 
mation regarding the state variables of the input and output 
filters and feeds a total error signal to the pulse modula- 
tor, whose transfer function F^ was developed earlier 
[2,7,8]. 


The power stage transfer functions Fjj 
using the following equations: 


T 11 V I 

+ 

T 12 d 

= V 

o 

T 21 V I 



II 

h>-> 

T 31 V I 

+ 

T 32 S 

rH 

< > U 
II 

T 41 V I 

+ 

T 42 5 



etc. are written 


(5-1) 
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Thus it can be seen that 



and 



H P 
A 


11 



(5-2) 


(5-3) 


5.2.1 Development of Power Stage Transfer Functions 

As can be seen from equations (5-1) T tl , T 2l , T 31 and T 41 

A A 

can be evaluated with d = 0 and the other four with Vj = 0 . 
The starting point for the evaluation of the transfer func- 
tions is the small signal equivalent circuit model for the 
buck regulator power stage developed in Chapter 4, Figure 22 

5. 2. 1.1 Evaluation of T u , T 3l , T 3l and T 41 

* 

These transfer functions are evaluated with d = 0 in Figure 
22 . The resulting circuit is shown in Figure 29 . In Fig- 
ure 29 the input filter has been replaced by its forward 
transfer function H(s) and its output impedance Z(s). 

1 + sClR cl 

s 1 2 L1 Cl + sCKRj^ + R cl ) + 1 


H(S) 


(5-4) 


1 


f "• 
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S f 


2 (s) = 


S^L1C1R C1 + SCIR^R^ + SLl + R^ 
S 2 L1C1 + SCI (R l1 + R cl ) + 1 


(5-5) 


From the equivalent circuit of Figure 29 the following can 
be derived: 


1 I 

i ! 

i 


,) H 


T 11 


T 21 ~ * = 


'31 


T 


41 


A "* 

V I 

A 

K „ 

D(1 + sCR^) H 

A. 

v 

I 

A 

’cl 

a^H 

V I 

A 

A 

1 L1 

4 - a i H 


v. 




(5-6) 


where 


a l = 
A 
D 


S LCR^^ + sCR^{R a + R c + CR ^) + Rj^ 


* a x + D' z Cl + sCR^J 

V 

- duty cycle - ™“, V T = supply voltage 

I 


(5-7) 


V 


*L + R * 


“ \ 

Using equations (5-3) and (5-6) the following are derived: 


1 

11 

- DHj.d 

+ SCR C ) 

21 

= D(1 + 

sCR^) 

31 

= a l 


; 

A/B - 

■ a i 

41 

+ 

SLl 


In the derivation of equations (5-6) the resistance R C1 has 
been assumed negligibly small. This is not an unrealistic 
assumption since the ESR of the input filter capacitor (R C1 ) 
can be assumed negligibly small compared with the other re- 
sistances; also in the derivation of T 41 the following is 
used: 


A A 



(5-9) 
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5.2. 1.2 Evaluation of T 12 , T aa , T sa and T 4 a . 

These transfer functions are evaluated with Vj = 0 in Figure 
22 . The resulting circuit is shown in Figure 30 . 

From the equivalent circuit of Figure 30 the following 
are derived: 


T 


12 


- V o _ V o (R L ~ DZ) (1 * sCR c ) 


DA 


22 


32 


42 


^ V*L " D Z)tl + sCR l J 

d D V 

14 - ~ Z V a l + *L fl + SCR L )] 

a ' v 

*£l _ z V a i + V 1 + sCr l )] 

d ~ (r L 1 * sL1 i H i 1 A 


(5-10) 


Using equations (5-1) and (5-10) the following are der- 
ived: 



I 

d* 

>° 

l 

D 2 2) (1 

+ SCR C ) 

12 


D 



i 

>° 

l 

D 2 Z) (1 

+ sCR^) 

22 


dr l 



_ ' z V a l 

+ R^a 

+ SCR L ) ] 

32 


*L 



(5-11) 
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, 2 V a i * V 1 * SC V 3 

42 “ (Kjj + SLl) R^ 


where £>, and a^ are as defined in equation (5-7) and 



(5-12) 


5.2.2 Feedback Transfer Functions 

Transfer functions F 3 , F AC , F M and F DC constitute the feed- 
back. A two loop standardized control module (SCM) cont- 
rolled [2,7,8] buck regulator was used to obtain experimen- 
tal results that are discussed later in this dissertation. 
The buck regulator used is shown in Figure 31 , and with re- 
ference to that figure the following are defined [2,7,8] : 

F 3 = snL 


P AC ~ sCJR 


14 


DC 

R 


rrCi 


sC l R 14 


S- — > + -L.) 

* R X Z C 


x R ll //R 12 ' g “ R 


R x 


11 


Z c R 13 * SC. 


(5-13) 


2R 4 C| 

F m =* — (Pulse Modulator Transfer Function) 

M - Constant depending on the type of control used. 
5.3 DESIGN OF THE FEEDFORWARD LOOP . 

The transfer function Vq/Vj is used to design the feedfor- 

A 

ward. With Vj = 0 the following equations are derived from 
Figure 28 : 


+ C 2 V C1 + S 


(5—14) 


and 


dF 


12 _ 
A 

V A =5 
Cl 


i A 
X LI 


O 

dF 

~K 

dF 


32 


42 


(5-15) 


Substituting for v Cl , i L1 and d from (5-15) in (5-14) re- 
sults in 


v 


o _ 


^12 


V 


X 


1 - ° 2 (!§V m - °3<^) F H 


(5-16) 


C-'d- 
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In the absence of feedforward i.e, with c 2 
be seen that: 

Vj, _ V^L - p2z > + 3CB e> p M 

v x a l + D 2 Z(1 + sCRjJ 

The effect of peaking of the output impedance of the input 

2 

filter 3 is to cause a reduction in the term (R L -D 2) and 
also an increase in the denominator, thus resulting in a 
substantial loss of loop gain. 

With feedforward the detrimental effect of peaking of 2 
could be avoided by a proper choice of the feedforward loop 
gains c a and c 3 . Choosing 


“ c 3 = 0 it can 


(5-17) 


_ „ -D 

2 " V_F„ (5-18) 

O M 

c_ = 0 


leads to 



v oV*l - ° 2z) (1 - 1 - »<»c ) 

DA_ 

D 2 „ go[V V 1 + SCR L M 1 


(5-19) 


which can be simplified to 


i„ d2z)( i + sC W 

“ — “ — — — — 

v x - D Z) 

Cancellation of the two terms leads to 

v VPJ1 + sC&„ JR- 
o _ o M c T» 



(5-20) 


(5-21) 


Thus a proper choice of the feedforward loop gains has re- 
sulted in the transfer function v Q /v x being completely inde- 
pendent of the input filter output impedance 2. St is also 
noted from equations (5-11) and (5-16) that at frequencies 
other than the resonant frequencies at which 2 peaks, the 
gain of $ 3Ql is fairly small since 2 would be small at those 
frequencies. Thus the addition of feedforward would not af- 
fect, in any noticable manner, the open loop gain and phase 
margin at any frequency other than' those at which Z peaks. 

The following points regarding the feedforward loop de- 
sign can be made: 

1. It has been shown analytically that a proper choice 
of feedforward loop gains results in eliminating com- 
pletely the effect of peaking of 2 on the loop gain. 


S3 


2. The gain c 3 = 0, thus the inductor current informa- 
tion is not needed, only the input filter capacitor 
voltage information is used. 

3. The feedforward loop gains are independent of the in- 
put filter parameter values, and are free of any fre- 
quency dependent term. 

4. The feedforward loop gains are independent of the 
type of feedback control used. The pulse modulator 
transfer function F M is, however, an integral part of 
the design and thus the compensation depends on the 
type of duty cycle control used. The feedforward loop 
design process is independent of the particular type 
of control used and thus the same design can be used 
for other types of control, for example for single 
loop control, current injected control and others. 

S . 4 IMPLEMENTATION OF FEEDFORWARD 

The buck regulw >r used to obtain experimental results that 
are discussed later is shown in Figure 31 . The feedforward 
circuit processes the small signal variation across the in- 
put filter capacitor and adds this processed information to 
the feedback signal . 

Two circuit implementations of the feedforward design 
were used in making measurements and they are discussed 


next. 




LOAO 
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5.4,1 Nonadaptive Feedforward Circuit 

It was shown earlier that the feedforward loop gain is 

c 3 (s) = -D 2 /V 0 F m , equation (5-18) 

The nonadaptive feedforward circuit was developed for the 
buck regulator of Figure 31 . The key parameters of the re- 
gulator are as follows - 

Input-Output Parameters 

Vj — 25-40 volts Vq = 2Q volts P Q = 4Q watts 
Power Stage Parameters 

L = 230 micro H C = 300 micro F R^ = 0.2 ohm 
R c = 0.067 ohm 'nominal) R L = 20 ohm (load) 


Pulse Modulator Parameters 

M = V x T on = 0.88 * 10 3 V-sec 

Control Circuit Parameters 


e e = 

6.7 volts 

R u 

'= 33.3 Kohm 

R 13 

- 16.7 

Kohm 

13 “ 

2 Kohm 

R 14 

- 47 Kohm 

r 4 

-40.7 

Kohm 

n = 

0.65 

=; 

= 5600 picoF 

c a 

= 0.01 

microF 


The buck regulator was operated in a predetermined duty cy- 
cle control mode (constant VjT 0n control) , [2,7,8]. Substi- 
tuting for F , [2,8] and for D leads to 


2 V R 4 C i 


where M = VjT qk is constant. 

For the nonadaptive design the input voltage was kept cons- 
tant at Vj - 30 volts. Substituting in equation (5-22) it 
is calculated that <= a (s) = -0.03. The nonadaptive feedfor- 
ward circuit implementation is shown in Figure 32 . 'The in- 
put to the circuit is the input filter capacitor voltage and 
a series capacitor (27 microF) blocks out the dc component. 
The input is then multiplied by the gain of 0.03 implemented 
by the 5.1 Kohm and 164 ohm resistances. The feedforward 
signal available at the potential devider' network is then 
subtracted from the feedback signal available at the output 
of the integrator in the feedback loop. The result is then 
fed to the pulse modulator. The capacitor voltage fed into 
the operational amplifier subtracting circuit consists of 
two components - a small signal variation and a component 
corresponding to the switching frequency. The feedback sig- 
nal is also at the switching frequency, but the amplitude of 
the feedback signal is large compared to the switching fre- 
quency information in the capacitor voltage, and thus the 
second component has negligible effect. It is to be noted 
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that the circuit of Figure 32 constitutes the feedforward 
circuit and also the summing junction shown in Figure 31 . 

The nonadaptive circuit has the advantage of being ex- 
tremely simple and easy to implement. The gain of the poten- 
tial devider in the feedforward circuit is, however, a func- 
tion of supply voltage and thus the circuit of Figure 32 
cannot be used at any other value of supply voltage. Mea- 
surements made using this circuit are discussed in the next 
chapter. 


5.4.2 Adaptive Feedforward Circuit 

The adaptive feedforward circuit is shown in Figure 33 . 
From equation (5-22) it is clear that change's in supply vol- 
tage Vj will change the gain of the feedforward loop, the 
circuit of Figure 33 implements the feedforward of equation 
(5-22) and adjusts the gain automatically as Vj changes. 

The input voltage in Figure 33 is allowed to vary bet- 
ween 25v and 40v. Xt is fed to a voltage devider and then 
squared. The input filter ' capacitor voltage consists of a 
large dc component and this is blocked out by the 27 microF 
capacitor in series with the feedforward path. The small 
signal variation and the small magnitude component at 
switching frequency are then devided by the squared input 
voltage. A pair of resistances provides the final gain; the 





! 


feedforward signal now available is then added to the feed- 
back signal at the output of the integrater. As for the no- 
nadaptive circuit the swi telling frequency component in the 
input filter capacitor voltage is small compared to the cor- 
responding component in the feedback signal and thus its ef- 
fect on the duty cycle implementation is t7,egligible. 

The gain of the feedforward circuit is thus a function 
of input voltage Vj and is made adaptive to changes in Vj . 
The feedforward is thus capable of tracking ai„y variations 
in supply voltage. 



Chapter VI 


ANALYTICAL AND EXPERIMENTAL VERIFICATION OF THE 
FEEDFORWARD DESIGN 

This chapter presents extensive measurement data that 
was made co verify experimentally the feedforward design 
outlined earlier. The feedforward control was designed for a 
buck regulator and the same regulator is used in obtaining 
measurements. 

Measurement data pertaining to the open loop gain and 
phase margin are presented first; data made using the adap- 
tive feedforward circuit designed earlier confirm the adap- 
tive nature of the circuit. The audio susceptibility and out- 
put impedance measurements presented next show that the 
feedforward significantly improves performance in both cate- 
gories. Lastly, measurements of transient response are in- 
cluded that show that the feedforward improves the transient 
response. 
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MEASUREMENTS OF THE NONADAPTIVE FEEDFORWARD CONTROL 


The buck regulator with feedforward used to obtain experi- 
mental results is shown in Figure 34 , and is the same as 
that presented earlier in Chapter V. The parameters of the 
regulator are the same as given earlier in Chapter V with 
the single-stage input filter parameters specified as: 


R L1 = 0.2 ohm LI = 116 micro-H 


Cl = 20 micro-F 


The feedforward circuit used was the nonadaptive feedforward 
circuit for a fixed input voltage discussed in Chapter V, 
with Vj “ 30 V. The small signal open loop transfer function 
of the multi loop controlled buck regulator of Figure 34 
without feedforward can be expressed as [2/7,81 

G T< S) = F DcVlZ * + F 3 F AC F 22 F M t 6 " 1 ) 

In equation (6-1) E ^ and F*^ are 1: * ie power stage transfer 

functions, F DC , F^ and F AC are the feedback control loop 

transfer functions, as discussed in Chapter 5 (section 5.2). 

The peaking of the output impedance of the input filter, 

2(s), affects the transfer functions F ^ and as is ® een 

below-. 


F 1Z = 


F 22 = 


V q (R l - D*Z) (1 + sCR c ) 

D[ D z 2( 1 + sCR L ) + aj] 

V r l ~ 1 + sCr l> 

DR l [ T) Z Z( 1 + SCR ) + a ] 


(6-2) 


(6-3) 




LOAD 
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The peaking of 2(s) reduces the gain of both and I ?£ % and 
thus reduces the open loop gain at the resonant frequency of 
the input filter. 

The addition of the feedforward loop modifies the open 
loop transfer function as shown below — 


WmW 1 * sCR c> 

0 T “ " 


D a. 


* * C *L> 

D a. 


(6-4) 


The addition of feedforward modifies the transfer functions 
Fj^ and and thus it can be. seen from equation (6-4) that 
the open loop gain with feedforward is not affected by the 
peaking of the input filter output impedance 2(s). G' T (s) is 
now independent of 2(s) and is a function only of the feed- 
back loop parameters and the power stage parameters, unlike 
G T (s) of equation (6-1) which is affected by the peaking of 
Z(s). Equation (6-4) is also the open loop gain of the buck 
regulator without input filter, as can be seen from equa- 
tions (6-1), (6-2) and (6-3) by setting Z(s) = 0, and it can 
thus be concluded that the addition of the feedforward loop 
should eliminate completely the peaking effect of the input 
filter output impedance and that the open loop gain with 
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feedforward should he identical with the open loop gain 
without the input filter fl0,llj, 

SINGLE STAGE INPUT FILTER 

Measurements of the open loop gain and phase margin of 
the buck regulator of Figure 34 were made with and without 
feedforward, and are presented in Figure 35 (a) and (b). The 
input filter resonates at around 3 KHs and results in dis- 
turbances in the open loop gain and phase margin at that 
frequency. The feedforward eliminates all these undesirable 
disturbances as is evident in the figures, thus providing 
close agreement with theory. It can also be seen from Figure 
35 that the characteristics with feedforward are almost 
identical to the gain and phase margin plots of the buck re- 
gulator without input filter, thus providing close agreement 
with the analytical prediction made earlier. 

TWO- STAGE INPUT FILTER 

The experiment was further extended to the buck regula- 
tor with a two-stage input filter. The two-stage filter of 
Figure 36 was used with the following parameter values: 

R 1 = 0.2 ohm L x = 325 micro H 0^= 200 micro F 

R 2 = 0.02 ohm L 2 = 116 micro H C a = 20 micro F 

R 3 = 0.075 ohm (ESR of C x ) . 
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The feedforward circuit used was the same nonadaptive feed- 
forward circuit used to obtain the measurements of Figure 35 
with the feedforward input being the voltage at capacitor C 2 . 

The other parameters of the circuit are the same as used to 
obtain Figure 35 . Measurements of the open loop gain and 
phase margin were made with and without feedforward and the 
results are shown in Figure 37 (a) and (b). The open loop 
gain and phase are affected at the resonant frequencies of 
the two stages of the input filter because the output impe- 
dance of the input filter 2(s) peaks at both resonant fre- 
quencies. The use of the feedforward circuit eliminates the 
detrimental effects of the input filter, as is evident from 
Figure 37 . 

REMARKS 

The following points regarding the above measurements 
are noteworthy: 

(1) Measurements of the open loop gain and phase margin 
show that the input filter output impedance causes distur- 
bances in the gain and phase margin at the filter resonant 
frequencies and the addition of feedforward eliminates these 
disturbances, providing close agreement with theory. 

(2) The feedforward compensation circuit is independent 
of the input filter parameter values. 



Figure 36: Two- stage input filter. 
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(3) The feedforward compensation scheme is independent of 
the input filter configuration- It was demonstrated above 
that the same feedforward compensation network is equally 
applicable to a single stage input filter and a two-stage 
input filter. These observations lead to a stronger conclu- 
sion that the feedforward can provide effective compensation 
for an unknown source impedance. For example, a preregula- 
tor which often has an unknown, dynamic output impedance can 
interact with a DC-DC converter downstream and result in 
system instability. The feedforward compensation scheme out- 
lined can be used to isolate the switching converter from 
the source thus preventing interaction between the switching 
converter and equipment upstream. 

6.2 ADAPTIVE FEEDFORWARD MEASUREMENTS . 

The buck regulator w?th feedforward used to obtain experi- 
mental results is shown in Figure 34 and the parameters of 
the circuit are as specified in section 6.1. The adaptive 
feedforward circuit for variable input voltage presented in 
Chapter 5 was used in obtaining measurements .The two-stage 
input filter of Figure 36 was used with the same parameter 
values as in section 6 . 1 with the feedforward input being 
the voltage at capacitor C a . The value of R 3 was changed 
to- 

R 3 = 0.2 ohm (ESR of C 2 plus external damping resistance) 


Measurements were obtained at four values of supply voltage 
using the same adaptive feedforward circuit in all cases - 
this was done to confirm the adaptive nature of the feedfor- 
ward circuit. 

A computer program was written to calculate the gain 
and phase margin of the open loop transfer function with and 
without two-stage input filter, at various input voltages. 
Equation (6-1) was used to calculate the gain and phase mar- 
gin; setting 2(s) = 0 in the equation gives the gain and 
phase margin without input filter. The following expression 
for the two- stage input filter was used: 


z(s) = 


z \ * * sL a 


where 


1 + sC 2 ( Z x + R a ) +s‘L £ C £ 
3j(S) - (Rj * sti^ ) // (Rg ’ 1/sC^) 


(6-5) 


( 6 - 6 ) 


Figures 38 - 41 show the computed values of open loop gain 
and phase margin with and without the two-stage input filter 
for input voltages V x = 25v, 30v, 35v and 40v. It can be 
seen that the two- stage input filter resonates at two fre- 
quencies and at each frequency the output impedance Z(s) 
peaks, thus causing sharp fluctuations in the open loop gain 
and phase margin. Measurements of the open loop gain and 
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phase margin at each of the above values of supply voltage 
were obtained with and without feedforward and are also 
plotted on the figures. It can be seen clearly that the ad- 
dition of feedforward removes the sharp fluctuations in open 
loop gain and phase margin caused by the input filter, pro- 
viding close agreement with theoretical prediction made ear- 
lier. The analytical prediction that the open loop gain and 
phase margin with feedforward are identical to the charac- 
teristics without input filter is also confirmed, as exami- 
nation Figures 38-41 show. 

The two stage input filter was modified so that 
R 5 = 0.075 ohm (ESR of Cg) and measurements of the open loop 
gain and phase margin with and without feedforward were 
made. Figure 42 shows the calculated values of open loop 
gain and phase margin together with the measured values at 
V r = 25v. With the external damping resistance set to zero 
the effect of the input filter is seen to be more pro- 
nounced. Measurements without feedforward shown plotted on 
the figure also show the pronounced effect of the input fil- 
ter. The. addition of feedforward effectively eliminates the 
sharp fluctuations in gain and phase margin caused by the 
undamped input filter. 
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Figure 38: (a) Open loop gain at Vt^Sv: Calculated 

values and measured values (A) without 
feedforward 
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(b) Open loop phase margin at 

Vj-2Sv: Calculated values and measured values 

(A) without feedforward 


Figure 38; 
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input filter 


calculated values with input 
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(A) measured values with input 
filter and with feedforward 
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(d) Open loop phase margin at 

Vjs25v: Calculated values and measured values(A) 

with feedforward 
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calculated values with input 
filter, without feedforward 


(a) measured values with input 
filter, without feedforward 
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(a) Open loop gain at Vj=30v: Calculated 
values and measured values (a) without 
feedforward 
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(a) measured values with input 
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Figure 39: (c) Open loop gain at Vj =30v: Calculated 

values and measured values (a) with feedforward 



calculated values with input 
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(A) measured values with input 
filter and with feedforward 
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Figure 39: (d) Open loop phase margin at 

Vi=30v: Calculated values and measured values 
(a) with feedforward 
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(A) measured values with input 
filter, without feedforward 
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Figure 40: 


(a) Open loop gain at Vj=35v: Calculated 
values and measured values (a) without 
feedforward 





FREQ (HZ) 


Figure 40: (b) Open loop phase margin at 

Vj^Sv: Calculated values and measured values 
t&) without feedforward 
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Figure 40: 


(c) Open loop 4 gain at Vj=35v: Calculated 
values and measured values (zs>) with feedforward 
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(d) Open loop phase margin ah 

Vj-35v: Calculated values and measured values 

(A) with feedforward 


Figure 40: 
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(a) Open loop gain at Vj=40v: Calculated 
values and measured values (A.) without 
feedforward 
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(A) measured values with input 
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Figure 41: (c) Open loop gain at Vj=40v: Calculated 

values and measured values (a) with feedforward 
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calculated values without 
input filter 


calculated values with input 
filter, without feedforward 


(A) measured values with input 
filter and with feedforward 
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(d) Open loop phase margin at 
V =40 v: Calculated values and measured values 
(a) with feedforward 
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The following observations regarding the measurements 
are made: 

(1) The peaking of the output impedance of the first 
stage is more pronounced and has a greater effect on the re- 
gulator than the peaking of the second stage. * 

(2) Both- analysis and measurement results indicate that 
the open loop gain is higher at lower values of duty cycle 
D. This is further manifested by examining eqation (6-4). It 
shows that a lower value of D results in a higher gain. 

(3) The effect of input filter peaking varies with the 

input voltage. This is explained by noting that both and 

F ZZ' equations (6-2) and (6-3) depend on the duty cycle D. 

The effect of peaking of Z(s) is to cause a reduction in the 
z 

term (R^ -D 3) and the amount of reduction would be greater 
if D is larger. Consequently it is expected that at small 
values of Vj , when D is larger / the effect of peaking would 
be more pronounced. This is confirmed by examination of Fig- 
ures 38-41 . 

(4) The addition of the ' feedforward loop effectively eli- 
minates the perturbation in open loop gain and phase caused 
by the input filter. The feedforward works effectively at 
all four values of supply voltage Vj . Since the same feed- 
forward circuit was used in making all the measurements of 
Figures 38 - 42 the adaptive nature of the feedforward cir- 
cuit is confirmed. 


GfilN(DB) 


120 


ORIGINAL PASS t® 
OF POOR QUALITY 



FREQ (HZ) 


Figure 42: (a) Open loop gain at Vj=25v; Calculated 

values and measured values *a) without 
feedforward 
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calculated values with input 
filter, without feedforward 


(A ) measured values with input 
filter, without feedforward 
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Figure 42: (b) Open loop phase margin at 

Vj-25v: Calculated values and measured values 
(<&) without feedforward 
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calculated values with input 
filter, without feedforward 


(A) measured values with, input 
filter and with feedforward 



Figure 42; (c) Open loop gain at V 2 =25v: Calculated 

values and measured values (A) with feedforward 



calculated values without 
input filter 


calculated values with input 
filter, without feedforward 




(A) measured values with input & 
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Figure 42: (d) Open loop phase margin at 

Vj =25 v: Calculated values and measured values 
( a ) with feedforward 
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(5) From the results presented in section 6.1 and from 
measurements presented in this section it is logical to con- 
clude that the adaptive feedforward circuit can provide ef- 
fective compensation for an arbitrary, unknown source impe- 
dance for a variable supply voltage. The adaptive 
feedforward -circuit has been shown to be able to track the 
supply voltage and adjust the gain in accordance with supply 
voltage changes. The adaptive feedforward circuit can effec- 
tively isolate the switching regulator from its source impe- 
dance, thus preventing any interaction between the switching 
converter and equipment upstream. 

6.3 MEASUREMENTS OF CLOSED LOOP INPUT-TO-OUTPUT TRANSFER 
FUNCTION ( AUDIOSUSCEPTIBILITY ) . 

The closed loop input- to-output transfer function (audiosus- 
ceptibility) of a switching regulator is an important char- 
acteristic. It refers to the regulator's ability in attenu- 
ating small signal sinusoidal disturbances propagating from 
the regulator input to its output. The gain of the closed 
loop input- to -output transfer function should be as small as 
possible; thus the regulator will effectively attenuate 
noise at the input so as not to affect operation of the re- 
gulator payloads. Unfortunately, as pointed out in Chapter 
2, the peaking of the output impedance of the input filter 





and the peaking of the forward transfer function of the in- 
put filter increase the audiosusceptibility. Measurements 
of audiosusceptibility with and without feedforward were 
made using the two- stage input filter of Figure 36 with a 
damping resistor so that Rg = 0.2 ohm. Measurements were 
made by injecting a small sinusoidal signal at the input to 
the converter and then using a HP network analyser to mea- 
sure . the audiosusceptibility [3]. The feedforward circuit 
used was the adaptive feedforward circuit of Figure 33 and 
the buck regulator used is shorn in Figure 34 , with its 
parameters as specified in section 6.1. 

Figures 43 - 46 show the measured values of audiosus- 
ceptibility with and without feedforward at four values of 
supply voltage Vj ~ 25v, 30v, 35v and 40v using the same 

feedforward circuit in all cases. The top trace in each fig- 
ure is the plot without feedforward and it can be seen 
clearly that the audiosusceptibility is degraded 
at the two resonant frequencies where the outp- 
ut impedance Z(s) and the transfer function H(s) of the in- 
put filter peak , with the first stage resonating around 
600 Hz and the second stage around 3 KHz. It is also evident 
from the figures that the audiosusceptibility is dependent 
on duty cycle D or the supply voltage. At higher values of 
supply voltage when the duty cycle is low the audiosuscepti- 
bility is lower , specially at the lower frequencies. 
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Figure 45: Audiosusceptibility with and without feedforward 

at Vj^SSv (measurements) 
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The addition of feedforward substantially improves the 
audiosusceptibility, specially at the lower frequencies, as 
is evident upon examination of Figures 43 - 46 . The peaking 
effect of audiosusceptibility with feedforward loop is how- 
ever, more pronounced at the two resonant frequencies of the 
two-stage input filter; this is explained by noting that in 
equation (2-6), Chapter 2, the audiosusceptibility is shown 
to be affected both by the peaking of Z(s) and by the peak- 
ing of the forward transfer function H(s) of the input fil- 
ter. The peaking of H(s) cannot be controlled in any fashion 
by the addition of a feedforward loop since the control loop 
is not affected by H(s); thus the feedforward loop is effec- 
tive in cancelling 2(s) while the peaking effect of H(s) is 
manifested in audiosusceptibility. Figures 4? and 48 show 
the transfer functions H(s) and Z(s) of the two-stage input 
filter used. The peaking of H(s) at the two resonant fre- 
quencies is clearly seen. 

The two- stage input filter was modified by removing the 
damping resistor so that = 0.075 ohm. Measurement data of 
the audiosusceptibility with and without feedforward are 
shown in Figures 49 - 52 . The top trace in all these fig- 
ures is the closed loop input- to- output transfer function 
without feedforward and it can be seen that the gain is 
higher than in the earlier case (with damping resistance). 


Figures 53 and 54 show the forward transfer function H(s) 
and the output impedance 2(s) of the input filter used. Com- 
paring Figures 47 and 48 with Figures 53 and 54 it can be 
seen clearly that both Z(s) and H(s) peak at significantly 
higher values when the external damping resistance is re- 
moved. 

The addition of feedforward substantially improves the 
audiosusceptibility as is evident from Figures 49 - 52 . The 
audiosusceptibility with feedforward peaks at the two reso- 
nant frequencies of the two- stage input filter as before, 
but in this case the peaks are higher. This is explained b.y 
noting that H(s) of the two-stage input filter without damp- 
ing resistance peaks at a significantly higher value, as 
shown in Figures 53 and 54 , than that with a damping resis- 
tance as shown in Figures 47 and 48 ; thus the effect of 
H{s) on the closed loop gain would be expected to be greater 
in the former case. 

It can therefore be concluded that the addition of 
feedforward significantly improves the audiosusceptibility, 
specially at the lower frequencies. The audiosusceptibility 
with feedforward is affected by the peaking of H(s) since 
the effect of peaking cannot be eliminated via any control 
means. The same adaptive feedforward circuit was used in 
making all the closed loop gain measurements mentioned 
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The closed- loop output impedance of a switching regulator 
should be as low as possible in order that the regulator be- 
have as much as an ideal voltage source as possible. Howev- 
er, as pointed out in Chapter 2, the peaking of the output 
impedance of the input filter increases the closed- loop out- 
put impedance of the regulator. 

Measurements of the regulator output impedance with and 
without feedforward were made, using the two -stage input 
filter of Figure 36 with Rg = 0.075 ohm (ESR of C a ) . Mea- 
surements were made by injecting a small signal sinusoidal 
disturbance at the output in parallel with the load of the 
regulator, and then using a HP network analyser to measure 
the corresponding voltage and current [8,9]. The feedforward 
circuit used was the adaptive feedforward circuit of Figure 
33 . Figures 55 - 58 show the measured values of output 
impedance with and without feedforward at four values of 
supply voltage Vj = 25v, 30v, 35v and 40v using the same 
feedforward circuit in all cases. 

It can be seen clearly from the figures that the output 
impedance is increased at the two resonant frequencies of 










OCDB) 


143 


0F p °or quality 



FREQ ( HZ ) 


Figure 57: output impedance with and without feedforward 

(X) at Vj-35v 
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figure 38 : Output impedance with and without feedforward 

(X) at Vj ”40v 
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the two- stage input filter. This is a consequence of the 
disturbances in the loop gain produced at those frequencies 
by the peaking of the input filter output impedance Z(s). As 
was seen earlier in sections 6,1 and 6.2 the effect of Z(s) 
on the open loop gain depends on the duty cycle D, at higher 
values of D i.e, lower supply voltages, the effect of 2(s) 
on the open loop gain is higher. Thus the effect of Z(s) on 
the output impedance would be greater at lower supply vol- 
tages, and this is experimentally confirmed as can be seen 
from Figures 55 - 58 . 

The addition of feedforward almost totally eliminates 
the undesirable perturbations in the output impedance char- 
acteristic at all supply voltages. 

6.5 MEASUREMENT OF TRANSIENT RESPONSE AND STARTING OF THE 
REGULATOR . 

In this section measurements of output voltage and other 
parameters for a step change in input voltage or load are 
presented with and without feedforward. 

6.5.1 Small Amplitude Transient Response Measurements . 
Photographs of the output voltage ripple and other parame- 
ters are presented with and without feedforward for two cas- 
es: 

(1) Step change in supply voltage. 
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(2) Step change in load. 

The step changes mentioned above are small enough so that 
the output voltage remains in regulation throughout the 
transient response period. 

6. 5. 1.1 Step Change in Input Voltage. 

The buck regulator used is shown in Figure 34 with the par- 
ameters as specified in section 6.1. The input filter used 
was a single-stage input filter with the following parame- 
ters: 

R u = 0*2 ohm LI = 325 micro H Cl = 220 micro F 
The adaptive feedforward circuit of Figure 33 was used in 
making the measurements. 

The input voltage was abruptly switched from Vj = 30v 
to 40v and photographs of output voltage ac ripple , input 
filter capacitor voltage,, output filter inductor current and 
control voltage without using a feedforward loop were taken 
as shown in Figures 59 (a) and (b) , 60 (a) and (b), respec- 
tively. The control voltage is the input to the pulse modu- 
lator; without feedforward it is the output at the integra- 
tor in the feedback loop while with feedforward it is the 
sum of the above signal and the feedforward signal. Figures 
61 (a) and (b) , 62 (a) and (b) show the photographs of the 
same variables with feedforward control. Comparing, for ex- 
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ample, the photographs of the output voltage ripple with and 
without feedforward. Figures 59 (a) and 61 (a), it can be 
clearly seen that the transient response is improved with 
the addition of feedforward. The amount of overshoot is 
less with feedforward. The magnitude of the oscillations in 
the output voltage caused by the interaction between the in- 
put filter and the regulator control are also lessened with 
the addition of feedforward. 

Comparison of the photographs of input filter capacitor 
voltage, output filter iivhistor current and control voltage 
do not reveal much difference between the two cases (with 
and without feedforward) . This may be explained by noting 
that the gain of the feedforward loop is fairly small ( 0.03 
for Vj = 30v ) and thus the feedforward signal would be 
fairly small in amplitude. The addition of such a small am- 
plitude signal to the fairly large amplitude waveforms re- 
corded on the photographs will not show very clearly. The 
output voltage ripple, however, is small in magnitude and 
the effect of adding feedforward shows clearly. A computer 
program was written to simulate the step change in voltage; 
results from the program are presented in the next chapter 
and show close agreement with the measurements. 

Thus it is concluded that feedforward improves tran- 
sient response for the case where the supply voltage is sub- 
jected to a step change. 
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Output voltage ripple (a) and input filter 
capacitor voltage (b) without feedforward 


( a ) Y — O.lv/div 

(b) Y — 5v/div 


X- — 1 msec/div 
X — 1 msec/div 
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(b) 

Figure 60: Output filter inductor current (a) and control 

voltage (b) without feedforward 

• (a) Y — 0.5 A/div X — 1 msec/div 

(b) Y — 0.5v/div X — 1 msec/div 
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Figure 61: Output voltage ripple (a) and input filter 

capacitor voltage (b) with feedforward 

(a) Y — O.lv/div X — 1 msec/div 

(b) Y — 5v/div X — 1 msec/div 
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6.5, 1.2 Step Change in Load. 

The buck regulator used is shown in Figure 34 # with the 
parameters as specified in section 6.1 with the following 
changes : 

Vt = 25v ' R l4 = 3.69 Kilo ohm 

A single stage input filter was used with the following par- 
ameters : 

R L1 = 0.2 ohm LI = 325 micro H Cl = 100 micro F 
The adaptive feedforward circuit of Figure 33 was used in 
making the measurements. 

The load was switched repetitively between =10 ohms 
and R l = 20 ohms using a transistor switch. Figures 63 (a) 
and (b) show the photographs of the output voltage ripple 
without and with feedforward# respectively# as the load is 
switched. The output voltage ripple without feedforward. 
Figure 63 # shows distinct oscillations caused by the inter- 
action between the input filter and the regulator control 
loop. The oscillation frequency coincides with the input 
filter resonant frequency. These oscillations are eliminated 
with the addition of feedforward, as is evident from Figure 
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6*5.2 Large Amplitude Transient Response Measurements . 
Photographs of the output voltage ripple and other parame- 
ters are presented with and without feedforward for two cas- 
es: 

(1) Large step change in supply voltage. 

(2) Starting of the converter. 

The cases mentioned above are large signal changes so that 
the voltage regulation is momentarily lost for part of the 
transient response period. 

6.5.2. 1 Large Step Change in Supply Voltage. 

The buck regulator used is shown in Figure 34 with the par- 
ameters as specified in section 6.1. The input filter used 
was the single-stage input filter of section 6. 5. 1.1. The 
adaptive feedforward circuit of Figure 33 was used in making 
the measurements with feedforward. 

The input voltage was abruptly switched from Vj = 40v 
to 25v and photographs of the output voltage were made with 
and without feedforward; Figures 64 (a) and (b) show the 
photographs . 

The step change in supply voltage is large enough to 
cause the regulator to lose regulation — the output voltage 
drops by about 1.5v before the regulator recovers. This may 
be explained by noting that for such a large change in sup- 
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ply voltage the input filter capacitor voltage drops down 
close to 20v , and since this value is lower than the design 
range of 25v - 40v for the regulator supply, the ‘regulator 
loses regulation. Further details are given in the next 
chapter which presents results from a computer program writ- 
ten to simulate this large step change. The results show 
close agreement with the measurements presented here. 

Figures 64 (a) and (b) show that the behavior of the 
regulator is similar with and without feedforward for such a 
large step change in supply voltage. This is expected since 
the regulator control loop momentarily loses its control 
function during this transient period. Since the feedforward 
is designed to compensate for the effects of input filter 
interaction via a duty cycle modulation scheme it is expect- 
ed that the feedforward does not contribute anything under 
these conditions. Computer based simulation results present- 
ed in the next chapter confirm the measurement results. It 
is also to be noted that the feedforward does not have any 
detrimental effect on the transient response. 

6. 5. 2. 2 Measurements of the Start Up Behavior of the 
Regulator 

This section investigates the start-up behavior of the 
switching regulator. The regulator used is shown in Figure 
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Figure 64: Output voltage without feedforward (a) and with 

feedforward (b) 


(a) 

Y — 0 . 5v/div 

X — 0 . 5 msec/di v 

(b) 

Y — 0 . 5v/div 

X — 0.5 msec/div 
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34 with the parameters as specified in section 6.1. The sin- 
gle stage input filter of section 6.5.1. 1 was used with the 
supply voltage set at 30v, The adaptive feedforward circuit 
of Figure 33 was used in making the measurements with feed- 
forward. 

Prior to starting, the output voltage and the output 
filter inductor current were both zero. Photographs of the 
output voltage and output filter inductor current were made 
with and without feedforward and are shown in Figures 65 (a) 
and (b) 66 (a) and (b) . The output voltage without feedfor- 
ward builds up from zero to 2Qv (regulated output) in about 
3 msec. The inductor current rises sharply at starting but 
is limited to about 6.0& by the peak current protection cir- 
cuit built in with the regulator. It settles down to its 
steady state value in about 4 msec. 

The photographs with feedforward show similar behavior 
- thus the feedforward does not contribute significantly to 
this transient period nor does it present any detrimental 


effects. 
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Figure 65: Output voltage (a) and output filter inductor 

current (b) without feedforward 

(a) Y — 5v/div X — 0.5 msec/di v 

(b) Y — 1 A/div X — 0.5 msec/di v 
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(a) 




(b) 

Figure 66: Output voltage (a) and output filter inductor 

current (b) with feedforward 

(a) Y — 5v/div X — 0.5 msec/div 

(b) Y — 1 A/div X — 0.5 msec/div 
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Cb) 

Figure 66: Output voltage (a) and output filter inductor 

current (b) with feedforward 

(a) Y — 5v/div X — 0.5 msec/di v 

(b) Y — 1 A/div X — 0.5 msec/div 







6.6 CONCLUSIONS 

Extensive measurements made to verify experimentally the 
feedforward design are presented in this chapter. Measure- 
ments of the open loop gain and phase margin, closed loop 
gain, output impedance and transient response confirm the 
effectiveness of the feedforward in improving performance, 
as predicted in the analysis. The following points are note- 
worthy: 

1. The feedforward eliminates the detrimental effect on 
open loop gain and phase margin of the output impe- 
dance of the input filter. 

2. The feedforward circuit is shown to be independent of 
input filter parameters and also independent of input 
filter configuration. 

3. The feedforward effectively eliminates the interac- 
tion between an unknown dynamic source impedance and 
the regulator control loop. 

4. The closed loop input- to- output transfer function 
(audiosusceptibility) and output impedance are both 
improved significantly by the addition of feedfor- 
ward. 

5. The addition of feedforward improves transient res- 
ponse in those cases where the interaction between 
the input filter and the control loop degrades, the 
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0 
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response. Examples of the above are small step chang- 
es in supply voltage and load, and results for these 
cases are included. For large signal transient behav- 
ior the feedforward does not in any way degrade the 
performance — examples of these are a large step 
change in supply voltage and the start up behavior of 
the regulator. 

No detrimental effects have been observed due to the 
use of the proposed feedforward compensation scheme 
through the course of study and through extensive ex- 
periments when the system was subjected to different 

r 

forms of small and large signal disturbances* 


Chapter VII 


DIGITAL SIMULATION OH BUCK REGULATOR FOR 
TRANSIENT ANALYSIS 

This chapter deals with a program developed to simulate 
the real time behaviour of the regulator with and without 
feedforward control under both steady state and transient 
conditions. Results for a step change in supply voltage are 
included that are in close agreement with the measurements 
presented in the last chapter. 

7.1 DESCRIPTION OF THE METHOD 

The starting point is the definition of the buck regu- 
lator in terms of state equations [12,13,14]. The switching 
regulator u^ed was the one used in Chapter 5 and 6 and is 
shown in Figure 34 The buck regulator without feedforward is 
described by a set of state equations. The six state varia- 
bles used are: 

x(l) - i L1 (current ininput filter inductor) 

x(2) « v Ci (input filter capacitor voltage) 

x(3) = e c (control voltage, input to pulse modulator) 

x(4) = i^ (output filter inductor current) 

x(5) = v_ (output filter capacitor voltage) 
c 

x(6) = (compensation loop capacitor voltage) 
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The control voltage is the output of the integrator in the 
feedback control loop if no feedforward is used and it is 
the sum of the integrator output and the feedforward signal 
if feedforward is used. The state variable- e_ is the voltage 
at the capacitor in the compensation loop. The pulse mo- 
dulator has as its input the control voltage e^ and it con- 
verts the voltage to a duty cycle signal d(t). The transfer 
function of the pulse modulator is a constant if the supply 
voltage is a constant-, [8], for the constant volt-sec. 
(VjT on ) control mode that was used. Thus it can be seen that 
the six state variables defined above completely describe 
the buck regulator system of Figure 34 . 

The states equations for the complete system without 
feedforward are developed [12,13,14] for the two time per- 
iods T ow and T OPF . 

During Tq^ the transistor switch SI is on while diode 
S2 is off and the system is characterised by the following 
state equations: 

x = FI x + Gl u (7-1) 


where 








In equation (7-6) \F ? is the supply voltage, E a is the refer- 
ence voltage, E Q is the r duration voltage drop of the power 
transistor and E D is the conduction voltage drop across the 
diode. 

During the transistor SI is off while diode S2 is 

conducting and the equations are: 

x = F2 x + G2 u (7-7) 

where 

F2 = FI except for 

F2 (2,4) a 0 

F2 (3,2) - o (7-6) 

F2 (4,2) = 0 


and 


G2 = 0 except for 
G2 (1, 1) = 1/Ll 


G2 (3, 2) 


C 1 R 14 


G2(3,4) = ■— 


( 7 -! 


G2 (4, 4) 



Since u is not a function of time the solution to the state 
equations (7-1) and (7-8) can be written as: 


During T 0Fp 

x(t 

+ « “ * p 

* £(t) + D P * u 

(7-10) 

During T 0M 

x(t 

+ T > - *1, 

* £(t) + D n * u 

(7-11) 


where are the state transition matrices defined for 

a small fixed step size T [12] and 


$ p 


* 


J e -F2 ‘ s ds 


G2 


°N = *N * 


T 


"FI • s . 
e ds 


G1 


(7-12) 


Equations (7-11) and (7-12) are used to simulate the behav- 
ior of the regulator without feedforward. Starting at the 
beginning of the ON time period equation (7-12) is used to 
propagate the state. The step size T is defined as half the 
ON and OFF times and this determies the matrices d , D. T/ 
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<j)p and D F . Substitution into equation (7-11) along with 
the initial values at the beginning of the ON period propa- 
gates the state through time T. Similar substitution propa- 
gates the state through the ON period £ 12, 13, 14]. For the 
type of duty cycle contol used, constant VjTq n control, the 
ON time is fixed whenever Vj is fixed. 

At the end of the ON period, equation (7-10) is used to 
propagate the state during the OFF period, with the initial 
condition being the state at the end of the ON interval. The 
end of the OFF interval is determined as the point when the 
control voltage e c equals the comparator voltage S^, which 
was 7.0v in this case. Newton's iteration is performed to 
find the exact point in time when e c equals E T £12,13,14], 

The equations describing the system with feedforward 
are very similar to those without feedforward. The nonadap- 
tive feedforward circuit of Chapter 5 was used in the simu- 
lation program and it is shown in Figure 67 . The output of 
the capacitor is included as a state variable. The state 
variables are thus defined as; 

x(l) = i LI 
x(2) - v C1 

x(3) - v f (output of C f ) 
x(4)= e c 
H(S)= i L 


(7-13) 
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x(6)= v c 
x{7)= e R 

The equations during T ON and T 0Wp are the sane as equations 
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(7-1) and (7-8) with the following definitions: 
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LI 
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C l**14 


0 o 

0 0 

0 0 



(7-15) 
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K3 = 


R f2 

R fi -5- R f2 


(7-16) 


In equation (7-15) kl, ki and u are as defined earlier. 


I? 2 ~ FI except for 
F2(2,5)= 0 
F2(3,5)= 0 
F2 ( 4 , 2 ) = 0 
F2(4,5)= -R c *kl/ C[ 
F2(5,2)= 0 

and 

G2= 0 except for 
G2 ( 1 , 1 ) = 1/LI 
02(4,2)= 1/C[*R J4 
G2 (4, 4)= n/Cj-R 4 
G2 ( 5 , 4 ) = -1/L 


( 7-f 17) 


(7-18) 


The solutions of the state equations and the procedure 
for propagation of the state from cycle to ccle is identical 
to the procedure for the case without feedforward. 

The discrete simulation technique outlined [12,13,14] 
lends i.tself easily to simulation of transient response for 
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a step change in supply voltage. The change in supply vol- 
tage is made to occurr at the start of a new ON period; thus 
all that is necessary to do is to change the value of ^ON 
before the start of a new cycle. 

7.2 DESCRIPTION OF THE PROGRAM 

A flowchart of the program [12,13/14] is shown in Figure 63 
while a listing of the programs is included in the Appendix. 
The flowchart is for the simulation without feedforward, the 
simulation with feedforward proceeds in an identical fash- 
ion. 

A minimum off time duty cycle control is implemented in 
the program. The regulator also has peak current protection 
and this capability is also programmed. During the ON period 
the current in the switch may rise above the set value of 
the peak current. Newton's iteration [12,13,14] is performed 
to find the exact point in time where this occurrs and the 
program will automatically terminate the ON time calculation 
when the current exceeds the set value. 



BEGIN 


Preset data, calculate 
Kl, K2. 

Initialize all matrices, 
vector to zero. 


Set up matrices PI, 
F2, Gl, G2. 

T OFF = T ON <V l“ E o } / E o 
T « FRAC * T™ 


Compute D p 



T = FRAC * T, 


Compute <j> N , D w 


Initialize x, output V c 
and u. IT = 0 


Flow chart for discrete 
simulation 
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Figure 63: continued 
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Figure 68: continued 
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7.3 SIMULATION RESULTS 

The simulation program was used to simulate the behavior of 
the system with and without feedforward for a small step 
change in supply voltage from 30v to 40v and also for a 
large step change from 40v to 25v. Simulation results pre- 
sented are in close agreement with the experimental results 
presented in Chapter 6. 

7.3.1 Simulation of steady state operation 
The computer program discussed can be used to simulate the 
steady state operation of the regulator system [12,13,14]. 
An accurate picture of the behavior of the system is ob- 
tained since the simulation method used is an exact method. 

The buck regulator of Figure 34 with feedforward was 
simulated with the following single-stage input filter par- 
ameters : 

Rkl= 0.2 ohm LI = 325 microET Cl =220 microF 
The supply voltage Vj was set equal to 30v and the other 
parameters were as presented in section 6.1, Figures 69 (a) 
- (g) show the plots of the output voltage, input filter in- 
ductor current, input filter capacitor voltage, output fil- 
ter inductor current, the feedforward voltage v f , the con- 
trol voltage and the voltage e r . It is noticed that the 
switching period is around 50 microsecs and that the output 
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voltage ripple is around 60 milivolts. The input filter in- 
ductor current has a small magnitude ripple component at the 
switching frequency, as does the input filter capacitor vol- 
tage. The voltage fed forward v f consists of the switching 
frequency component, however its magnitude is around 50 mi- 
livolts. The switching frequency component of the control 
voltage is seen to be around 800 milivolts and thus the ef- 
fect of adding the switching frequency component from the 
input filter capacitor voltage to the control voltage would 
be negligible as discussed in Chapter 5 ; since the magni- 
tude added would be around 1.5 milivolts (v f multiplied by 
the feedforward gain) . 

7.3.2 Simulation for a small step change in s upply voltage 
The same buck regulator was used with its parameters as spe- 
cified in section 7.3.1. A value of 0.6 ohm for was 

used. The supply voltage Vj was abruptly switched from Vj = 
30V to 40V. Experimental results for this step change are 
given in Chapter 6 where a value of 0.2 ohm was specified 
for R Li * A value of 0.6 ohm was used in the simulation as 
this represents the combined winding resistance and source 
impedance. The values of the other parameters used are: 

T 0 jpf>( minimum) = 5 microseconds 

I Q { maximum) = 6 amps, (set value of peak current) 

E q = 0.2 volt E d = 0.7 volt 


(7-19) 
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Figure 69 : 


(a) Output voltage during steady state 
operation. 
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Figure 69: 


(b) Input filter inductor current during 
steady state operation 
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Figure 69: (c) Input filter capacitor voltage during 

steady state operation 
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Figure 69: (d) Output filter inductor current during 

steady state operation 
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Figure 69: (£) Control voltace during steady state 

operation 
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Figures 70 (a) - (d) are plots of the computed values 
of output voltage, input filter capacitor voltage, output 
filter inductor current and the control voltage, all without 
feedforward. 

The nonadaptive feedforward design of Figure 67 was 
used in the simulation of transient response with feedfor- 
ward. The feedforward circuit parameters used at Vj = 30V 
were! 

C f = 27 microF R^j =5.1 Kohm R f2 - 164 ohm 
when the supply voltage is switched to 40V the value of R fa 
in the program is changed to 90.97 ohm at the switching ins- 
tant thus changing the feedforward circuit gain in accor- 
dance with the feedforward design. 

Figures 71 (a)-(d) are plots of the computed values of 
the same variables with feedforward. Experimental measure- 
ments for the same step change in supply voltage were pre- 
sented in Chapter 6 and are repeated here for convenience; 
Figures 72 (a)-(d) are the measurement waveforms without 

feedforward loop while Figures 73 (a)-(d) are for the case 
with feedforward loop. 

Comparing the plots of the computed values with and 
without feedforward it is clearly seen that the analytical 
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Figure 70: (b) Input filter capacitor voltage simulation 

without feedforward 
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Eigure 70: (c) Output filter inductor current simulation 

without feedforward 
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Control voltage simulation without 
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(c) Output filter inductor current simulation 
with feedforward 
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(b) 

Figure 72: Measurement without feedforward: output voltage 

(a) and input filter capacitor voltage (b) 

(a) y_o.lv/div x — 1 msec/div 

(b) Y 5v/di v x — 1 msec/div 
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Figure 72: Measurement without feedforward: output filter 

inductor current (c) and control voltage (d) 

(a) Y— 0.5 A/div X— 1 msec/div 

(b) Y 0 . 5v/div x— 1 msec/div 
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Measurement with feedforward: output voltage (a) 
and input filter capacitor voltage (b) 

( a ) Y — 0 . lv/div X — 1 msec/di v 


Figure 73 


X — 1 msec/div 


Y — 5v/div 



• * 


. ••••' 




Measurement with feedforward 
inductor current (c) and co: 

( a ) Y — 0.5 A/div X — 1 
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results agree fairly closely with the experimental data. The 
analytical results show that the feedforward is effective in 
reducing the output voltage oscillations — the amplitude of 
the oscillation is reduced (Figures 70 (a) and 71 {a))j this 
is confirmed by the experimental results (Figures 72 (a) and 
73 (a)). -The plots of the input filter capacitor voltage, 
output filter inductor current and the control voltage are 
also in excellent agreement with the experimental data. Also 
in excellent agreement with the measured data is the obser- 
vation that the plots of filter capacitor voltage, inductor 
current and control voltage do not show any noticable dif- 
ference between the two cases (with and without feedfor- 
ward) . Ah explanation for this observation was provided in 
Chapter 6. 

It can thus be concluded that the addition of feedfor- 
ward improves the transient response for the step change in 
Vj , as demonstrated both from the analytical result as well 
as experimental data. 

7.3.3 Simulation for a Large Step Change in Supply Voltage 
The regulator used was the one used in Section 7.3.1 with 
the same parameters. The single stage input filter parame- 
ters used were: 

R u = 1 ohm LI = 325 microH Cl = 220 microF 



The supply voltage was abruptly switched from Vj = 40V to 
25V. Experimental results for this step change are given in 
Chapter 6 where a value of 0.2 ohm was specified for R Ll . A 
value of 1 ohm was used in the simulation as this represents 
the combined winding resistance and source impedance. The 
other parameters are as specified in Section 7.3.1. 

Figures 74 (a-d) are plots of the computed values of 
output voltage, input filter capacitor voltage, output fil- 
ter inductor current and the control voltage, all without 
feedforward. 

The nonadaptive feedforward design of Figure 67 was 
used also in the simulation of transient response. When the 
voltage is switched from 40V to 25V the value of R fa is 
changed from 90.97 ohm to 230.27 ohm in order to change the 
gain of the feedforward circuit in accordance with the feed- 
forward design. Figures 75 (a-d) are plots of the computed 
values of the same variables with feedforward. Experimental 
measurements for the same step change in supply voltage were 
presented in Chapter 6 and are repeated here for conveni- 
ence; Figure 76 shows the measured output voltage ripple 
without and with feedforward. 



Comparing the plots of the output voltage ripple with 
the measured data it is clearly seen that the analytical re- 
sults agree closely with the experimental data. It can be 
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(a) Output voltage 
feedforward 
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Figure 74: (c) Output filter inducto 

without feedforward 
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Figure 75: (a) Output voltage simulation with feedforward 
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Figure 75: (d) Control voltage simulation with 

feedforward 
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(b) 

Figure 76: Measured output voltage ripple without (a) and 

with (b) feedforward 

(a) Y — 0.5v/div X — 0.5 msec/di v 

(b) Y — 0.5v/div X — 0.5 msec/di v 
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seen that the regulator loses regulation momentarily and the 
output voltage drops down to around 18.3V before regulation 
is regained. It should be noted when comparing the analytic 
cal plots without feedforward (Figures 74) with those with 
feedforward (Figures 75 ), the transient waveforms are almost 
identical. This is expected, since during this large signal 
transient the regulator has momentarily lost its voltage re- 
gulation. Therefore the ability of feedforward compensation 
via duty cycle modulation is also lost momentarily. The re- 
gulator is operating in an open loop fashion (under minimum 
off time control). This can be verified by observing the 
voltage waveform input to the regulator. 

The input filter capacitor voltage in Figures 73 (b) 
and 75 (b), which is also the input voltage to the regulator 
goes down close to 20V in both cases. For such a low input 
voltage to the regulator the duty cycle would have to be 
close to unity to maintain regulation, which is only possi- 
ble if the off time is decreased sharply. However a minimum 
"off t ^ ine 5 microseconds is implemented for various mag- 
netic reset purposes and thus the regulator loses regulation 
till the input filter capacitor voltage rises. 

The analytical plots of the control voltage without and 
with feedforward both show the large disturbance in control 
voltage as a result of the loss of regulation. The plot with 


feedforward does show a somewhat smaller peak value of the 
overshoot than the one without feedforward but this differ- 
ence is about 9%. Negligible difference is noticed between 
the analytical plots of the inductor current without and 
with feedforward. The large drop in current as a result of 
loss of regulation is noticed in both plots, as is the over- 
shoot when regulation is regained. 

The analytical results thus strongly support the exper- 
imental data and show that the regulator momentarily loses 
its control function. It is also confirmed through simula- 
tion and testing that the feedforward does not have any de- 
trimental effect under large transient conditions. 


Chapter VIII 


DISCRETE TIME DOMAIN STABILITY ANALYSIS OE BUCK 

REGULATOR 

An important concern is the interaction of an input filter 
with a switching regulator which often results in system 
instability. As described in detail in Chapter 2, the input 
filter interaction with the regulator control loop causes a 
reduction in the loop gain which may result in system inst- 
ability. When a switching regulator is acquired as a 'black 
box' for use in a system and an input filter is put in ex- 
ternally, the improper choice of filter parameters can cause 
the system to be unstable [4,5]. Work presented earlier 
[4,5] showed that the input filter design should be incorpo- 
rated in the regulator design itself in order to avoid loop 
instability and other input filter related problems. 

In this Chapter the stability of a buck regulator is 
examined by varying the input filter parameter values. Sys- 
tem instability is predicted analytically, and backed up 
with experimental data. It is then shown that the addition 
of the feedforx^ard control enables one to stabilize the sys- 
tem. This is, verified both analytically and experimentally. 
One can conclude that the addition of feedforward , in fact, 
isolates the regulator from the input filter; any input fil- 
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ter can be used with the buck regulator without causing loop 
instability. 

8.1 ANALYTICAL INVESTIGATION OF STABILITY 

An exact time domain analysis is performed to calculate the 
eigenvalues of the system with a single stage input filter; 
it was used to study how the eigenvalues move with a change 
in the input filter parameter (inductance). To correlate the 
’exact 1 discrete time domain analysis results with the ’ap- 
proximate’ continuous time average analysis presented in 
Chapters 2 throuyh 5, a computer program which calculates 
the closed loop poles of the continuous system was also 
written. The closed loop poles as calculated using the pro- 
gram track closely the eigenvalues derived from the exact 
time domain analysis. System instability occurs for some va- 
lues of filter inductance but the addition of feedforward 
eliminates the problem. 

8.1.1 Calculation of Eigenvalues . 

The method used in calculating the eigenvalues has been de- 
scribed earlier [15,16,17] . The state equations of the sys- 
tem are written for the buck regulator of Figure 34 . The 
state variables describing the system without feedforward 

x(l) = i LJ (input filter inductor current) 


are: 
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x(2) = v cl (input filter capacitor voltage) 
x(3) = i L (output filter inductor current) 
x(4) = v c (output filter capacitor voltage) 
x(5)=e E 
x(6) = e c 

As discussed in Chapter 7 e R is the voltage at the capacitor 
Cg, in the feedback loop while e c is the control voltage (in- 
put to pulse modulator) . Except for the order in which they 
appear, these are the same state variables as used in Chap- 
ter 7, where it is noted that the system is completely de- 
scribed by the above six state variables. 

The state equations without feedforward are written as 
in Chapter 7; 

During T n (SI on S2 off) 

UN (8-1) 

x == Fl»x + Gl*u 

During T rtc1 _ (SI off S2 on) 

0FP (8-2) 

x = F2*x + G2*u 
■■■* 

The matrices FI and F2 are identical to the ones in Chapter 
7 except for the following reordering; the third rows of 
equations (7-2) and (7-8) are moved to the sixth place with 
the other rows being moved up one place to get FI, F2 of 
equations (8-1) and (8-2), while the third columns of equa- 


t. 




2: 


tions (7-2) and (7-8) are moved into sixth place with the 
other columns being moved up one place in a manner similar 
to the rows. Matrices G1 and G2 of equations (8-1) and (8-2) 
are identical to the ones in Chapter 7 except that the third 
row of Chapter 7 is moved into sixth place with the other 
rows being moved up one place. The vector u of equations 
(8-1) and (8-2) is defined exactly as in Chapter 7. 

Figure 77 serves to establish notation regarding time 
instance t E etc. and the solutions to the state equations 
are written as in Chapter 7: 

During T Qpp 

x(t K + T) = $ p • x(t K ) + D p • u 
During T QN (8-3) 

S“K + T OFF + T) = • S (t K + T OFF> 


+ D. t • u 

The matrices , D F and D N are defined, as is time step 

T, exactly as in Chapter 7. 

From (8-3) the following is derived [15,16,17] : 

x(t K+ i } » * <f> p c — ^r) + Hi * d f ’ H (8 -4) 

+ D w • “ 

For the constant volt-sec. (VjT on ) duty-cycle control mode 
used in the stability investigation, the T 0 jj is fixed if Vt 
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is fixed and T OJ?F is determined by the point at which e c in- 
tersects the threshold voltage E T : 

6 4 

E_ = t $-,(6,1) * x(I) + l D_(6,J> • u{J) (8-5) 

1 1=1 * J=1 r 


Tqpf is thus a function of x(tg ) and (8-4) is a nonlinear 
equation [15,16,17]. Equations (8-1) - (8-5) thus represent 
exactly the nonlinear buck regulator system without feedfor- 
ward. 

For the purpose of stability analysis the steady state 
or equilibrium state is necessary [15,16] and this is der- 
ived as 

^ (t K+l J ~ £ (fc K } = -* 

(8-6) 

x* = $ K • $p ■ x* + ( 4 > K • D p + D n )u 

The nonlinear system is linearised for a small perturbation, 
around the steady state operating point [15,16,17]: 


. Sx(t R ) = x(t K ) - x* (8-7) 

and 

a S (t KH-l)> - K • 4 [*F • S (t K> + °F • S £ (t K> ' 8 - 8 > 

X* 
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Equation (8-8) is derived assuming that the small perturba- 
tion in the state does not change the value of Tq FIP [15,16]. 
The expression in the curly brackets of equation (8-8) is 
used to investigate stability: 

= ^ (8-9) 

Equation (8-9) thus describes a linear system which will be 
stable if and only if all the eigenvalues of if? are abso- 
lutely less than unity -- 

|X^(^)| <1 i = 1,2,3,... (8-10) 

or in other words if all the eigenvalues lie inside the 
unit circle [15,16]. 

The buck regulator system. with feedforward is described 
by a set of seven state variables as in Chapter 7. The nona- 
daptive feedforward circuit of Figure 67 was used in the 
analysis and the state variables then are: 

x(U = i L1 
x(2) - v G1 

x(3) = v f (voltage at capacitor C f ) 

x(4) = i h (8-11) 

x(5) = v c 

x(6) = ejj 




The state equations during T Qn and T 0ff are written as: 
During T 0N 

x = FI . x + Gl.ii (8-12) 

During T opp 

x = F2.x + G2.u (S-13) 

Matrices FI, F2, G1 and G2 of equations (8-12) and (8-13) 
are similar to the ones in Chapter 7 and are defined exactly 
as the ones in equations (8-1) and (8-2). The solution of 
the state equations, the definition of the equilibrium state 
and the matrix ijp are identical to the case without feedfor- 
ward. 

A computer program was written to calculate the matrix 
and its eigenvalues both without and with feedforward. It 
can be easily seen that by changing the single stage input 
filter parameters a set of eigenvalues can be calculated 
from which inferences regarding system stability can be 
drawn. 

8.1.2 Description of the Program Used in Calculating 
Eigenvalues 

A flow chart for the program is given in Figure 78 and a 
computer listing of both programs is included in the appen- 



There are essentially three steps involved in the 
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calculation of eigenvalues. 

First the approximate -steady state values are computed 
[15,16,17] , using an approximate value of the off time T 0FF : 

T OFF ■ T 0N ' (V I - E o )/E o (8-14) 


Equation (8-6) is used to calculate the approximate steady 
state values x* with ISML routine LEQT2F being used to solve 
the system of linear equations. It is to be noted that equa- 
tion (8-6) is not valid for the last state variable e c since 
both ^(6,6) and 0^6,6) equal unity, and equation (8-5) is 
used to calculate e„. 

The approximate steady state calculated is used to find 
the exact steady state values [15,16]. The best way of det- 
ermining the exact steady state is to determine the exact 
off time by iterative linearization { Newton's method) on 
the cycle-to-cycle matching condition for e c [15,16,17]. The 
vector ZT is defined : 


,K 


ZT = x(t K + T£ FF ) 


ZT - 4» f * £(tj,) + D 


U 


(8-15) 


with ZT( 6 ) = E t 

The cycle-to-cycle matching condition is expressed as: 




SMAT ~ x(6) 


4> n C6,D • ZT(I) + E t 


(8-16) 




+ l D W (6,J) 
J=1 W 



The value of SMAT is calculated using the available value of 
if SMAT is smaller than a certain error tolerance then 
the value of used is accurate. If SMAT exceeds the to- 

lerance then Newton’s method is applied to calculate a bet- 
ter estimate of T 0 p P . The steps involved in the calculation 
are i 

1. Calculate ZT and SMAT. 

2. If SMAT ,is not less than specified tolerance, then 
perturb T OFF by the amount DT OFjn . 

3. Calculate values of state variables corresponding to 
the perturbed value of off time, using equations 
(8-5) and (8-6). 

4. Calculate SMATN. 

5. Calculate a better estimate of T QFF using Newton's 
iteration (15,16,17). 

6. Calculate a better estimate of the state variables. 

7. Calculate SMAT and go to step 2. 

The last step is to calculate the matix and its ei- 
genvalues using the exact steady state values computed 


0 

D 
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Figure 78: Flowchart for stability analysis 
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above; Figure 79 is a flowchart for the calculation. The 
following are defined: 


FUN1 (x* + Dx(I)) = 


• X + 

Dp • u 

X* + Dx(I) 

(8-17) 

FtM2(x* - Dx (I) ) = 


• X + 

d f * £ 

x* - Dx (X) 

(8-18) 

RfUNI - F0N2) 

TEMP 2 = — ~ 

(_ 2 • Dx (1) 

* * * 

(FUN1 - FUH21 
2 • Dx (6} J 6X6 

(8-19) 


ESI = $ • TEMP 2 


( 8 - 20 ) 


TEMP2 is thus the matrix of partial differentials used in 
equation (8-8). 

The matrix TEMP2 is evaluated one column at a time. For 
the I'th column a perturbation DX(I) is made in X(I). Two 
new sets of steady state values XA and XB are defined, as is 
shown in the flowchart. For the new set of values XA it is 
necessary to calculate an exact value of off time TFFC1 ir 
order to obtain ^ , Dp. The error function SETA is defined: 

6 4 

ZETA = t <M6,I) • XA { t J + l D^fCvJ) u(J} ~ E-, (8-21) 

1=1 * J=1 F T 
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Figure 79: continued 
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Calculate <j> pT , D pT 
and 2ETAN 
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SLOPE = {ZETAN-2ETA)/ 

dt off 

TFFC2 = !TFFC2 - 

(2 ETA/SLOPE) 




Calculate <j> p / D p 

IT=IT+1 




Figure 79: continued 












The value of 2ETA is calculated using the available value of 
TFFC1, if it is within the specified tolerance then the va- 
lue of TFFC1 is the value needed. If SETA is not within to- 
lerance then Newton* s iteration [15/16] is used to obtain a 
better estimate of TFFC1/ as is evident from the flowchart. 


The values of 
is evaluated. 


, D F are obtained and then the vector FUN1 


A similar procedure is used to obtain the exact off 
time TFFC2 for the other set of values XB. After calculation 
of FUN2 an estimate for the I'th column of TEMP2 is ob- 
tained. 

In the calculation of the partial derivatives it was 
observed that a better estimate of the derivatives can be 
obtained using the Richardson's extrapolation method [18]. 
The derivative f 1 (a) is sometimes approximated as: 


f (a + h) - f (a - h) 


8-22 


where h is the step siae, 


A better estimate may be obtained by evaluating A(h) and 


A(rh) 


f (a + rh) - f (a - rh) 


r*' 


| 





where r is usually 0.5 [18]. 
Combining A(h) and A(rh) thus: 


B (h) — ? 2 &(h\ (8*24) 

1 - r 2 

gives a better estimate, B(h), of the derivative f'(a) [18]. 
Thus by combining the value of A(h) obtained from equation 
(8-22) with two values for the step size h, it is possible 
to get a more accurate estimate of the derivative. 

In the program four values of step size were used. The 
perturbation DX(I) is assigned four values and for each va- 
lue the vectors FUMI and FUN2 are calculated. This results 
in four estimates for the derivatives in the I ’ th column of 
TEMB2 . Combining these four estimates as in the program 
[18] gives an accurate estimate of the I'th column. 

The whole procedure is then repeated for the next co- 
lumn, thus an accurate set of values for the derivatives is 
obtained. Multiplying TEMP2 by results in the matrix 
[15,16,17], Finally the eigenvalues X of the matrix are 
obtained using the IMSL routine EIGRF. 

The eigenvalues with feedforward are calculated in an 
identical fashion, the only real difference being the change 
in the order of the system from six state variables to sev- 


D 

0 


0 


0 

0 

D 


0 

D 

D 


231 


8.1.3 Calculation of the Closed Loop Poles of the System 
The closed loop poles of the system of Figure 34 can be cal- 
culated. The expression for the closed loop gain v 0 /Vj. can 
be easily written using Figure 33 for the case without feed- 
forward: 




1 + F DC- P M- P 12 + F 3- P AC- P M- F 22 


(8-25) 


with Fjj , etc. as defined in Chapter 5. 

Equation (8-25) can also be used to calculate the closed 
loop poles of the system without input filter; all that is 
necessary to do is to set H = 1 and 2 = 0, as discussed ear- 
lier. 


Expressions for Fj^, A etc. are substituted into equa- 
tion (8-28) for the case without input filter. The resulting 
expression in the denominator is then a polynomial of the 
fourth degree: 


p 1 (s) = s 4 (k 7 k 6 lc) + 

S 3 (K 7 K 6 Kg + K ? CjLC + KgK^CR^) + 

S 2 (K 7 CiK 8 + K 7 K 6 + K 3 K 2 CE c + K 5 K 4 
+ K 5 nL C 1^) + 

S{K 7 C[ + K $ K 2 + K 3 K x C R c + nL K 5 ) 


(8-26) 
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The constants k ? , k 3 etc. are defined as functions of the 
regulator parameters like inductance L etc. The appendix 
contains a listing of the program used and the constants are 
defined there. An IMSL routine ZPOLR is used to calculate 
the roots of Pj which are the closed loop poles of the sys- 
tem without input filter. 

For the case with input filter the same equation, 
(8-25), is used and substitution leads to a sixth degree po- 
lynomial in the denominator: 



P 2 (s) * s 6 (k ? k 6 y 8 ) + 

s s <K 7 C'Y B ♦ k 7 k 6 y 9 + k 5 k 4 y 5 ) +• 

S 4 (K 7 C-Y 9 + K 7 K 6 Y 10 + K 3 K 2 Y 1 + KgllLYj 

+ k 4 K 5 Y 6 ) + 

S 3 (K 7 c i Y 10+ K 7 K 6 Y u + K 3 K 2 Y 2 + K 3 K 1 Y 1 + 
K 5 nLY g + K 5 K 4 Y ? ) + 

s 2 (« 7 0iY 11 + k 7 k 6 y 12 + k 3 k x y 2 + k 3 k 2 y 3 
+ K 5 nLY 7 + k 5 k 4 y 4 ) + 

S(C*K 7 Y 1 2 + K 3 K 1 Y 3 + K 3 K 2 Y 4 + nLK 5 Y 4 ) 

+ K 3 K 1 Y 4 


(8-27) 
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The constants Yj etc. used in equation (8*27) are again de- 
fined in the program. The same routine ZPOLR is used to ob- 
tain the six closed loop poles. 

8.1.4 Location of Eigenvalues and Poles as a Function of 
Input Filter Parameters . 

The computer programs written were used to locate closed 
loop poles and eigenvalues as function of input filter in- 
ductance. The buck regulator of Figure 34 was used and the 
parameters of the regulator were as defined in section 6.1 
except for the following changes: 

Vj = 25 v R 13 = 200 Kohras 

C a = 100 pf R l = 10 ohm 

The values of R ^ and C ^ were changed so that the compensa- 
tion loop is largely ineffective, [7,8]. This was done in 
order to facilitate making the measurements discussed in 
section 8.2. It was found difficult to achieve instability 
in the high loop gain system used in making the measure- 
ments. Making the compensation loop largely ineffective re- 
duced the loop gain so that with a large input filter the 
system was made unstable. 

A single stage input filter with the following parame- 
ters was used: 

R = 0.2 ohm LI variable Cl = 220 microF 

l 


First the eigenvalues and closed loop poles without input 
filter and feedforward were calculated using the eigenvalue 
and the closed loop programs. Figures 80 and 81 are the 
plots of the eigenvalues and closed loop poles respectively. 
It is noticed that there are four poles and four eigenvalues 
since without input filter the system is of the fourth ord- 
er. 

The eigenvalues and poles are: 


a 0.844 

* 10' 7 

s 2 = -55,000 + jO 

= 0.12 


s z = -50,000 + jO 

= 0.98 + 

jO. 103 

s 3 = -346 + j2240 

= 0.98 - 

jO. 103 

s* = -346 - J2240 


The eigenvalues and the closed loop poles are related 
[15] : 

z = e sT (8-28) 

s = o + jw 

Thus 

| z | - e aT (8-29) 

^z = u)T 

Knowing the location of the closed loop poles it is thus 
possible to locate the eigenvalues if a value of T is known. 
In this case T is the switching period and its value was 



found to be 41.85 microsec.; thus equations (8-28) and 
(8-29) could be used to crosscheck the results from the ei- 
genvalue program once the closed loop poles are known. Equa- 
tion (8-29) was used to calculate the eigenvalues and the 
calculated values were found to be very close to the values 
from the eigenvalue program. The two poles on the real axis 
correspond to the two eigenvalues close to the origin while 
the two complex poles result in the two complex eigenvalues. 

Next the single stage input filter was put in and the 
eigenvalues and closed loop poles calculated using the res- 
pective programs. Figures 82 and 83 are the plots of the ei- 
genvalues and the closed loop poles respectively. The values 
of the filter inductance LI used were — 

50, 150, 325, 450, 650, 800, 1000, 1425 and 1800 microH 
Table 1 lists the eigenvalues and closed loop poles, for 
different values of LI. Starting at 50 microH it is seen 
that the complex eigenvalues move to the right, closer to 
the unit circle as inductance is increased. The two eigen- 
values on the real axis z x and z^, which correspond to z x 
and z a on Eigure 80 , do not move with change in inductance . 
The complex eigenvalues with a real part of 0.98 z s and z e , 
which correspond to s 3 and on Figure 80 , are pushed -out 
of the unit circle for a value of LI = 800 microH. For va- 
lues above 800 microH the eigenvalues do not move much; the 





Figure 80: Eigenvalues without input filter 
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real part stays at around 1.003 for the ones outside the 
unit circle and around 0.97 for the ones inside. This may 
be used to explain the bunching up of the eigenvalues for 
large values of inductance. Thus the eigenvalue plot of Fig- 
ure 82 predicts that the system will be unstable for values 
of filter inductance above 800 microH. 

The plot of t' s closed loop poles. Figure 83 , shows 
similar behavior. The two poles on the real axis s l and s a 
are not changed as inductance increases, however two complex 
poles s s and s e move into the right half plane, indicating 
instability. The poles move across the origin at 800 microH 
which checks well with the eigenvalue plot. Equation (8-29) 
was used to calculate the eigenvalues from the closed loop 
poles and the calculated eigenvalues agreed closely with the 
results from the eigenvalue program for all the values of LI 
that were used, thus providing a crosscheck. 

The other pair of complex poles s 3 and s^ in Figure 83 ini- 
tially move closer to the right half plane, then move away 
from it for values of LI greater than 450 microH. This be- 
havior should be reflected in the eigenvalue plot for eigen- 
values z 3 and z however it is noticed than in equation 
(8-29) the value of T is so small that for a change in the 
real part from -300 to -500 the change in the eigenvalues 
will be about 1%, which is too small a change to show up in 
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TABLE 1 


Eigenvalues and Closed Loop Poles for various values of LI, 
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TABLE 1 
continued 



Eigenvalues 


-0.003 + j0 
0 . 1233 + j 0 
0.971 * 
j 0.132 ■ 
0.9897 ± 
j 0.1026 


-0.003 + jO 
0.1233 + jO 
0.97 ± 

3 0.1162 
0.9975 ± 
j 0.0954 


-0.003 *i* jo 
0.1232 + jO 
0.97 ± 

3 0.111 
0.9996 ± 
j 0.089 


Closed Loop Poles 

-55,045 + 

30 

-SO, 005 + 

30 

-323 A 


3 3165 


-214 a 


3 2237 


-55,045 + 

30 


-50,003 + jO 
—419 & 
j 2684 
-49.8 ± 

3 2189 


-55,045 + jO 
-50,005 + jO 
-483 ± 
j 2546 
43 . 6 £ 
j 2069 
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TABLE 1 
conti vrned 



Eigenvalues 

Closed Loop Poles 


-0.003 + jO 

-55,046 + jO 


0.1233 + jO 

-50,005 + jO 

LI = 1000/* H 

0.971 £ 

-513 ± 


j 0.1078 

j 2455 


1.0002 ± 

98.3 ± 


j 0.082 

j 1909 


-0.003 + jO 

-55,046 + jO 


0.1233 + jO 

-50,005 + jO 

LI = 1425/t H 

0.972 ± 

-517 ± 


j 0.1052 

j 2371 


1-003 + 

133 ± 


j 0.071 

j 1650 


-Q.0Q3 + jO 

-55,046 + jO 


0.1233 + jO 

-50,005 + jO 

LI - 1800 /AH 

0-972 ± 

-511 ± 


j 0.1055 

j 2339 


1.004 ±' 

141 ± 


j 0 . 063 

j 1487 











Figure 82 


Eigenvalues with input filter, without 
feedforward 
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Figure 83: 


Closed loop poles with input filter, without 
feedforward 



the eigenvalue plot. The bunching up of the closed loop 
poles for higher values of LI is also noticed; this provides 
another crosscheck with the eigenvalue plot. 



Thus it is seen that the eigenvalues and the closed 
loop poles match both qualitatively, as examination of Fig- 
ures 82 and 83 show, and quantitatively through the use of 
equation (8-29) to crosscheck the eigenvalues. It is also 
noticed that system instability is predicted for values of 
LI above 800 microH. 

Lastly the eigenvalues were calcul^t^d for the case 
with feedforward using the program. The nonadaptive feedfor- 
ward circuit of Figure 67 was used with the following param- 
eters: 

= 27 microF — 5100 ohms 

R fz variable from 210 ohm to 300 ohm. 

Figure 84 shows the plot of the eigenvalues with feedforward 
for the following values of R f2 , with the value of LI kept 
fixed at 1425 microH: 

3£f2 = 210, 220, 230, 235, 240, 260, 270 and 300 ohms 
The value of R fa was changed to see the effect of feedfor- 
ward on the eigenvalues. Table 2 lists the eigenvalues for 
different values of R f a . It is noticed from Figure 84 that 
there are seven eigenvalues and that the three on the real 
axis z-j, , and z 3 do not change as R f2 , changes in value. 
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As R f2 is increased one pair of cc-nplex eigenvalues and 
z s move outside the unit circle while the other pair zq and 
z 7 move inside. For the value of LI used the closed loop 
poles and eigenvalues without feedforward clearly predict 
instability, whereas with feedforward it is seen that the 
system will be stable for some values of Rf. a . For the regu- 
lator under consideration the value of is calculated to 
be 220 ohm from the feedforward design as obtained in Chap- 
ter 5. A value less than the design value causes the eigen- 
values 2 6 and z 7 to be outside the unit circle while it is 
noticed that a value higher than the design value also caus- 
es instability in the system. For the design value of 220 
ohms all the eigenvalues are inside the unit circle thus 
clearly indicating that the addition of feedforward stabi- 
lises a system that was unstable due to the interaction bet- 
ween the input filter and the regulator control loop. 

As R fa is increased further the real parts of the com- 
plex eigenvalues do not change much but eigenvalues s ^ and 
Zg do move outside the unit circle. For R fg between 220 and 
260 ohm all the eigenvalues are inside the unit circle, with 
the optimum value being around 230 ohm. 

In order to see clearly the effect of changing the 
feedforward loop gain the open loop gain and phase margin 
were calculated for various values of R f2 . 


The parameters of 


o 
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the buck converter used in the calculation were the same as 
used in the eigenvalue calculation, with LI = 1425 raicroH. 
Figure 85 (a) - (f) shows the plots of the open loop gain 
and phase margin for three values of R fa . The solid lines 
on the figures are the plots without feedforward and with 
feedforward (R fa = 220 ohms), whereas the triangles repre- 
sent values with R fa either lower than or higher than the 
design value of 220 ohm. Exact compensation of the input 
filter interaction is achieved with R f2 = 220 ohm, with any 
value other than 220 ohm the compensation is not complete. 
Decreasing the feedforward loop gain (210 ohm) results in a 
loss of open loop gain at the filter resonant frequency. In- 
creasing the feedforward loop gain results in a loss of 
phase margin at the resonant frequency. For a value of 240 
ohm the open loop gain obtained is close to the gain with 
220 ohm but the phase margin is decreased. Increasing the 
feedforward gain further (300 ohm) makes the situation worse 
as Figure 85 (f) shows. The open loop gain in Figure 85 (e) 
is decreased to around 8 db at a frequency of about 300 Ha 
while at that frequency the phase margin is almost aero, in- 
dicating that the syst*' n is close to instability. This is 
confirmed on examination of Figure 84 which shows that the 
system is unstable for R fa = 300 ohm. it should be noted 
that the open loop gain with an input filter and R.^220 ohm 
is identical to that without input filter. 


It is thus concluded that exact compensation is ob- 
tained at the design value of the feedforward loop gain. In- 
creasing or decreasing the feedforward gain from the design 
value results in incomplete compensation. There may be no 
compensation achieved at all if the feedforward gain is far 
different from the design value. 

It is noticed in Figure 84 as in the other eigenvalue 
plot that the eigenvalues do not seem to move very much, for 
example • in Figure 84 the real part of the pair z ^ and z s 
changes from about 0.9885 to 0.9975 while from Figure 82 it 
changes from about 0.98 to 1.004. It is noted however that 
from Figure 83 the above change in real part corresponds to 
a change in the real part of the poles from about -323 to 
+140. This may be explained by noting as before that T in 
equation (8-29) is fairly small so this change in real part 
will correspond to a very small change in the real and imag- 
inary parts of the eigenvalues. 

From the results presented in this section it is seen 
that the eigenvalue program is a powerful tool for the ana- 
lysis of stability. It is also seen that for the regulator 
used, a high value of input filter inductance will cause 
instability, but the addition of feedforward restores sta- 
bility. Experimental data that strongly support this analyt- 
ical prediction are presented next. 
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Figure 85: 


(a) Open loop gain: calculated values without 
feedforward, with Rj? a =220 ohm and Rp a =210 ohm 
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Figure 85: (c) Open loop gain: calculated values without 

feedforward, with R f2 =220 ohm and R^ a =240 ohm 
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Figure 85: (d) Open loop phase margin: calculated values 

without feedforward, with =220 ohm and 
K„„=240 ohm 
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Figure 85: (e) Open loop gain: calculated values without 

feedforward, with R f ^=220 ohm and =300 ohm 
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TABLE 2 

Eigenvalues for different values of , 



Eigenvalues 

R^ - 210 ohm 

0*00116 + jO 
0.1232 + jO 
0.9998 + jO 
0.9974 * j 0.0738 
0 . 9883 ± j 0.101 

R f2 = 220 ohm 

0.00017 + jO 
0.1282 + jO 
0.9998 + jO 
0.9963 ± j 0.0742 
0.9907 A j 0.1006 

R fa = 230 ohm 

0.00122 + jO 
0.1282 jO 

0. 9998 + jO 
0.9955 A j 0.0743 
0.9922 ± j 0.10121 






Eigenvalues 
-0.00025 + jO 
0.1282 + jO 
0.9998 + jO 
0.9952 ± j 0.0742 
0.9931 + j 0.1006 
-0.00011 + jO 
0.1282 + jO 
0.9998 + jO 
0.9948 ± j 0.0743 
0.9936 ± j 0.10049 
0.00039 + jO 
0.1283 + jO 
0.9998 + jO 
0.9946 + j 0.0744 
0.9941 + j 0.1011 


D 

D 

D 






D 
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TABLE 2 
continued 


R f2. ~ 270 ohni 


R._ = 300 ohm 
f 2 , 


Eigenvalues 
0.00023 + jO 
0.1283 + jO 
0.9998 + jO 
0.9933 + j 0.0742 
0.9964 + j 0.1008 
0,00029 + jO 
0.1285 + jO 
0.9998 + jO 


0.98952 + j 0.0724 
1.0054 + j 0.1036 
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8.2 EXPERIMENTAL VERIFICATION OF STABILITY ANALYSIS . 

In this section experimental data pertaining to the stabili- 
ty of the system are presented. The buck regulator of Figure 
34 was used with the parameters as defined in section 8.1.4 
with the exception that the input filter inductance LI had a 
constant value of 1.425 mH. The adaptive feedforward circuit 
of Chapter 5 was used, however it is noted that if the input 
voltage is kept constant, as it was in the following mea- 
surements, the nonadaptive and the adaptive feedforward cir- 
cuits provide the same gain. 

The regulator was set up with the HP network analyser 
to measure open loop gain exactly like in Chapter 6 
The open loop gain and phase margin were measured without 
input filter, with and without feedforward. Figure 86 shows 
the measured values of the open loop gain and phase margin. 
It can be seen that without the input filter the gain dips 
down around 320 Hz — this may be explained by noting that 
the values of R l3 and C 2 are such that the compensation loop 
[7,8] is largely ineffective. Thus the open loop gain and 
phase margin correspond to a value of of - 0.355 as shown in 
the literature [7,8]. 

The. addition of the input filter causes the gain to 
fall below 0 db and when the gain is close to 0 db the phase 
margin is only about 10 degrees at a frequency of about 250 
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Eigure 86: (a) Open loop gain.: measured values without: 

input filter, with (») and without feedforward 
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Figure 86: 


(b) Open loop phase margin: measured values 
without input filter, with («) and without 
feedforward 





261 


Hz. The system is thus marginally stable and a small signal 
250 Hz sinusoidal disturbance is injected into the loop to 
induce unstable operation and it was observed that the out- 
put voltage has oscillations. Figure 87 shows the output 
voltage ripple -- it is clearly seen that the oscillations 
are around 1.2 v in magnitude. It is concluded that the sys- 
tem is clearly unstable. Analysis presented earlier showed 
that the syatem would be unstable for LI - 1425 microH while 
experimental results show a marginally stable system,, thus 
lending support to the analysis., 

Feedforward was added to the circuit and measurements 
of the open loop gain and phase margin were made. It is seen 
clearly from Figure 85 that the addition of feedforward re- 
moves the instability — this is confirmed by recording the 
output voltage ripple. Figure 88 shows the output voltage 
ripple and it can clearly be seen that the oscillations that 
were present earlier without feedforward have been eliminat- 
ed, making the system stable. The 250 Hz oscillation in Fig- 
ure 88 is the injected disturbance. The analytical predic- 
tion that feedforward elxminates instability caused by input 
filter interaction is thus strongly supported by experimen- 
tal data. 
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(b) 

Figure 87: Output voltage ripple without feedforward: (a) 

Y: 0 . 2v/div X: 5 msec/div. (b) Y:0.2v/div X:10 
msec/div. 
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Output voltage ripple with feedforward: (a) 

Y: 0 . 2v/div X: 5 msec/div. (b) Y:0.2v/div X:10 
msec/di v. 
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S. 3 CONCLUSIONS 

It is seen that the eigenvalue analysis program is a power- 
ful tool for the analysis of stability. Using a buck regula- 
tor system a stability analysis was performed for various 
values of input filter inductance and it was observed that 
system instability was predicted for large values of induc- 
tance. A valuable crosscheck was provided by a program that 
calculated the closed loop poles and it was observed that 
the closed loop poles and the eigenvalues track very closely 
for all the values of inductance used, including those va- 
lues for which system instability is predicted. 

The addition of feedforward removes the instability at 
large values of inductance as the eigenvalue analysis pre- 
sented shows. It is noticed that complete compensation of 
the input filter interaction is obtained for the design va- 
lue of feedforward gain; any other value of feedforward gain 
results in incomplete compensation. It was demonstrated both 
from the eigenvalue analysis and the open loop gain and 
phase margin plots that a too large or too small value of 
feedforward gain would not be able to remove the instability 
caused by the input filter interaction. 

Experimental measurements were made to confirm the ana- 
lytical results presented. It was noticed that without feed- 
forward the output voltage has large magnitude oscillations 
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indicating system instability. These oscillations occur at a 
value of input filter inductance for which the eigenvalue 
analysis also predicts instability thus tying in the analy- 
sis with experimental measurements. The addition of feedfor- 
ward (with the gain set to its design value) completely re- 
moves these oscillations as measurement data show. Since the 
eigenvalue analysis also predicted stability with the addi- 
tion of feedforward, it is thus concluded that the addition 
of feedforward does result in removing the system instabili- 
ty caused by input filter interaction with the regulator 
control loop. 
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Chapter IX 


EXTENSIONS OF THE CONCEPT OF INPUT FILTER 
COMPENSATION TO OTHER CONVERTER TYPES 

The concept of feedforward loop compensation to mitigate 
converter — input filter interaction is extended to several 
other cases, three of which are discussed in this chapter. 

9.1 EXTENSION TO OTHER TYPES OF REGULATOR CONTROL 

The concept of feedforward compensation for the input filter 
has been developed with reference to a buck regulator. Ex- 
tensive experimental data, together with analysis, confirm 
the validity of the concept. This was done exclusively on a 
buck converter with multiple- loop control and constant 

volt-sec duty cycle modulation. Use of other types of con- 
trol are discussed in this section. 

9.1.1 Feedforward c ompen s at i on for multiple- loop control 
using other types of duty- cycle control 

The feedforward design was presented, in Chapter 5 where it 

was noted that the gain of the feedforward loop depends on 

the transfer function, of the pulse modulator , F^ . The 

transfer function F M is dependent on the type of duty-cycle 

control used [2,8]. In this dissertation the constant volt- 

sec. type of duty-cycle control was used. However, it is 
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easy to see that for other types of duty-cycle control/ for 
example for constant frequency control or constant c'*£ time 
control, the feedforward design procedure outlined in Chap- 
ter 5 is still valid. The feedforward design process pre- 
sented earlier is general since it treated the pulse modula- 
tor in terms of its transfer gain only. Thus other types of 
duty-cycle control can be easily incorporated in the feed- 
forward compensation after the transfer gain of the duty-cy- 
cle control is identified. 

As an example for multi-loop constant frequency control 
the following pulse modulator gain is identified [8]: 


where 


P = 
M 


2H 4 C' 

nM 


M = V t C1-2D)T 


(9-1) 


In equation (9-1) T p is the constant switching period. Sub- 
stitution into the feedforward loop gain expression from 
Chapter 5 results in: 

C = ■ VnM 

C 2 V *2R„C' Z} 

O 4 1 

For a constant supply voltage D and M would be constants and 
thus the feedforward gain would be a constant. An adaptive 




feedforward circuit that tracks the supply voltage can be 
realized using equation (9-2). 

For multi-loop constant off time control F^ is identi- 
fied as [8] : 


2R 4 q 


V T 
I OFF 


(9-3) 


where T OF pis the constant off time. Sustitution into the 
feedforward loop gain expression results in* 


c 2 = 


2 

~ D nV I T OFF 

Wi 


(9-4) 


For a constant supply voltage -D would be a constant and thus 
the gain would be a constant. An adaptive feedforward cir- 
cuit that tracks the supply voltage variations can be real- 
ized from equation (9-4) . 


9.1.2 Feedforward compensation for single-loop control 


A multi- loop feedback controlled buck regulator was used in 
the feedforward design discussed in Chapter 5. However, it 
is noted that in the design • procedure the feedback loop 
transfer functions played no part in the feedforward design. 
This is evident since the feedforward gain is independent of 
the feedback loop transfer functions. 
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For single loop feedback control the regulator can be 
represented by a generalized block diagram as shown in Fig- 
ure 89 . The feedback loop has as its input the output vol- 
tage ripple and it can be represented as F S/ as shown in 
Figure 89. Expressions for and F m are given in the litera- 
ture [8], and it is noted that the transfer function F B es- 
sentially processes the small signal variation in the output 
voltage and feeds it to the pulse modulator exactly like F Dc 
of Figure 28 . A feedforward loop that will sense the input 
filter state variables/ process this information and feed it 
to the pulse modulator can thus be designed in a manner ex- 
actly as presented in Chapter 5. The transfer function F M 
would be a constant [6] and would be a part of the feedfor- 
ward gain, exactly as it is in the design presented. 

As an example the following is identified for single 
loop constant frequency control [6]: 

P = (1/A T ) (9-5) 

where A,. T„ are constants, 
o p 

The error processor consists of an amplifier and a lead-lag 
compensation network and its transfer function F E is availa- 
ble in the literature [61 . Since F M is a constant the de- 
sign of feedforward loop is relatively simpler as substitu- 
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Figure 89: Generalized small signal model for single loop 

control 
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tion into the feedforward gain expression of Chapter 5 
shOWS: 


C 2 


-VAT 
o op 



(9-6) 


For constant supply voltage the gain is a constant. An 
adaptive feedforward circuit that tracks the supply voltage 
variations could be developed exactly as shown in Chapter 5. 

It is thus seen that an extension of the feedforward 
concept to single loop control is easily affected. 


9.2 FEEDFORWARD COMPENSATION FOR THE BUCK-BOOST AND BOOST 
REGULATORS 

Extension of the f eedforward concept to the two other types 
of regulators — boost and buck-boost can be obtained follow- 
ing the outline presented in Chapter 5. A feedforward design 
for the buck-boost regulator is presented. 

The state space model of the buck-boost regulator is 
shown in Figure 90. In the figure the feedforward loop has as 
its input the input filter state variables (inductor current 
and capacitor voltage). These two are multiplied by the 
transfer functions c^ and c s whose properties have, yet to be 
determined. This state space model is based on the model 
for the buck-boost power stage developed in Chapter 4. The 
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feedback loops form the standardised control module (SCM) 
developed earlier [2,7,8] and discussed in Chapter 5. The 
pulse modulator has as its transfer function F^ [8] and has 
as its input information regarding the state variables of 
both the input filter and the output filter. The total 
state feedforward/feedback error signal is thus used to mo- 
dulate the duty cycle of the switch for loop gain correc- 
tion. 

The generalized small signal model for the buck-boost 
regulator is developed next and used to write the transfer 

A A 

function v /v . 

o x 


9.2.1 Generalized small signal model for the buck-boost 
regulator 

The generalized small signal model is shown in Figure 91. The 
regulator is modeled according to the three basic functional 
blocks ; power stage, error processor and duty cycle pulse 
modulator, as in Chapter 5. It is noticed that in this case 
the flux is a state variable, because it is continuous un- 
like the inductor current. The transfer functions % '®ac 
F dc play no part in the feedforward design process, as ear- 
lier in Chapter 5. 

The power stage transfer functions F^ etc, are written 
as in Chapter 5, using the following equations: 


T 11 V I + T 12 d = v o 


T 21 v I + T 22 d = * 
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‘ r 31 v I + T 32 d " V CI 


T 41 v I + T 42 d = *11 


19 - 7 ) 


Thus it is seen that 


and 



T 

T 


HF 


11 


11 

12 


12 

A 


(9-8) 


(9-9) 


The transfer functions '^ 21 ' T 31 and 1 44 are developed 
with d «0 and the others with Vj - 0 . The starting point for 
the derivation is again the small signal equivalent circuit 
model for the buck-boost regulator developed in Chapter 4. 

A 

Evaluation of transfer functions with d = 0. 

ft 

The equivalent circuit model with d = 0 is shown in Figure 
92 In Figure 92 the input filter has been replaced by its 
forward transfer function H(s) and its output impedance 
Z(s), the expressions for H and Z being the same as in Chap- 
ter 5, From the equivalent circuit of Figure 92 the follow- 
ing are derived, assuming that R C1 - 0 : 


H*D 1 *D' *1^(1 + s CR c ) 
T 11 = A 

H-L^D^l + sCRjJ 

t 21 = ’ O 

s 


(9-10) 






















41 ~ (R^ + sLl) A 


where: 

a l = s2l s ck l + sC V R s + D ' R o + 

+ R s DD'R c + D ,2 R l 

A = + D 1 2 Z{1 + sCRjJ (9-11) 

DN s 

D 1 = — ' w ^ ere D is the duty cycle. 

P 

*L + R C * *L 

Using the above equations the following can be derived: 


F H = D 1 D’^d + SCRJ 

L s D 1 <1 + SCR^) 

P 21 I 

s 

F 31 = a l 

A/H - a 1 
F 41 (R^ + S &1) 


0 - 12 ) 


The ESR of the capacitor R Cl has been neglected in the deri" 
vation for the same reasons as outlined in Chapter 5, while 
equation (5-9) is used in the derivation of T A1 . 
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Evaluation of the transfer functions with 


The small signal equivalent circuit with v £ =0 is sh< 


n 


in Figure 93 The following are derived from the equival 


)<?> 
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p = V q (X + sCR l )D 1 2 2(DR + DSL S - D' 2 ^) + (DR - D-^Z) 

dd ' 2 d i r l 

F 42 = “ F 32^^ R L1 + SL ^ 

(9-15) 

where R = R^+DD 1 R c (9-16) 

The feedback transfer functions Eg, E^ c and F DC consti- 
tute the two loop SCM feedback control scheme reported ear- 
lier (2,8). The transfer functions are defined in the liter- 
ature [2,7,8], and it is noted that they do not play any part 
in the design process, as in Chapter 5* The pulse modulator 
transfer function F^ depends on the type of duty cycle con- 
trol used [7,8] and it will be seen that the feedforward 
loop gain depends on as earlier in Chapter 5. 

9.2.2 Design of the Feedforward Loop 

A 

With Vj - 0 the following equations are derived from 
Figure 91 : 

♦ Vci + Vli 1 F M= <* (9 - 17 > 

and after substitution 




F 

(112, F 
v a } m 


V 


X 


F 32 

1 “ C 2<~F> F m - C 3HT )P M 


(9-18) 


In the absence of feedforward i.e. with c a = 0 = c 3 it can 
be seen that 



’o' 1 


+ SCR ) 
c 


j -3 


2 Z + d' 2 r^ - DsL £ 


Pm 


DD' A 


(9-19) 


The effect of peaking of Z is to cause a reduction in the 
gain of the transfer function and thus also in the open loop 
gain. It is also noticed that equation (9-18) shows the po- 
sitive aero term in the transfer function. 

With the addition of feedforward the effect of peaking 
of Z could be avoided by a proper choice of feedforward loop 
gains c a and Cg. Two approximations are made to facilitate 
the design: 

1. The effect of the positive zero is assumed to be neg- 
ligible at the frequency at which Z peaks. Thus it is 
assumed that: 

D f 2 R l - DsLg « (9-20) 


This is not an unrealistic assumption since Z peaks 


at a much lower frequency. 
2. It is also assumed that • 
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at the frequency at which 2 peaks. Again this is not 
an unrealistic assumption given the low value of R. 
Using the above assumptions, the following choice is made: 


c 3 


= 0 


-DDi 


c 2 " V F,, 

O M 


(9-22) 


Substitution into equation (9-18) leads to: 



V Q ( l+SCBj.) {- Dl 2 Z+D ' 2 E l }p m 
DD'i 



a l DD l 



r lV 

i 2z+d 



(9-23) 


which leads to 


A 


v o V o (1 + sc V d ’Vm 



(9-24) 


It is seen from equation ( 9 - 24 ) that the transfer function 
is now independent of the peaking effect of 2. It is also 
noted that equation (9-24) is valid only at low frequencies 
where the two assumptions made earlier hold. The feedforward 
loop gains are very similar to the gains developed for the 
buck regulator, and are identical if N g = N p . It is also 
noted from equations (9-15) and (9-18) that the gain of F^ 
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would be fairly small at frequencies other than the resonant 
frequencies at which S peaks. Thus the addition of feedfor- 
ward would not affect/ in any noticable manner/ the open 
loop gain and phase margin at any frequency other than those 
at which 2 peaks. Thus the feedforward compensation is ef- 
fective for most input filter designs because the filter re- 
sonant frequencies are relatively low in comparison to the 
positive zero. 

9.2.3 Analytical and Experimental Verification 
The buck-boost regulator shown in Figure 94 was used to ob- 
tain verification of the feedforward design. A multi-loop 
standardized control module (SCM) type feedback control was 


used [2/7/8]. The parameters 

of 

the converter are 

as fol 

lows ; 




Power Stage Parameters: 




V r - 20 v V 0 = 28 v 

D 

= 0.5833 P Q = 28 

watts 

L s = 220 microH 

C 

« 300 microF 


R c = 0.05 ohm 


= 28 ohm (load) 

■ 

N s = N p = 33 

R S 

= 0.087 ohm 



Pulse Modulator Parameters ; 

-3 

M “ v i T oi r 0 .5*10 v-sec. 
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Control Circuit Parameters: 


E e = 7V 

Rg = 47.5 Kohm 

n = 0. 667 


R u =43.2 Kohm 
R* = 40 Kohm 
c| = 5600 pf 


R 2 = 15 Kohm 
R s =1.1 Kohm 
c' 2 = 32000 pf 


Input Filter Parameters 

Rj =0.2 ohm Lj = 325 microH Cj = 200 microF 

R 2 , = 0.02 ohm L 2 = 116 microH C 2 =20 microF 

R 3 = 0.075 ohm ( ESR of C A ) 


The regulator was operated in a predetermined duty cycle 
control mode (constant Vj T QN control ) [2,7,8]. 

The open loop gain can be written from Figure 91 as 


G t (s) 


{ 


12 

F ( + p F ( 
*D<r A ' A 


s§4 


M 


2_ c ( 32^ p 
1 c 2 { A ; M 


(9-25) 


It is noted that in deriving the above equation c* 3 is set 
equal to zero since that is the design value. Equation 
(9-25) was used to plot the open loop gain and phase margin 
without input filter, with and without feedforward. The 
transfer functions F fc- , F, , F A/ , and F M constitute the feed- 
back and are defined in the literature [2,7,8], while , 
F 22 and F ^ were defined earlier in this chapter. Equation 
(9-25) is also valid for the case without input filter, all 
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that is necessary to do is to set 3=0, while for the case 
without feedforward c 2 is set equal to 2 ero. Figure 95 shows 
the calculated values of open loop gain and phase margin. It 
is important to note at this point that equation (9-25) is 
an exact equation and that the assumptions made earlier in 
the feedforward design process are not made in deriving this 
equation. 

It is seen from Figure 95 (a) and (b) that the two 
stage input filter causes disturbances in the open loop gain 
and phase at the two resonant frequencies, but the addition 
of feedforward removes these disturbances. The feedforward 
loop is seen to be effective at both the resonant frequen- 
cies though at the higher frequency there is a drop in phase 
margin. The two assumptions made in the feedforward design 
process have not been used in calculating the open loop gain 
with feedforward, and thus it is seen that Figure 95 indi- 
cates that the feedforward is effective in eliminating input 
filter interaction with the regulator control loop. 

The same regulator circuit was used to obtain experi- 
mental results using a single stage input filter: 

Rj.jSsO.2 ohm L =116 microH Cl = 20 microF 
The regulator was set up to measure the open loop gain and 
phase margin [8J. The nonadaptive feedforward circuit of 
Figure 67 was used in making the measurements. The feedfor- 
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Figure 95: (a) Open loop gain: Calculated values without 

input filter, with (A) and without feedforward 
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Figure 95: (b) Open loop pi 

without input fill 
feedforward 


without feed 
^ forward 



margin: Calculated values 
and (A) and without 


ward gain was calculated from equation (9-22) as -0.009 and 
this led to the following choice of feedforward circuit par- 
ameters : 

C f = 27 microF R fl = 5.6 Kohm R f2 = 50 ohm 
The input of the feedforward circuit is the input filter ca- 
pacitor voltage, and the output is subtracted from the out- 
put at the integrater in the feedback loop, as was done ear- 
lier. 

Figure 96 shows some preliminary measurements obtained. 
It is noticed that the input filter causes a disturbance in 
open loop gain and phase margin at the resonant frequency 
but the feedforward compensates for the input filter inter- 
action to some extent. The approximation of equation (9-20) 
is valid at low frequencies, thus if the input filter reso- 
nates at high frequencies the feedforward as designed may 
not be effective in completely compensating for input filter 
interaction. Measurements of Figure 96 were obtained using 
an input filter with a somewhat high resonant frequency 
which can be used to explain the fact that the feedforward 
is effective to a somewhat lesser degree in conpensating for 
the input filter interaction. Due to lack of time extensive 
experimental verification of the feedforward design was not 
possible; more measurement data is necessary to check the 
effect of feedforward on open loop gain and phase, and other 
performance specifications. 
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Figure 96; 
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The feedforward design process for the boost regulator 
would be identical to the one presented above; similar as- 
sumptions may be necessary to accomodate the positive aero 
in the boost regulator transfer function. 


9.3 REMARKS 

The following observations are made: 

1. The feedforward design process presented in Chapter 5 
is general; it can be used for any type of duty cycle 
control and for single loop or multi- loop feedback 
control schemes. 

2. The problem of input filter interaction with the re- 
gulator control loop is present in the buck-boost re- 
gulator also and the concept of feedforward compensa- 
tion can be extended to design a feedforward loop 
that eliminates the interaction. 




Chapter X 

CONCLUSIONS AND FUTURE WORK 



P:; 

t . 




10.1 C ONCLUSIONS 

A novel scheme for input filter compensation is presented in 
this dissertation. The input filter, while a necessary com- 
ponent of a switching regulator, interacts with the regulator 
control loop and this interaction can cause serious perfor- 
mance degradation such as loop instability, degradation of 
transient response, audiosusceptibility and the output impe- 
dance characteristics [3,4,5,6,10]. These problems are 
caused mainly by the peaking of the output impedance of the 
input filter and its interaction with the control loop. 

Conventional input filter design techniques for single 
stage and two stage input filters [3,5,6] minimize the peak- 
ing effect but this often results in a penalty of weight or 
loss increase in the input filter. A novel feedforward com- 
pensation scheme is developed for buck regulators that eli- 
minates the undesirable interaction between the input filter 
and the regulator control loop. The compensation scheme is 
implemented by a feedforward loop that senses the input fil- 
ter state variables. Extensive analytical and experimental 
data confirm that the regulator is made immune to the peak- 
ing of the input filter output impedance. 


y 


D 


D 

D 
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The following points regarding the feedforward 
compensation scheme presented are noteworthy: 

1. The feedforward design process is straightforward and 
it is seen that the feedforward loop gain is indepen- 
dent of the type of feedback control used and of the 
feedback loop transfer functions. 

2. An important characteristic of the feedforward con- 
trol scheme presented is its ease in implementation. 

3 . Extensive experimental data supported by analytical 
results show that significant improvement in perfor- 
mance is achieved with the use of feedforward in the 
following categories — 

(a) open loop gain and phase margin; 

(b) audiosusceptibility; 

(c) output impedance; and 

(d) transient response. 

The data shows clearly that the peaking effect of the 
input filter output impedance is eliminated with the 
use of feedforward. 

4. It is shown that the feedforward loop is independent 
of the input filter parameters and also . independent 
of the input filter configuration. 

Examination of the extensive experimental data con- 
firm that the use of the feedforward compensation 


5 . 
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scheme does not have any detrimental effect on regu- 
lator performance. 

6. If a switching regulator is acquired as a 'black boss’ 
for use in a system, providing an inadequate input 
filter could lead to system instability [4,5]. This 
is shown analytically using a program that calculates 
the eigenvalues of the system, and confirmed experi- 
mentally. It is then shown that the use of the feed- 
forward compensation scheme can stabilize a system 
that was made unstable due to input filter interac- 
tion with the regulator control loop. Experimental 
verification of the stabilizing action of the feed- 
forward loop is presented and this shows that an ar- 
bitrary input filter may be used with the feedfor- 
ward- loop controlled regulator. 

7. From the results presented it is logical to conclude 
that the feedforward loop can provide effective com- 
pensation for an unknown source impedance. For exam- 
ple, a preregulator which often has an unknown, dy- 
namic output impedance can interact with a switching 
regulator downstream and result in system instabili- 
ty. The feedforward circuit effectively isolates the 
switching regulator from its source thus preventing 
interaction between the switching converter and 
equipment upstream. 
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8. The use of the feedforward compensation scheme re- 
moves some of the input f ? Iter design constraints 
thus making the input filter design process simpler 
and allows the input filter to be optimized. A high 
performance regulator system with an optimum input 
filter can thus be realized with the use of the feed- 
forward compensation scheme. 

9. The feedforward design concept presented can be easi- 
ly extended to types of control other than those used 
in this dissertation. 

10. Extension of, the concept of feedforward compensation 
to other types of switching regulators is possible. 
Such a scheme for a buck-boost regulator is presented 
and preliminary results indicate the validity of the 
compensation scheme. 

10.2 SUGGESTIONS FOR FUTURE WORK 

The following suggestions for future work are made: 

1. The concept of feedforward compensation for buck- 
boost switching regulators needs to be verified ex- 
perimentally . Data indicating the effect of using 
the feedforward scheme on performance categories like 
audiosusceptibility, output impedance, transient res- 
ponse and stability is also needed. 
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2 . An investigation into the use of the feedforward com- 
pensation scheme for the boost converter should be 
made . 

3. The use of the feedforward loop is shown to isolate 
the regulator from its source, and this could be used 
to advantage in situations where regulators are put 
in series/parallel for load sharing. 

4. Current injected control is an important mode of con- 
trol and an investigation of the feasibility of using 
feedforward compensation for this type of control 
could be made. 

5. Single loop control using the constant frequency duty 
cycle type control is very widely used and test re- 
sults using feedforward compensation for this type of 
control are needed. 

6. It may be possible to simplify the adaptive feedfor- 
ward circuit by using the dc component of the input 
filter capacitor voltage instead of the supply vol- 
tage. This may have an effect on the transient res- 
ponse for a step change in supply voltage since then 
the duty cycle would change in accordance wildi chang- 
es in input filter capacitor voltage. 


Appendix A 


SOFTWARE FOR 

DISCRETE TIME SIMULATION OF A BUCK REGULATOR 


C * DISCRETE TIME SIMULATION OF A BUCK CONVERTER. * 

c 

C THIS PROGRAM SIMULATES THE STEADY-STATE BEHAVIOR AS WELL AS 
C THE TRANSIENT PERFORMANCE FOR A BUCK REGULATOR 
C WITH INPUT FILTER, WITHOUT FEEDFORWARD. 

C 

DIMENSION X(6,9Q00) , Fl(6,6),F2(6,6),G1 (6,4) ,02(6,4) , TIME (9000), 
1U(4),PHIF(6,6},PHIN(6,6), DF(6,4),DN(6,4), TVEC1(6),TSWIT(2), 

2V I SWI T(2), V0( 9000 ),PHIFT(6.6),DFT(6,4), TVEC2( 9000) 

DIMENSION F3(6,6J,G3(6,4 ) ,TVEC3(6) » PH. FT3( 6,6) ,DFT3( 6,4 ) , 

1 PHI F3(6,6) , DF3(6,4 J 
C 

REAL L1,L,K1,K2, JQMAX 
C 

DATA RL1, LI, Cl. L, RO, RN, RC, C/1.0, 325. E-6.20.E-5, 230. E-6, 0.2, 0.63, 
20. 067, 300. £-6/ 

DATA RL,R1 1 ,R12,R13, R14,C2, Rl|,CPl/20.0, 33. 3E3, 16. 7£3, 2. E3, 

1 4? . E3 , 0 . 0 1 E-6 , 40 . 7E3 , 5600 . E- 1 2/ 

DATA ER, ET, VI ,£0/6.7,7.0,40. ,20./ 

DATA TON, TOFMIN,TSWIT,VISWIT,VITON/2. POE-5, 5. E-6, 0.003, 10. , 
125.0, 50., 0.88E-3/ 

DATA NIT, FRAC, E PS, TF/ 10, 0.5, 5. £*’6,0.01 2/ 

DATA IQMAX.EQ, ED/8.0,0.2,0. 7/ 

CL1S=0. 5291665 
VC1S=39. 48866 
ECS=7.0 

CLS=0. 09291243 
VCS=20. 05715 
ERS-20. 05712 
VOS= 19. 99638 
C 

TON=V I TON/VI 

K1=(R12/(R14 # (R1 1+Rl2)))+(1./R13)-( RN/R4) 

K1“K1-( RN#RO*( RL+RC} )/( R4*RL*RC) 
K2=(RN/R4)-(1./R13J-(R12/{R14*<R11+R12) )) 

DO 1 1=1.6 
DO 80 J=1 ,4 
Gif I * J |=0. 

G2( J,J)=0. 

DF( l,J)=0. 

DN{ I , J )=0. 
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80 DFT( I , J )=0. 

DO 1 J=1 , 6 
F1( l.J)=0, 

1 F2{ I , J )=0. 

C 

DO 2 (=1,6 
DO 2 J=1 , 6 
PHIFTf l,J)=0. 

PHI F( I , J )=0. 

2 PHIN(l,J)=0. 

C 

F1(l,1)=(-T. *RL1 )/Ll 

F1{1,2)=-1./L1 

F1{2,1)=l./C1 

F1(2,4)=-1./C1 

F1(3,2)=(-1.*RN)/(CP1*R4) 

FI ( 3, 4 )=( -1 . *RL ft RC*K1 )/( CPI *( RL+RC) ) 

Fl(3,5)=(fiL*K2)/{CP1#<RL+RC)) 

F1<3,6)=1./(CP1+R13) 

F1(4,2)=1./L 

F1(4,4)={-1.«RO)/L+(-1.*RC*RL)/(L w (RC+RL)) 
F1(4,5)=(-l. tt RL)/(L*(RC+RL)) 
F1(5,4)=RL/(C*(RL+RC) ) 

Fl(5, 5) =- l. /(C*{ RL+RC) ) 

FI { 6, 4 )-< RL*RC ) / ( C2*R 1 3*( RL+RC > ) 

FI ( 6, 5 )=RL/( C2*R1 3*( RL+RC) ) 

F1(6,6>=- 1./(C2#R13) 

G1(t,1 I-1-/L1 
G1(3,2)=1./(CP1*R14) 

G1 ( 3,3 )=RN/(CP1*R4) 

G1{4,3)=-T./L 
DO 3 1=1,6 
DO 3 J=1 ,6 

3 F2( I , J )=F1 ( I , J ) 

F2{2,4)=0. 

F2{ 3,2)=Q. 

F2(4,2)=0. 

G2{1,1)=G1(1,1) 

G2(3,2)=G1 (3,2) 

G2(3,4)=RN/(CP1*R4) 

G2{4,4)=-1./L 
DO 530 1=1,6 
DO 531J=1,6 
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F3< i,J)=F2( l,J) 

DO 532 J=1,4 
532 G3{ l,J)=G2( I , J ) 

530 CONTINUE 

00 535 1=1.6 
F3(4,l )=0.0 

535 F3(l,4)=0.0 
DO 536 1 = 1, J* 

536 G3(4,t)=G.O 
C 

NSW I T=0 
1TEMP=NSW1T+1 
TXX=TSWIT( I TEMP) 

NT=0. 1+1 . /FRAC 

TOFF=TON*(V»~EO)/EO 

TJME( 1 )=0.0 

ND=1 

M=0 

IT=0 

T=FRAC*70FF 

C 

CALL STRAN{T,PHIF,DF,F2,G2, 6, 4,6, 4, 1,0, EPS. DlFMAX, HER, (FLAG) 
C 

T=FRAC*TON 

GALL STRAN(T,PHiN,DN,F1.G1,6,4,6,4,1 ,0,EPS,Di FMAX, ITER, I FLAG) 
WRlTE(6,704) 

704 FORMAT(//,' VALUES OF X ON THE NEXT PAGE* ) 

C 

X(1,t)=CL1S 
X(2, 1 }=VC1S 
X( 3, 1 )=ECS 
X[4,1 )=CLS 
X(5,1 )=VCS 
X( 6, 1 )=ERS 
V0( 1 )=V0S 
U{1)=VI 
U{ 2 )=ER 
U( 3 )=EQ 
U(4)=ED 


1F(X<3,ND).GE.ET) GO TO 60 
TEMPA=ABS( X{ 3 , ND ) -ET ) 






t-itm! i* U "- ''j j rnCT h oai* | p— HI 


CT3 a 




ORIGINAL P&Q1E IIS 
OF POOR QUALITY 


uoooo 




mn urn n cn 


-.i 








1 F(TEMPA. LT. EPS) GO TO <10 
20 CONTINUE 

I F( X( 3, ND) . LT . ET ) GO TO 20 

MINIMUM OFF TIME CALCULATION * 

T=T0FMJ N 

CALL STRAN) T, PH I FT, OFT, F2,G2,6,4, 6,4, 1 ,0, EPS, 01 FMAX, ITER, I FLAG) 

00 6 1 = 1,6 

6 TVEC1( I )=X( l,ND) 

CALL FUTURE( TVEC1 ,U, PHI FT, DFT, 6 , 6 , 6 , 4, 6, 6 , 6 , 4, 1 ) 

TEMP=T1ME(ND) 

ND=ND+1 

TIME(ND)=TEMP+T 
DO 7 1=1,6 

7 X( I ,ND)=TVEC1( I ) 

VO { ND ) = ( RL*RC ) *X( 4 , N D ) / ( RL+RC ) +RL*X( 5 , HD ) / ( RL+RC } 

I F{X( 4, ND) .GT.0.0 ) GO TO 570 
I T 1=0 

571 I F(ADS(X{4, NDn.LT.EPS) GO TO 579 
ITI=ITI+1 

1 F{ NlT.GE.NI T) GO TO 100 

SLOPE3=F2(4,4)*X(4,ND)+F2(4,5)*X(5,ND)+G2(4,4)*U(4) 

T3=-X{ 4, ND ) /SL0PE3 

CALL STRAN(T3 , PH I FT3.DFT3 , F2, G2 , 6 , 4,6, 4, 1 , 0, EPS, D I FMAX, I TER, I FLAG) 
DO 572 1=1.6 

572 TVEC3) I }=X( I ,ND) 

CALL FUTURE ( TVEC3,D, PHI FT3, DFT3, 6,6, 6,4,6,6,6,4, 1 ) 

TIME(ND)=TIME(ND)+T3 

DO 573 1=1,6 

573 K{ l,ND)=TVEC3( I ) 

VO(ND)=RL*RC*X(4,ND)+RL*X{5,ND) 

VO(ND)=VO(ND}/( RL+RC) 

GO TO 571 
579 CONTINUE 

CALL STRAN( T3 , PHIF3,DF3»F3,G3,6,4,6,4,1,0,EPS,Dt FMAX, ITER, IFLAG) 

DO 574 1=1,6 

574 TVEC3( I )=X{ I ,ND) 

X( 4, ND)=0.0 

GALL FUTURE) TVEC3,U, PHIF3,DF3,6,6,6,4,6,6,6,4,1) 

TEMP=TIME(ND) 
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ND-ND+7 

TIME! ND )=TEMP+T3 
00 575 1=1,6 
575 X( I ,ND)=TVEC3( I ) 

V0(ND)={RL*X{5,ND))/{RL*RC) 

570 CONTINUE 
GO TO 40 
C 

C OFF TIME CALCULATION • * 

c 

28 T=FRAC*T0FF 
IT=0 

i 25 CONTINUE 

DO 8 1=1,6 

8 TVEC1 ( I )=X( I ,ND) 

IF(ND.LT.300) GO TO 5000 
WRITE! 6, 706) 

706 FORMAT!//, 1 STATE VECTOR IN OFF TIME 1 } 
WRlTE(6,9004){TVEC1( I ), 1=1,6} 

9004 FORMAT! 6fT5. 11) 

WRJTE(6,*) VO(ND) 

5000 CONTINUE 

iF!X(4,ND).LE.0) GO TO 502 
C 

v CALL FUTURE! TVEC1,U, PH IF, OF, 6, 6, 6, 4, 6, 6, 6, 4,1) 

I C 

■f. TEMP=T IME(ND) 

1 ND=ND+1 

> TIME! ND)=TEMP+T 

DO 9 1=1,6 

i 9 X{I,ND)=TVEC1( I ) 

; VO(ND)=RL*RC*X(4,ND)+RL*X{5,ND) 

‘ VO(ND)=VO(ND)/(RL+RC) 

' • GO TO 520 

! 502 CONTINUE 

C 

C * CALCULATION TO FIND WHEN X[4)=0 * 

! ■ c 

IT 1=0 
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503 I F{ABS(X{ 4, ND) J.LT.EPS) GO TO 509 
ITI=ITI+1 

1 F( ITI.GE.NJT) GO TO 100 

SLOPE3=F2(4,4)*X[4,ND)-i-F2(4,5)*X(5,ND)+G2{4,4)*U(4) 

T3=-X{4,ND)/SL0PE3 

CALL STRAN f T3 , Pti l FT3 , DFT3 , F2 , G2, 6 , 4 . 6 ,4 , 1 , 0 , EPS, 0 1 FMAX, 1 TER, t FLAG ) 
DO 504 1=1,6 

504 TVEC3( 1 )=X( I , ND) 

CALL FUTURE{TVEC3, U, Pill FT3,DFT3, 6,6, 6, 4,6,6, 6,4,1) 

TIME| ND)=TlMEvHD)+T3 - 

DO 505 1=1,6 

505 X{ l,ND)=TVEC3( I) 

V0{ND)=RL*RC*X(4,ND)+RL*X(5,ND) 

VO ( ND )=V0 ( HD) /t RL+RC ) 

GO TO 503 
509 CONTINUE 

CALL STRAN(T, PH I F3,DF3, F3,G3,6,4,6,4, 1,0, EPS, D I FMAX, ITER, I FLAG) 

DO 5401=1,6 
540 TVEC3( l)=X{ I , ND) 

X(4,ND)=0.D 

CALL FUTURE( TVEC3, 0, PH l F3, DF3 , 6, 6, 6, 4, 6, 6, 6,4, 1 ) 

TEMP=T IME( ND) 

ND=ND+1 

TIME(ND)=TEMP+T 
DO 511 1=1,6 
5H X( l,ND)=TVEC3( I ) 

V0( ND)=( RL*X( 5, ND) )/( RL+RC) 

520 CONTINUE 

IF(X(3,ND).GE.ET) GO TO 30 
GO TO 25 
C 

G CALCULATION TO HIT THRESHOLD WHEN * 

C X{3)=ET. * 

c 

30 I F( ABS(X(3, ND)~ET) . LT. EPS) GO TO 40 

IT= I T+1 

IF( IT.GE.NIT) GO TO 100 

SL0PE=F2{3,4)*X{4,ND)+F2{3,5PX{5,ND}+F2t3,6)ttX(6,ND} 

SL0PE=SLOPE+G2(3,2)«U(2) 

T={ ET-X(3, ND) ) /SLOPE 
C 
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CALL STRAN{T, PH I FT, OFT, F2,G2,£,4,6,4, 1,0, EPS, D1 FMAX, ITER, I FLAG) 
C 

DO 10 1=1,6 

10 TVEC1 ( I )“X( I , ND) 

CALL FUTURE (TVEC1,U, PHI FT, DFT, 6,6, 6,4, 6,6,6, 4,1) 

-TtMEf ND)=TIME{ND)+T 
DO 11 1=1,6 

11 X{ l,ND)=TVEG1( I ) 

V0( ND )=RL*RC*X( 4, ND)*RL*X! 5, ND ) 

VO C ND J =V0 ( ND ) / ( RL+RC ) 

GO TO 30 

: 40 CONTINUE . . 

IF1ND.LT. 300) GO TO 5001 
WRITE! 6, 440) TIHE(ND) 

440 FORMAT!/, 'SWITCH ON TIME=*, F15.11 ) 

5001 CONTINUE 

JF(TIME(ND).GE.TF) GO TO 120 
C 

C CHANGE VI, TON, TOFF AT THE SWITCHING 
G INSTANT. ALSO CHANGE PH IF, DF, IF * 

C NECESSARY. , * 

c 

I F(TIME! ND). LT.TXX) GO TO 60 
Mswrr=NswiT+i 
ITEMP=HSWIT+1 
TXX=TSWIT( l TEMP) 

V I =VI SWI T( NSW IT ) 

TON=VITON/VI . 

T0I r F=T0N*! VI -EO)/EO 

u(i)=vi 

T=FRAC«TDFF 

C 

CALL STRAN! T, PH I F, OF, F2, G2, 6,4, 6,4, 1 ,0, EPS, D I FMAX, I TER, I FLAG) 

C 

T=FRAC«TOH 

C 

CALL STRAN(T, PHIN,DN, FI, G1, 6, 4, 6, 4, 1,0, EPS, D I FMAX, ITER, I FLAG) 

C 

60 CONTINUE 
M=0 

T=FRAC*TON 

C 
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ON TIME CALCULATION * 

65 CONTINUE 

DO 12 1=1 >6 

12 TVEClf I)=X( j,ND) 

IFIND.LT. 300} GO TO 5002 

: WRITE{6,709) ’ 

709 FORMAT ( //, 1 STATE VECTOR IN ON TIME 1 ) 

WRiTE(6,9005}(TVEC1( I), 1=1,6) 

9005 FORMAT( 6F15. 11 ) 

WRITE(6,«) VO(ND) 

5002 CONTINUE 
C 

call FUTURE(TV£C1,U, 2111^04,6,6,6,4,6,6,6,4,1 ) 

c 

TEMP=TIME(ND) 

ND=ND+1 

T I ME( MD )=TEMP+T 
DO 13 1=1,6 

13 X( I,ND)=TVEC1( I ) 

VQ(ND)=RL^RG lt X(J| J ND)'FRL«'X{5,ND) 

VO(ND)=VO(ND)/{RL+RC) 

M=M+1 

IF(X(4,ND).GE.IQMAX) GO TO 70 
IF(M.LT.Nr) GO TO 65 

G .. ■ 

GO TO 999 
C ' 

C ■ -tfrtH «««»«- 

C CALCULATION TO FIND WHEN INDUCTOR * 

G CURRENT HITS THE THRESHOLD VALUE. * 

c 

70 CONTINUE 
IT=0 

71 I F(ABS(X[N, ND)-IQMAX) . LE. EPS) GO TO 999 

IT=!T+1 

1 F( IT. GE. NIT) GO TO 100 

SL0PE2=X(2,ND)*F1{4,2)+X[4,ND)*n(4,4)+X{5,ND)*F1(4,5) 

T=( I QMAX-X( 4, ND ) )/SL0PE2 

CALL STRAN(T, PHI FT,DFT, Fl ,G1 ,6,4,6,4,1 ,0,EPS, DI FMAX, ITER, IFLAG) 
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15 

16 

999 

C 

5003 

TOO 

204 

120 

103 

104 

105 

106 

107 

14 

C 

c 
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DO 15 1=1,6 
TVEC1 ( I }=X{ t,ND) 

CALL FUTURE{TVEC1, U, PHI FT, DFT, 6,6, 6, 4, 6, 6, 6, 4, 1 ) 
T IHE( ND)=T I ME( NDJ+T 
DO 16 1=1,6 
X( I,ND)=TV£C1( 1 ) 

VO ( UO ) =RL#RC*X t 4, ND ) + RL*X( 5 , ND ) 
VO(ND)=VO(ND)/(RL+RC) 

GO TO 71 
CONTINUE 


IF(ND.LT,300) GO TO 5003 
VmiTE{6,4l!2j TIME(ND) 

FORMAT (/>' SWITCH OFF T1ME=’ , F15. 11 J 
CONTINUE 

IF{T1HE(ND).LT.TF) GO TO 20 

GO TO 120 

CONTINUE 

WRITE (6,204) TIME(ND) 

FORMATE/, 'MAXIMUM ITERATION AT TIME=' , F15. 11 ) 
CONTINUE 
DO 104 1=1, NO : 



CONTINUE 
DO 105 1=1, ND 

WR I TE( 8, 103) T I ME ( 1 ),X(4, 1 ) 

CONTINUE 

DO 106 1=1, ND . 

V/RITEf 8, 103 ) TIHE( l),VOf I ) 
CONTINUE. 

DO 107 1=1, ND 

V1RITE(8,103) TIKEE I ),X{2, I } 

CONTINUE 

DO 14 l=1,ND 

TVEC2( I )=VOf f ) 

CALL ABPL0T(ND,TIME,TVEC2) 
STOP 
END 
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SUBROUT I NE STRAN{ TAU, PH I , THETA, A, B, NA, NB, MA, MB, MODE, NTERMS, 
TOL, D I FMAX, ITER, I FLAG ) 
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REAL PHI (NA, NA), THETA(NA, NB),A( NA,NA),B{NA, MS}, 

1 WORK! ( 1 0, 1 0 ) , WORI<2{10, 10 ) , DUMMY! 1,1), FAC1 , CONI , FAC2, COM2 
C1=1 
TS=TAU 
DO 4 l=T,MA 
DO 2 J=1,MA 
WORK1 { I , J)=0 
W0RK2(l,J)=O 
PHI ( 1 ,J)=Q 
2 CONTINUE 
W0RK1 (1,1 )=1 
W0RK2C I, l)=TAU 

piun, d=i 
4 continue 
DO 6 1=1, MA 
DO 6 J=1,MB 
6 THETA( l,j)=0 
DirMAX=1.E8 
ND0=50 

IF(NTERMS.GT.O) NDO=NTERMS 
FAC1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMUL( DUMMY, WORK! , A, TS, 7,1,10, 1 0, MA, NA, MA,MA,MA, MA. 1 } 

FAC1=FAC1*I 

CON1=1./FAC1 

IF( NTERMS. EQ.O .AND. 1 . GE.4 ) CALL SERRORf PH l , WORK! , CONI, HA, NA, 

1 10,10,MA,MA,DIFMAK) 

CALL MXADD( DUMMY, PHI , K0RK1 , Cl , CONI ,1 , 1,NA,NA, TO, 10, MA,MA, 1 ) 

1F{M0DE. EQ.2) GO TO 500 

FAC2=FAC1*(1+1) 

C0N2=TAU/FAC2 

CALL HXADDf DUMMY, W0RK2, WORK T, Cl ,C0N2, 1 , 1 , 1 0, 10, 10, 10, MA, MA, 1 ) 
500 CONTINUE 
JTER=I 

I F( NTERMS, GT. 0 . OR. J . LT. 4 ) GO TO 10D0 
IF(DIFMAX.LE.TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

1 F( MODE. Eq, 1 ) CALL MXMUL (THETA, WOR K2 , B, Cl , NA, NB, 10,10, NA, NB, 

1 MA, MA, MA, MB, 2 ) 

I FC ITER. Eq.NDO.AND. NTERMS. EQ.O) I FLAG=1 
RETURN 
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END 

C 

G . 

1 SUBROUTINE SERROR(AMX,BMX,CCC, IA,JA, IB,JB, I DO, JDO, D I FMAX) 
DIMENSION AMX( IA,JA),BMX( IB,JB) 

Dl FHAX=1 . E-30 
DO 100 1=1, IDO 
DO 50 J=1,JD0 

I F( AMX( I , J J . EQ.0.0 ) GO TO 50 
CHANGE^ ABS{ BMX( I , J )*CCC/AMX( 1 ,J ) J 
I F( CHANGE. GT.DIFl-IAX) DI FMAX=GHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 
C 

c : 

1 SUBROUTINE FUTUREJX, V, PHI , THETA, LP,MP, LT,HT, 1 P, JP, IT,JT,M0DE) 
DIMENSION X{ MP) , V( MT J , PH I ( LP, MP ) , THETA{ LT , MT) , TEMP( 20 J 

DO 20 1=1, IP 
SUM=0 

DO 10 J=1,JP 
1 0 SUM=S r JM+PH I ( I , J )*X( J ) 

20 TEMP( l }=SUM 
DO 30 1=1, IP 
30 X( I )=TEMP( I ) 

IF(MQDE.Eq.2) RETURN 
DO 60 1=1, IT 
SUM=0 

DO 50 J=1,JT 

50 SUM=SUM+THETA{ I,J)*V(J) 

X{ I )=X£ I J-i-SUM 
60 CONTINUE 
RETURN 
END 


1 SUBROUTINE MXADD(RMX,AteX,BMX,ACC,BCC, IR,JR, IA,JA, IB,JB, ID0,JD0, 
1 MODE) 

DIMENSION' AMX£ IA,JA),BMX£ IB,JB),RMX( IR,JR1 

J F£ IA. LT, 1 DO.OR.JA. LT. JDO) GO TO 999 
IF( IB.LT. IDQ.OR.JB.LT.JDO) GO TO 999 
GO TO (10, 100), MODE 
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10 CONTINUE 

DO 50 1=1, IDO 
DO 50 J=1 , JDO 

50 AMX( t,J)=AHX( I , J)*ACC*BMX{ 1,J)*BCC 
GO TO 300 
100 CONTINUE 

I F( IR.LT, 100. OR.JR.LT. JDO) GO TO 9 99 
DO 20G 1=1, IDO 
DO 200 J=1 , JDO 

200 RMX( I, J)=AMX( I , J )#ACC-rBMX( I , J)*BCC 
300 RETURN 
99 9 CONTINUE 
RETURN 
END 


I^SUBROUTJNc. MXHUL(RMX,AMX,BMX,CCC, IR,JR, !A,JA, !B,JB, IDA, J DA, IDB,JDB 

'DIMENSION AMX( IA,JA),BMX( IB,JB),RMX( IR,JR),TEMP{20) 

1 F{ 1DA*IDB^DA*JDB,GT. 1A*IB*JA*JB) GO TO 999 
IF(JDB.GT.JDA) GO TO 999 
GO TO (10, 210), MODE 
10 "DO 100 1=1, IDA 
DO 20 L=1,J0A 
20 TEMP(L)=AMX( i,L) 

DO BO J=1 , JDB 
SUM=0 


DO NO |(=1,JDA 

NO SUM= SUM+TEM P { K ) *BMX( K, J ) 

AHX{ J,J)=SUM*CGC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 
210 CONTINUE 

IF( IR.LT. IDA. OR.JR.LT, JDB) GO TO 999 
DO i|00 1 = 1, IDA 
DO 300 J=1, JOB 


SUM=0 

DO 3N0 K=1 , JDA 

3R0 SUM=SUM+AMX( 1 ,K)*BMX( K, J ) 
RMX( i,J)=SUM*CCC 
380 CONTINUE 
NOO CONTINUE 
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600 RETURN 
999 CONTINUE 
RETURN 
END 
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C DISCRETE TIME simulation of a buck converter * 

C THIS PROGRAM SIMULATES THE STEADY-STATE BEHAVIOR AS WELL AS 
C THE TRANSIENT PERFORMANCE FOR A BUCK REGULATOR 
C WITH INPUT FILTER AND WITH FEEDFORWARD. 

C 

DIMENSION X(7,9000), FI (7, 7), F2(7,7},G1 (7,4) ,G2{7,4>, T!ME( 9000), 
1U[1|), PHIF(7,7),PH1N(7,7),DF{7,4J,DN{7,4),TVEC1(7),TSWIT(2), 
2ViSWIT(2},VO(90D0),PIMFT(7,7),DFT(7,4),TVEC2(90(Hn 
D‘ MENS JON F3(7,7},G3(7,4j,TVEC3{7),PHIFT3(7,7),DFT3(7,4), 

1PH I F3{7,7),DF3(7, 7) 

C 

REAL L1,L, K1 ,K2, IQMAX.K3 
C 

DATA RLl,L1,C1,L,R0,RN,RC,C/1.0,325.£-6,22.E-5,230.E-6,0.2,O.65, 
20.067 P 30G.E“6/ 

DATA RL,R 11, R12,R13,R14,C2, R4, CPI/20. 0,33. 3E3, 16. 7E3,2.E3, 

1 47 . E3 , 0 . 0 1 E-6, 40 . 7E3 , 5600. E- 1 2/ 

DATA ER , ET, V I , EO/6 .7,7.0,40. ,20./ 

DATA TON, TOFMIN, TSWIT, VISWIT, VITON/2.20E-5>5. E-6, 0.003, 10.0, 

125. ,50.,0.88E-3/ 

DATA N I T , FRAC , EPS, TF/ 10, 0 . 1 , 1 0 . E-6 , 0 . 01 20/ 

DATA I QMAX, EQ, ED/6 . 0 , 0 .2, 0. 7/ 

DATA RF1 ,RF2,CF,VFS/5.1E3,90.970,27.E-6,0.0507E-1/ 

CL1S=0. 69552390 
VC1 3=39.87245 
ECS=7.0 

GLS=0.40158130 
VCS-20.06141 
ERS=20. 06143 
VOS=20. 02125 
VFS=+0. 38948170 
C 

TON=VITON/VI 

K1 = (R12/(R14*(R1 1+R12) ))•*•( 1 ,/R13 )-(RN/R4) 

KI-KI-tRN^RO^JRL+RC) )/(R4*RL*RC) 

K2=tRN/R4)-[ 1 ./R13 )-(R12/(RT4 tf (R11+R12) ) ) 
l(3=RF2/{ RF1+RF2 J 
DO 1 1-1,7 
DO 80 J=1,4 
G1 ( I , J J= 0 . 

G 2 { l,J)= 0 . 


LO 

W 


o 




2 

G 



DF(l,J)=0. 
DNf I >J 3-0. 
DFT£I,J)=0. 
00 1 > 1,7 
Fl( I ,J)=0. 
F2( I > J)=0. 


DO 2 1=1,7 
DO 2 J=1,7 
PHIFT(1,J)=0. 
PHI F[ l,j'}=0. 
PHlN(l,Jj=0. 


FI ( 1,1 )=£"*1„ # RL1 )/Ll 

F1[1,2J*-»i*/L1 

F1(2,1)=1./C1 

FI £2, 3)=-1 „/(Cl*(RF1+RF2) ) 

F1(2,5}=-1./C1 

FU3,1J=1./G1 

FI { 3,3)=(-1 .*(CF+C1 ) )/(CF*C1*(RF1+RF2) } 

FI { 3, 3 )=0 

F1(3,5J=-1./C1 

F1£4,1}=£-1-*K3)/C1 

Flf4,2)={-*1.*RN)/(CPl*R4i 

F1{4,3)=(K3*£CF+C1))/(CF*C1*(RF1+RF2)) 

F 1 ( 4 , 5 )= ( K3 /Cl )- { £ RL*RC# K1 ) / { C PI * £ RL+RC ) } ) 
FT(4,6]=(RL*K2)/{CP1*(RL+RC)) 

FI ( 4,71=1 - /£CP1*R13) 

F1(!>,2J=1./L 

F1(5,5)=({“1. w R0J/L)+(£-1.«RC*RL)/{L*(RC+RL))) 

FT(5,6)=(-1.m)/£L«(RC+RL)) 

FI ( 6, 5 }- RF/ { C w f RL+RC ) ) 

Ft £ 6,6)=** 1 ./fC#f RL+RC) ) 
F1£7,5)=(RL*RC)/(C2*R13*{RL+RC)) 

Ft £ 7,6 J=RL/(C2*R1 3*£ RL+RC) ) 
F1{7,7)=-W£C2*R13) 

Cl £ 1> 1 )=1 ./LI 
Gt(4,2)=1 ./fCPl^Rl^) 

G1 (4,3)=RN/(CP1 i *R4J 
G1(5,31=*tT./L 
DO 3 1=1,7 
DO 3 J— 1 ,7 
F2(),J)=F1( F,J) 
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F2(2,51=0. 

F2(3,5)=0. . 

F2(3,3)=0 

F2(4 t 2)=0. 

F2 ( 4 , 5 ) = £ ~1 , *RL*RC*K1 ) / ( CPI * ( RL+RC ) } 

F2(5,2}=0. 

62(1, 1}=G1 (1,1) 

G2(4,2)=G1(4,2} 

G2(4,4)=RN/{CPV»R4) 

G2(5,4}=“1./L 
DO 530 1=1,7 
DO 531 J=1 ,7 

531 F3{ I , J)=F2( I , J ) 

DO 532 J=1 ,4 

532 G3{I,JJ=G2( i,J) 

530 CONTINUE 

DO 535 1=1,7 
F3(5,l 1=0.0 

535 F3{ 1,51=0.0 
DO 536 1=1,4 

536 G3(5, 1 )=0.0 
C 

NSWIT=0 

TEMP=NSWIT+1 . 

TXX=TSV/1T{TEMP} 

NT=0. 1+1. /FRAG 
TOFF=TON w (Vt-EO)/EO 
TlHEtl )=0.0 
ND=1 
*1=0 
1 T=0 

T=FRAC*TOFF 

C 

CALL STRAN(T, PHI F,DF, F2,G2,7,*4,7,4, 1,0jEPS, DJFHAX, ITER, 1 FLAG) 
C 

T= FRACTION 

CALL STRAN( T, PH 1 N, DN, FI , G1 , 7, 4, 7,4, 1 , 0, EPS, D I FHAX, 1 TER, I FLAG) 
WR!TE( 6,704) 

704 FORMAT{ //, 1 VALUES OF X ON THE NEXT PAGE 1 ) 

C 

X(t,1)=CL1S 

X(2,1)=VC1S 

X(3,1}=VFS 
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X{4,1)=ECS 
X(5,1}=CLS 
X(6,1)=VCS 
X(7,1 J-ERS 
V0( 1 }=VQS 
UU)=Vi 
' Uf2)=ER 
U(3)=EQ 
U(4)=ED 
C 
C 

lF(X(4,ND).GE.ETJ GO TO 60 
TEMPA=ABS{X(4,ND)-£T) 

I F( TEWPA.LT. EPS) GO TO 40 
20 CONTINUE 
C 

1F(X(4,ND).LT.ET) GO TO 28 
C 

C MINIMUM OFF TIME CALCULATION * 

C 

T=TOFMIN 

CALL STRANt T, PHI FT, DFT, F2,G2, 7,4,7»4, 1,0, EPS, 01 FMAX, ITER, J FLAG) 

DO 6 1=1,7 

6 TVEC1 ( 1 )=X{ I , NO) 

CALL FUTURE(TVEC1 ,U, PHI FT, DFT, 7, 7, 7, 4, 7, 7, 7, 4, 1 ) 

TEMP=T!ME(ND) 

: ND=ND+1 

TIME[ND)=TEMP+T 
DO 7 1=1,7 

7 X[ |,ND)=TVEC1{ I ) 

VO£ND)=(RL*RC)*X{5>ND)/(RL*RC)+RL#X(6,ND)/{RL+RC) 

I F(X(5, ND) .GT .0,0) GO TO 570 
I T I =0 

571 I F{ABS(X{5,ND) ) .LT.EPSJ GO TO 579 
ITI=iTI+1 

!F( IT! . GE. NIT) GO TO 579 

SL0PE3=F2(5,5)*X(5,NDJ+F2(5,6)*X{6,ND)+G2(5,4)*U(4) 

T3=-X(5,ND)/SL0PE3 

CALL STRAN(T3, PH! FT3,DFT3, F2,G2,7,4,7,4, 1,0, EPS,DI FMAX, ITER, IFLAG) 
DO 572 1=1,7 

572 TVEC3{ I )=X( I/ND) 
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CALL FUTUREJ TVEC3, U, PH t FT3, DFT3,7,7,7, 4, 7,7, 7,4, 1 ) 

TIME( ND)=TIM£l ND)+T3 
DO 573 1=1,7 

573 X(l,ND)=TVEC3( I) 

VO(H0)=RL*RC«X(5,NDJ+RL*X(6,ND) 

VO(ND)=VO{ND)/{RL^RC) 

GO TO 571 
579 CONTINUE 

CALL STRAN(T3, PHI F3,DF3,F3,G3, 7,4, 7,4, 1,0, EPS, D1FMAX, ITER, I FLAG) 
DO 574 1=1,7 

574 TVEC3U }=X( l,ND) 

X(5,ND)=0.0 

CALL FUTURE ( TVEC3 , U, PH I F3, DF3, 7, 7, 7, 4, 7, 7, 7, 4, 1 ) 

TEHP=TIME£ND) 

ND-ND-M 

T1ME(ND)=TEMP+T3 
DO 575 1=1,7 

575 Xn,ND)=TVEC3£ I } 

. VO(ND)=(RL*X(6,W)n/(RL+RC) 

570 GOUT l NUE 
GO TO 40 

OFF TIME CALCULATION * 
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T=FRACttT0FF 
I T— 0 

CONTINUE 



DO 8 1=1,7 

8 TVECl { I )=X{ I ,ND) 

I F£ND, LT.300) GO TO 5000 


WRITE(6,706) 

70fi FORMAK//,' STATE VECTOR IN OFF TIME') 
WR1TE(6,9004>(TVEC1(IJ, 1=1,7) 

9004 FORMAK 7F15. Iti) 

WR!TE{6,*) VO (NO) 

5000 CONTINUE 


! F{X{ 5, ND }.«_£. D.0) GO TO 502 


CALL FUTURE( TVECl, U, PHI F,DF,7, 7,7,4, 7,7, 7,4,1) 

c . 

TEMP-TIME(ND) 



oooon 




9 

502 


ND-NO+1 

TfHE(NDJ=TEMP-5-T 
DO 9 1=1 , 7 
X( l,ND)=TV£GU I) 

V0{ ND )=RL#RC*X{ 5, ND)-*-RL tt X{ 6, ND ) 
V0(ND)=V0(ND)/(RL+RG) 

GO TO 520 
CONTINUE 


* CALCULATION TO FIND WHEN X(5)=0 * 


I T I =0 

503 J F[ ABS(X( 5, NO) ) . LT. EPS) GO TO 509 

m=m+i 

I F ( I T I , GE . N ! T ) GO TO 100 

SL0PE3=F2(5,5)*X(5,ND)+F2(5,6)*X(6,ND)+G2{5,4)*U{4) 

T3=-X(5, NDJ/SLQPE3 

CALL STRAN ( T3 , PH I FT3 , DFT3, F2 , G2, 7, 4, 7 , 4, 1 , 0, E PS , D I FMAX, I TER, I FLAG J 
DO 504 1=1,7 

504 TVE03J I J=X( !,ND) 

CALL FUTUR£{TVEC3,U, PHI FT3,DFT3,7,7,7,4,7,7,7>4,1 ) 

TIHE(ND)=TIHE(ND)+T3 

DO 505 1=1,7 

505 X( I,ND)=TVEC3{ I ) . • 

VO( ND j=RL*RC*X{5, ND )+RL*X( 6, NO) 

V0( ND )=V0( ND)/( RL+RC) 

GO TO 503 
509 CONTINUE 

CALL STRAN{ T, PH I F3, DF3, F3, G3, 7, 4, 7, 4,1 , 0, EPS, D I FMAX, ITER, 1FLAG) 

DO 540 1=1,7 
540 TVEC3( I )=X( l,ND) 

X(5,ND)=0.0 

CALL FUTURE(TVEC3,U,PHIF3,DF3, 7,7, 7,4, 7,7,7, 4,1) 

TEMP=TIME(ND) 

ND=ND+1 

TIME(ND)=TEMP*T 
DO 511 1=1,7 
511 X{ |,ND)=TVEC3( I ) 

VO ( ND )= ( RL*X( 6 , ND ) ) / ( RL4RC) 

520 CONTINUE 

I F(X(4,ND) ,GE. ET) GO TO 30 
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C 


GO TO 25 

CALCULATION TO NIT THRESHOLD WHEN * 

X( 3 )=ET. * 

30 IF{ABS(X(4, NDJ-ET) .LT.EPS) GO TO HO 
1T=IT+1 

I Ff IT.GE.N1T) GO TO 100 

SL0PE=F2{4,1)*X(1,ND)+F2(4,3)*X(3.H0, l *F2(4,5)*X{5,ND)+ 
1F2£I|,6)*X(6,ND)+F2(4,7)*X(7,ND) 
SL0PE=SL0PE+G2{4,2)«U(2}+G2{4,4)*U(4) 

T=( ET“X{ 4, NO) )/SLQPE 

CALL STRAN f T, PH I FT, DFT, F2, G2, 7, 4, 7, 4, 1 , 0, EPS, D I FMAX, I TER, I FLAG) 


0 

C 


DO 10 1=1,7 
10 TVECrt J)=X( l,ND) 

CALL FUTURE(TVEC1 , U, PHI FT, DFT, 7,7, 7,4, 7, 7,7,4, 1 ) 
TfHE{ND)-T IHE{ND)FT 
DO I) 1=1,7 
■It X( J,ND)=TVEC1 ( 1 ) 

V0( ND )=RL^C*X( 5, ND) *RL*X( 6, ND) 
VO(ND}=VO[ND)/(RLTRC) 

GO TO 30 
40 CONT I HUE 

IF(ND.LT. TOGO) GO TO 5001 
WRJTEf 6,440) T1ME{ND) 

440 FORMAT!/, ’SWITCH ON TIME=* , FT5. 11) 

5001 CONTINUE 

IFCT1ME{ND).GE.TF) GO TO 120 

C ... 

C CHANGE VI, TON, TOFF AT THE SWITCHING * 

c instant. ALSO CHANGE PHI F, OF if *■ 

C NECESSARY. * 

c 


IF(T!ME(ND).LT.TXX) GO TO 60 
NSWIT=NSWIT+1 
TEMP=NSWiT+1 
TXX=TSWIT(TEMP) 
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VJ =V JSVf I T ( NSW ! T ) 

TQN=V I TON/VI 
TDFF=TON*!Vl-EO)/EO 
Uf 1 )=V1 
RF2=230.27 
K3=RF2/{ RF1+RF2 ) 

Ft (2j 3 )=~1 ,/{Cl*(RF1+RF2) ) 

FI (3>3 )={-1 .*{CF+G1 ) )/[CF # C1*( RF1+RF2) J 
Fl(4,1)=(-1.*K3)/Gt 

FI ( 4, 3 )=( K3*( CF+C1 ) )/( CF*C1*{ RF1+RF2 ) ) 

F1(4,5)=(K3/C1)-{(RL«RC 4 W)/(CP1*(RL+RG))) 

F2{2,3)=F1!2,3) 

F2( 3,3 )=F1 (3,3) 

F2{4, 1 3— FI (4,1) 

F£(4,3)=F1{4,3} 

T=FRAC#TOFF 

C 

GALL STRAN( T, PH I F, DF, F2, G2, 7, 4, 7, 4, 1 , 0, EPS, D I FMAX, ITER, I FLAG) 

C 

T=FRAC*TON 

G CO 

CALL STRAN (T, PHIN,DN, FI ,G1 ,7, 4, 7, 4, 1 ,0, EPS, 01 FMAX, I TER, I FLAG ) £ 

C vo 

60 CONTINUE 
M=0 

T=FRAC (t TON 


ON TIME CALCULATION * 


65 CONTINUE 

DO 12 1=1,7 
12 TVEG1 ( I )=X( 1 , NO) 

IF(ND.LT.IOOO) GO TO 5002 
WRITE! 6, 709) 

709 FORMAT (//,’ STATE VECTOR IN ON TIME 1 ) 
WR I TE( 6,9005 ) ( TVEC1 ( I }, !=’,7) 

9005 FORMAT! 7 FI 5. 11 ) 

WRITE( 6,*) VO(ND) 

5002 CONTINUE 


CALL FUTURE! TVEC1 , U, PR I N, DN, 7, 7, 7,4, 7, 7, 7,4, 1 } 
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rm 


cm cm cm cm cm 


EC 


TEMP=TIME(ND) 

ND=ND*1 

T!ME(ND)=TEMP+T 
DO 13 1=1,7 
13 X( J,NDJ=TVEC1( I ] 

V0(ND)=RL*RC*X(5,ND)+RL»X(6,ND) 

VO(ND)=VO(ND)/( RL+RC) 

M=M+1 

|F(X{5,ND).GE, iqMAX) GO TO 70 
IF(M.LT.NT) GO TO 65 

c • 

GO TO 999 

C . )i n u lt „ ] t „ !t u ic i i n it )i : t i t i; n i t i t tr ,ii J; ; i i i u Jt n i t l . fJ . 

C CALCULATION TO FIND WHEN INDUCTOR * 

C CURRENT HITS THE THRESHOLD VALHE. * 

c 

70 CONTINUE 

I T=0 ■ 

71 I F( ABS(X{5>ND)- IQMAX) , LE.EPSJ GO TO 999 
!T=|T+1 

IF(IT.GE.NIT) GO TO 100 

SL0PE2=F1(5,2)*X(2,ND)+F1(5,5)*Xf5,ND)*n(5,6}*X(6,ND)+ 

1G1{5>3)^U(3) 

T=( IQHAX-X(5,ND))/SLOPE2 

CALL STRAN [ T, PH I FT, DFT, FI , G1 , 7, 4, 7, 4, 1 , 0, EPS, D i FMAX, I TER, I FLAG) 
DO 15 1=1,7 

15 TVECi ( I )— X{ I >ND) 

CALL FUTURES TVEC1 , U, PH ! FT, DFT, 7, 7, 7, 4, 7, 7, 7, 4, 1 ) 

T I ME( ND )=T I ME( ND )+T 

DO 16 1=1,7 

16 X{ I,ND)=TVEC1( I) 

V0(ND)=RL*RC*X(5,ND]+RL*X(6,ND) 

VO(ND)=VO(ND)/[RL+RC) 

GO TO 71 
999 CONTINUE 

G ‘ 

IF(ND.LT.IOOO) GO TO 5003 
1VRITE(6,442) TIME(ND) 

442 F0RMAT(/, * SWI TCH OFF TIME= ,; , F15. IT) 

5003 CONTINUE 

I F(TIME( ND), LT.TF) GO TO 20 - 
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GO TO 120 
100 CONTINUE 

WR[TE(6,204) TIME(ND) 

204 FORMAT! /,* MAXI HUH ITERATION AT TIME=',F15.n ) 
120 CONTINUE 

111 DO 105 1=1, ND 
YmiTE(8,103) TIHEJ I ),X( 1,1) 

105 CONTINUE 

DO 106 1=1, ND 

HR I TEC S, 103) TIME(l),X{2, !} 

106 CONTINUE 

112 DO 104 1=1, NO 

WRITE! 8, 103) TIHE( I),X(3, t) 

103 FORMAT (2E20. 10) 

104 CONTINUE 

DO 107 1=7, NO 

WIUTE£8,103) TIME( l),X(ft, I) 

107 CONTINUE 

DO 122 1=1, ND 

WRITE! 8, 103) TIME! I ),X(5, I ) 

122 CONTINUE 

DO 123 1=1, ND 
WRITE[3,103) TIME! I ),V0{ I ) 

123 CONTINUE 

DO 124 1=1, ND 

WRITE(8,103) T1ME( I ),X(7, I ) 

124 CONTINUE 
109 CONTINUE 

DO 14 1=1, ND 
14 TVEC2( I )=VO( I ) 

STOP 
END 


SUBROUT I ME STRAN( TAU, PH l , THETA, A,B,NA, NB,MA,MB , MODE, NTERMS, 
1 TOL, DI FMAX, ITER, I FLAG) 

REAL PHI ( NA,NA$, TI!ETA( NA, NB) , A{ NA, NA) , B( NA, NB ) , 

1 V/ORK1 ( 1 0, 10 > , W0RK2( 10,10), DUMMY! 1,1), FACT , CONI , FAC2, G0N2 
01=1 
TS=TAU 
DO 4 1=1, MA 
DO 2 J-'l , MA 
W0RK1 ( I , J )=0 
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W0RK2C l,J)=0 
PHI (I , J )=0 
2 CONTINUE 
WORK1U,i)=1 
W0RK2( J, l)=TAU 
PHI ( l, l)=1 
4 CONTINUE 
00 6 1=1, HA 
DO 6 J=1,MB 
6 THETA(|,JJ=Q 
D I FMAX=1. E8 
NO 0=50 

IF(NTERMS,GT.O) NDQ=NTERMS 
f AC 1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMULtDt M Y, W0RK1 , A,TS, 1 , 1, TO, 10,NA,NA,MA,MA,MA,MA,1 ) 

FAC1=FAC1*H 

C0N1 = 1,/FAC1 

IFfNTERMS.EQ.l AND. I.GE.4J CALL SERROR( PHI ,Y?0RK1 , CONI ,NA, HA, 

1 TO , TO, HA, MA, D I FMAX ) 

CALL MXADD ( DUMMY, PHI , WQR K1 , Cl , CON 1 , 1 , 1 , NA, NA, 1 0, 1 0 , MA, MA, 1 ) 

I F(MODE. EQ.2) GO TO 500 
FAC2=FAC1*( 1+1) 

CQN2=TAU/FAC2 

CALL MXADQ( DUMMY, WQRK2, WORK! , Cl , C0H2, 1,1,10, IQ. 10, 10, MA, MA, 1 J 
500 CONTINUE 
JTER=I 

I F(NT£RMS.GT.O,OR, I , LT.4) GO TO 1000 
I F( D I FMAX.LE.TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

I F( MODE. EQ, 1 ) CALL MXMUL(THETA,W0RK2,B,C1,NA,NB,10 > 10,NA,NB, 

1 MA, HA, MA, MB, 2 ) . 

t F{ I TER . EQ , NDO . AND . NT ERMS . EQ . 0 ) I FLAG= 1 
RETURN 
END 


SUBROUTINE SERRQRf AMX, BMX, CGC, IA.JA, IB,JB, IDO, JDO,DI FMAX) 
DIMENSION AMX( IA.JA ),BMX( IB, JB) 

D(FMAX=1.E-30 
DO 100 1=1, IDO 
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DO 50 J=7,JD0 

tF(AMK{ l,J).EQ.0.0 ) GO TO 50 
CHANGE= ABS(8HX{ [ , J )*GCC/AMX( I , J ) ) 

I F( CHANGE. GT.D I FHAXJ D ! Ff<!AX=CHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 

■r> 
v- • 

c 

SUBROUT I NE FUTUREfX, V, PHI , THETA, LP, HP, LT.HT, I P, JP, 3T, JT.MQDE) 
DIMENSION X{MP),V(MT), Pill (LP*MP),THETA( LT,HT),TEMP{20) 

DO 20 1=1, IP 
SUM=0 

DO 10 J=1,JP 
70 SUM=SUM+PHi{ l,J)*X(J) 

20 TEMP( I )=SUH 
DO 30 1=1, IP 
30 X( I )=TEMP( I ) 

I F(MODE. EQ.2J RETURN 
DO 60 1=1, IT 
SUM=0 

DO 50 J=1 , JT 

50 SUM=SUM+THETA( I , J )*V( J ) 

X{ I )=X{ I )+SUM 
60 CONTINUE 
RETURN 
END , 

C 

G 

SUBROUTINE MXADD(RMX, AMX, BMX,ACC,BCG, IR,JR, IA,JA, IB, JB, IDO,JDO, 
. 1 MODE) 

DIMENSION AMX{ IA,JA),BMX( IB,JB),RMX( IR,JR) 

IFf IA.LT, IDO.OR.JA.LT.JDO) GO TO 999 
J Ft IB. LT. T DO.OR.JB.lt, JDO) GO TO 999 
GO" T 0 (10, 100), MODE 
10 CONTINUE 

DO 50 1=1, IDO 
DO 50 J=1 , JDO 

50 AMX( J , J )=AMX( I , J )^ACC+BMX( 1 , J }*BCC 
GO TO 300 
100 CONTINUE 

I F( IR.LT, IDO.OR.JR.LT, JDO) GO TO 999 
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DD 200 1=1, IDO 
DO 200 J=1,JDO 

200 RHXf I , J )=AMX( I , J )*ACC+BMX( I , J )*BCG 
300 RETURN 
999 CONTINUE 
RETURN 
END 


SUBROUTINE HXMUL(RMX,AHX,BMX,CCC, IR,JR, IA,JA, IB,JB, IDA, JDA, IDB, JDB 
1 , MODE) 

DIMENSION AMX( IA,JA),BMX( IB,JB),RMX{ IR,JRJ,T£MP{20) 

IE( IDA*IDB*JDA«JDB.GT. IA*IB*JA*JB) GO TO 999 
I F{ JDB.GT.JDA) GO TO 999 
GO TO (10 ,210), MODE 
10 DO 100 1=1, IDA 
DO 20 L=1 , JDA 
20 TEMP{ L)=AMX( l,L} 

DO 80 J=1,JDB 
SUM=0 

DO 1(0 K“1 , JDA 

40 SUM=SUM+TEMP(K)*BMX(K»J) 

AHX( I , J )=SUM*CCC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 
210 CONTINUE 

!F{ IR.LT. IDA.OR.JR.LT.JDB) GO TO 999 
DO HOD 1=1, IDA 
DO 380 J=1, JOB 
SUM=0 

DO 340 K=1 , JDA 

340 SUM=SUM+AMX( I , K)*BMX( K, J ) 

RMX( I , J )=SUM*CCC 
380 CONTINUE 
NOO CONTINUE 
600 RETURN 
999 CONTINUE 
RETURN 
END 


2 o 

**1 TO 

-u 55 


§2 

§?Z 

to 


to 

C »** 
^ US 

>&• 


Appendix B 


SOFTWARE FOR 

STABILITY ANALYSI S OF A HOCK REGULATOR 


325 


c 

G * THIS PROGRAM COMPUTES THE CLOSED-LOOP POLES OF A BUCK CONVERTER * 

C * WITH A SINGLE-STAGE INPUT FILTER, AS A FUNCTION OF THE INPUT * 

C * FILTER PARAMETERS, * 

C 

DIMENSION A(7)yBl6),CCl(7),CC2{5) 

REAL: NXEG, IXR,M,K7,K2,K3,K4,K5,R6,K7,L,L1,K8 

1 nteger NDEG, I ER 

COMPLEX Z(6),Zlf5),Z2(6),Z3(4) 

DATA VI, L,RQ,RC,C,RL, VO/25. 0,230. 0E-6,0.2,0.0G7,30G.0EH5, 10.0,20.0 
V 

DATA RN-Rll >R12,R73,R74,R4,C2,GP7 ,M/0.65,33.3E3,76.7E3,2.0E5, 
147.0E3,40.7E3,1.0E-70,5600.0E-12,0.88E-3/ 

C ' 

DO 2000 1=1,7000 
READ (5,*) RL1,L1,C1 
I F( RL1 . LT. 0. ) GO TO 2007 
0=0.8388579 
C 

Y1=RL*L7*C1*C*RC 

Y2=RL*L7*C7+C*RCtt{C1*RL1*RL-{D«*2)*Ll ) 

Y3=C7*RL1*RL-( p5W2)-*L1*C*RC»{ RL~( D**2)«RL1 ] 

YI|=RL-(D**2)*RL1 

Y5=RL*L1«C7*C*RL 

Y6=RL*L7*C7+C«RL*(C7*RL'l*RL~£Dtf i *2)*L7) 
Y7=C1*RL7«RL-fD**2)*L7‘fC*RL*{RL-(D**2)*RL1 J 
Y8=L*C*RL*L7*C7 

Y9= (' U*C*RL#C1*RL1 }+( C*RL*L1*C7*{ RO+RC+( L/(C*RL) J) ) 

Y1 G=L*C»RL+C*RL*C7mi«(R0+RC+( L/(C*RL> 5 }+L7*C7*RL*£ D»*2)*C#RL*L1 
Y1 1= C*RL* ( RO+RC+ ( L/ { C*RL J ) }+C1 *RL7*RL+ ( D ) *( L7*FD*RL' K RL7 ] 

Y72=RL+(D**2)*RL7 
RX={R1T*R72)/(R77+R72) ■ 

G=RX/R11 

K7=G/(R74+RX) 

K2=C2*(7.0*K1»R73) 

FM= ( 2 . 0*R4*C P7 ) / ( RN ft M ) 

K3=V0*RL*R4*FM 
K4=RH*L«C2*R13 
K5-V0*FM 
K6=CP1«C2»R13 ■ 

K7=D«RL*Rlj 
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K8=C*(R0*RC+(L/t C*RL) ) ) 

C 

C * THIS CALCULATION IS FOR THE CLOSED LOOP POLES WITH # 

C *' INPUT FILTER. * 

c 

A( 1 )=K7*K6*Y8 

At 2 )=K7*Cp1*Y8*K7*K6*Y9+K5#K4*Y5 

AC3)=K7*CP1*Y9+K7*K6*Y10+K3*K2*Y1+K5*RN*L*Y5+K5«K4*Y6 

A{ 4 J =K7«C PI *Y1 0+K7*K6*Yl 1 +K3*K2*Y2+K3*K1*Y1 +K5*RN*L*Y6+K5*K4»Y7 

At5)=K7*CP1*Y11+K7*K6*Y12+K3*K1*Y2+K3*K2«Y3+K5*RN*L*Y7*K5*K4*Y4 

At 6 )= CPI *K7*Y1 2+K3* K1 *Y3+K3*K2*Y4+RN*L*K5*Y4 

A(7)=K3*K1*Y4 

NDEG=6 

CALL ZPOLR(A*NDEG,Z, IER) 

C 

V/RITE(6, 101 } RL1,L1,C1 

101 FORMAT(/,10X, * RL1=', FI 5. 10, r L1=' , F15. 10, 1 C1=*,F15.10) 

WRITE (6, 103) 

103 FORMAT { /, 20X, 1 THE CLOSED LOOP POLES WITH [.FILTER ARE-') 

WRITEt 6,102}(Zt IX), IX-1,6) 

102 FORMAT* /,20X,2E15. 8) 

C 

B{ 1 )=K5*K4*Y5 

B { 2 )=K5*RN»L*Y5*K5 tt K4*Y6*K3*l<2*Y1 

B( 3 H<5*RN*L«Y6*K5*K4m+K3*K2*Y2+K3*Kmi 

B 1 4 } -K5*RN*L*Y7+K5*K4#Y4+K3*K1 *Y2+K3*K2*Y3 

B { 5 } =RN*L*K5*Y4+K3*K1 *Y3+K3*K2»Y4 

B(6)=K3*K1«Y4 

N0EG=5 

CALL ZP0LR(B,MDEG,Z1, IER) 

WRITEt 6, 104) 

104 FORMAT {/,20X, ' THE OPEN LOOP ZEROES ARE *) 

WRlTEt6,102)tZ1( IX), IX=1,5) 

C 

CG1(1)=K7**6*Y8 

CCl (2}=K7*CP1*Y8+K7*K6*Y9 

CC1(3)=K7*CP1*Y9+K7*K6#Y10 

CCl (4HK7*CP1*Y10+K7*K6*Y11 

CCl ( 5 )=K7*CP1»Y1 1+K7*K6*Y12 

CC1{G)=CP1*K7*Y12 

CC1f7)=0.G 
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CALL ZP0LR(CC1,6,Z2, IER) 

WRITE(6, 105) 

105 FORMAT [/, 20X, ' THE OPEN LOOP POLES ARE 1 ) 

WRITE(6,102)(Z2{ IX), IX-1,6) 

C 

c I t I I » Hit H It it » ft«- 

C ■* THIS CALCULATION IS FOR THE CLOSED LOOP POLES WITHOUT INPUT * 

C * FILTER. * 

C ■*»#»«*«*«»#*****»■»*«»»#»«»**»»»*«#*«*#»«***»#«*****»»«#»*«****#»*» 

c 


CC2£l)=K7*K6*L*C 

CC2(2)=K7*K6*K8+K7*CP1#L*C+K5 tt K4*C*RL 

CC2( 3 )=K7»CPl*K8+K7^S+K3*K2*G*RC+K5#K4+K5*RN*L*CaRL 

CC2(4)=K7*GP1+K3*K2+K3*Kl*CttRC+RN*L*l<5 

CC2(5)=K3*K1 

CALL ZP0LR(CC2,4,Z3, IER) 

WRITE{6,7B5) 

785 FORMAT ( / , 20X, * THE CLOSED LOOP POLES W/0 1. FILTER ARE-* 
WRITE(6,102)(Z3{ IX), IX=1,4) 

2000 CONTINUE 

2001 CONTINUE 
STOP 
END 


) 
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STABILITY ANALYSIS OF A BUCK CONVERTER WITH INPUT FILTER * 

THIS PROGRAM COMPUTES THE EIGENVALUES OF A BUCK REGULATOR 
SYSTEM WITH INPUT FILTER. 

DIMENSION FH6,6),F2(6,6},G1(6,4),G2{6,4LU(4bPHiF{6,6J, 
1PHIN(6,6),DF{6,4),DN{6,4),X{6),ZT(6),AA(6,6),BB(6),TEMP1{6,4), 
2PH1FT(6,6),DFT£6,4),XA(6),ZN1£6),XB(6) 

DIMENSION PSI(6,6),XN1(6}/FUN1(6),FUN2C6),C0N{4) J RVAL(6,9) 
DIMENSION TEMP2(6,6),DX(6) 

DIMENSION WKAREA£ 70 ) , WK( 70 ) 

REAL L1,L,K1,K2 
COMPLEX W(6),Z(6,6),ZN 

DATA RL1 J L1,C1,L J RO,RN,RC,C/0.2000,50.0E-6 J 220.E-6,230.E-'6,0.2EO 
l,0.65EO,0.067EO,300.E-6/ 

DATA RL, R1 1 , R12, R1 3, R14,C2, R4, CPI/10. OEO, 33 . 3E3, 1 6.7E3,2. £5, 

147. E3,0.01E-8,40.7E3, 5600. E-12/ 

DATA ER, ET, VI, EO/6 . 7E0, 7 . EO , 25 . EO, 20 . EO/ 

DATA VI TON, EQ, £D/0.88£-3,Q.2E0,0.7E0/ 

DATA N I T, EPS/90, 1 0 . E-6/ 

DATA EPS1, ERROR/2. E-6, 0.5E-6/ 

DATA CON/0 . 1 E-l ,0 . 05E-1 , 0 . 025 E-1 , G * 0125E-1 / 

DATA TOLL/0. 01/ 

Kl={ R12/( R14*( R1 1+R12) ) )+{ 1 ./R13)-(RN/R4) 

K1=K1-(RN*R0*'RL*RC) )/(R4*RL*RC) 

K2= £ RN/R4 ) - ( 1 . /R1 3 ) - ( R1 2/ ( R1 4* £ R11+R 1 2 ) ) ) 

TON=V I TON/VI 
, TOFFfi=TON*( ( V I -ED J/EO ) 

DO 1 1=1,6 
00 2 J=1 , 4 
G1( l,J)=0. 

G2{ l,J)=0. 

DF( l,J)=0. 

DN{ l,J)=0. 

CONTINUE 
DO 3 K=l, 6 
FI ( 1 , KJ=0. 

F2( l,K)=0. 

AA( I , !<)=0. 

PHI F( 1 ,K)=0. 
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PH I W ( I , K>=0. 

3 CONTINUE 

I CONTINUE 

DO 5 1=1,6 
ZT( J )=0.: 

BB{ I )=0 . 

FUNK I )=0. 

FUN2( 1J=0. 

KA£ 1 )=Q. 

KB( I )=0. 

XN1( I )=0. 

ZN1 ( I )=0» 

5 CONTINUE 

)=£-1.*BLl)/U 
n(1,2)=-T.00/L1 
FI (2,1 )=1. 00/C1 
F1(2,3)=-1.00/C1 
F1(3,2)=1 ,00/L 

FI ( 3, 3 )={ -1 , 00*RO) /L+( -1 . 00*RC*RL)/{ L»( RL+RC) ) 
FI { 3,ii)=“1 .OQ w RL/(L w (RL+RG) ) 
FK4,3)=RL/(C*(RC+RL)) 
F1(4,4)=“1.00/(C*(RC+RL)) 

F.T(*», 3 )=( RL*RC)/( C2*R1 3*£ RC+RL) ) 
f 1 (5,4 )=RL/£ C2*R13*( RC+RL) ) 
F1(5,5)=-1.0Q/(C2*R13) 

FI (6,2)={ -1 .00*RN)/(CP1«R4) 

FI { 6, 3 )=( -1 . 00*RL*RC*K1 )/( CP1*( RC+RL ) } 

F1{ 6, 4 )=( RL*K2 )/( CP1*( RL+RC) ) 
F1(6,5)=1.00/(CP1^R13) 

G1'( 1,1 J=1 .,00/LI 
G1 ( 3, 3 }=-l .00/L 
G1(6,2)=1.00/(CP1«R14) 

G1(6,3)=RiV(CP1*R4) 
ucu-vi 
U{2)=ER 
U( 3 )=EQ 
U(4)=E0 
DO 4 1=1,6 
DO 4 J=l,6 
F2( I ,J)=F1 ( I , J ) 

F2(2,3)=0.00 
F2(3,2)=0.00 
F2{6,2)=0.00 


4 
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G2(3,4)=»1.00/L 
G2(6,4}=RN/(GP1*R4) 

C 

T=TON 

CALL STRANf T, PH ) N, DN, Fl , G1 , 6, 4,6, 4, 1 , 0, ERROR, D 1 FHAX, ITER, I FLAG) 
T=TOFFB 

GALL STRAN ( T, PH I F, DF, F2, GS , 6 , 4, 6, 4, 1 , 0, ERROR, D I FMAX, t TER, 1 FLAG ) 

CALCULATION OF THE APPROXIMATE STEADY STATE * 

DO 5000 1=1,6 
DO 5000 J=1 , 6 
AA{ I , J )=0. 

DO 5001 K=1,6 : 

5001 AA( ) , J )=AA( I , J }+PHl’N( 1 , K)*PH i F(K, J) 

AA( I , J )=~1 .*AA( ! , J ) 

5000 CONTINUE 

DO 5002 1=1,6 

5002 AA( I, I )=1 . +AA( I , I ) 

DO 5003 1=1,6 

DO 5003 J=1 ,4 
TEMPI ( !,J)=0, 

DO 5004 K=1 , 6 

5004 TEMPIt l,J}=TEMP1( I,J)+PHIN( l,K)*DF(K,J) 

5003 TEMPI ( l,J)=TEMP1( l,J)*DN( t,J) 

DO 5005 1=1,6 

BB(1)=0. 

DO 5006 J=1,4 

5006 BB( I )=BB( 1 )+TEMP1 ( l,J)4UfJ) 

5005 CONTINUE 
MM-1 
NN=5 
IAA-6 
IDGT=0 

CALL LEQT2F(AA,MM,NN, I AA, BB, 1DGT,WKAREA, IER) 

DO 5007 1=1,5 

5007 X( I )=BB( I } 

X(6)=ET-PHlF(6,1)*X{1)-PHIF(6,2)*XE2)~PHiF(6,3)*X(3)~PHlF(6,4)*X{ 

l4)-PHIF{6,5)*X(5}-DF{6,1)*U{1)-DF(6,2)*U(2)-DF(6,3)m(3)-DF{6,4)* 
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211(4) 

WR1TE(6,5008) 

5008 FORMAT (1 HI , 37X, ' APPROXIMATE STEADY STATE VALUES,X= *) 
WR|TE(6,5009)(X( f ), J=1,6) 

5009 FORMAT (1H ,6F15.10) 

WRITE(6,5010) TOFFB 

5010 FORMAT(1HO,37X, * APPROXIMATE VALUE OF OFF TIME— *,Fl5.10) 

C 

C CALCULATION OF THE EXACT STEADY STATE VALUES * 
c 

IT1-0 

900 DO 26 1-1,5 
ZT( l )=0 . 00 
DO 27 K=T, 6 

27 ZT(l)=Zt(l)+PHIF(l,K)*X{K) 

DO 28 J=l,4 

28 ZT( I )=ZT( I J+DF( I ,J 
26 CONTINUE 

ITT-IT1-M ; ■: : 

SMAT-X(6)-PH1 N( 6, 1 }*ZT( 1 )-PH I N( 6,2)*ZT( 2)-PH I N{ 6, 3)*ZT{ 3 )-PH IN(6, 4 
iy*ZT{4)-PHIN(6,5)*ZT{5)-PHIN{6,6)*ET-DN(6,1)*un)~DN(6,2)*U(2)- 
2DN( 6, 3)*U(3)-DN(6,4)*U(4) 

IF(ABS($MAT).LE.EPS1) GO TO 70 
I F ( IT1.GT.60) GO TO 70 

dToff-o.oi*toffb 

T=TOFFBFDTOFF 

CALL STRAN(T, PHI FT, DFT, F2,G2,6,4,6,4, 1,0, ERROR, D I FMAX, ITER, I FLAG) 
DO 29 1=1,6 
DO 29 J-1,6 
AA( I, J )-0.00 
DO 56 K=1, 6 

56 AA( I,J)=AA( i,J)+PHIN( I ,K)*PHI FT( K, J ) 

AA( l,J)=-1.00*AA( I , J ) 

29 CONTINUE 

DO 31 1=1,6 

31 AA{ J, I )=1.00+AA( I, I ) 

DO 32 1=1,6 
DO 32 J=1 ,4 
TEMPI ( !,J)=0. 00 
DO 34 K=1,6 

34 ■ TEMPI: ( l ,J)=TEMP1( i ,J)+PHIN( l»K)*DFT(K,J) 
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32 TEMPT ( 1 , J )=TEMP1 ( 1,J)*0N( f,J) 
DO 35 1=1,6 
BB{ I J=0. 00 
DO 36 J=1,4 

36 BBf I )=BB[ I )+TEMP1{ J,J)*U{J) 

35 CONTINUE 
MM=1 


MN=5 
! AA=6 
IDGT=0 

CALL LEQT2F(AA,MM,NN, lAA^BB, 1DGT,WKAREA, )ER) 

DO 37 1=1,5 
XN1 ( I )=BB{ I, J 

XN1(6) = ET-PHtFT{6 5 1)*XNin)~PHiFT{6,2)*XN1{2)-PHIFT(6,3)«XfJl(3)- 
1PHJFT[6,4)*XNT(4)--PHIFT(6,5)*XN1(5)*-DFT(6,1)*Un)“DFT(6,£)*U(2)-- 
2DFT(6,3)' tt U(3)-DFT(6,4)' a U(4) 

DO 38 1=1 , 5 
ZN1 ( i )=0. 00 
DO 39 K=1,6 

ZN1( I )=ZN1{ I )+PHlFT{ l,K)*XN1(K) 

DO 40 J=1,4 

ZN1 1 i )=ZN1 ( I)+DFT( I ,J) W U( J J 

CONTINUE ’ 

SMATN=XN1 f 6)**PHIN(6, 1 )*ZN1( 1 J-PHIN(6,2)*ZN1 (2)-PHlN{6,3 )*ZN1{3}" 
1PHIN(6,4)*ZNi(4)-PHIN(6,55'«ZN1{5}--PHlM(6,6)'*ET-DN{6,1}*U(1)~DN|;6,2 
2PUL2)-DN(6,3)*U(3)-DN(6,4}*IJ(4) 

D SMAT=SMAT N - SMAT 

TQFFB=TOFFB+( (-1.00*SMAT)/{DSMAT/DTOFF) } 

T=TOFFB 

CALL STRANfT, PH I F, DF, F2, G2, 6, 4, 6, 4, 1 , 0, ERROR, D I FMAX, I TER, I FLAG) 

DO 42 1=1,6 
DO 42 J=1 , 6 
AA(1,J)=0.00 
DO 43 K=1 ,6 

AA{ !,J)=AA( 1,J)+PHtN{ i , K)*PHI F(K, J) 

AA( l,J)=-1.00*AA( !,J) 

CONTINUE 

DO 55 1=1,6 

AA( I , I )=AA( l> I )+1 . 

DO 44 1=1,6 
DO 44 J=1,4 
TEMPI ( I, J )=0, 00 
DO 46 K=1>6 
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46 TEMPI ( J,J)=TEMP1( l,J)+PHIN( l,K)*DF(K,J) 

44 TEMPI ( I , J )=TEMP1 { I , J)+DN£ 1 , J J 

DO 47 1=1,6 

BB{ 1 )=0 . DO * 

DO 48 J=1,4 • 

48 BB( I )=BB( l )+TEMP1 ( l,J )*U(J) 

47 CONTINUE 
MM=1 
NN=5 

1AA =6 
1 DGT=0 

CALL L£QT2F(AA,MM,NN, IAA, BB, I DGT, WKAREA, IER) 

DO 49 1=1,5 

49 Xf 1 }=BB( 

X{6)=ET-PIIIF{6,1)«Kn)-PHIF(6,2)*X[2)-PH1F{6,3)*X£3)-PHlF[6,4)* 

1X(4)-:PHIF(6,5)*X{5)-DF{6,1)*Un)*-DFf6,2)*U{2)-DF(6,3)->HJ(3)- 

2DF(6,4PU(4) 

GO TO: 900 
70 CONTINUE 

WRITE(6,50) 

50 FQRMAT( lllO, 37X, * EXACT STEADY STATE VALUES, X= ') 

WRITE£6,199)(X( I ), 1=1,6) 

199 F0RMAT( III ,6F15.10) 

WRITE? 6, 198) TOFFB 

198 FORMAT? 1 HO, 37X, * EXACT OFF TIME= *,F15.10) 

DVAL=TON/ ( TO N-i-TO FFB ) 

WRITE (6, 624 ) OVAL 

624 F0RMAT(1H0,37X, ’ VALUE OF 0=*,F15.10) 

WRITE? 6, 51) ITS 

51 FORMAT? 1 HO, 37X,' NO. OF ITERATIONS REQD.= *,14) 

C 

C ■* CALCULATION OF THE MATRIX PSI AND ITS EIGENVALUES * 

C ■ ■' 

DO 2000 1=1,6 
DO 2010 IRICH=1,4 
DO 61 12=1,6 
XA{ IZ)=X( IZ) 

61 XB(IZ)=X( IZ) 

DX( I )=CON( IRICH)*ABS(X( I)) 

IF(ABS£X£ l)).LE.TOLL) DX£I)=0.01 
XA£ I )=X( I )+DX( I ) 


hi 

UJ 


XB( I )=X( I )-DX( I ) 

TFFG1=T0FFB 

IT=0 

T8=T0FFB 

CALL STRANt T8, PHI F,DF, F2,G2,6,4,6,4, 1 ,0, ERROR, Dl FMAX, ITER, I FLAG) 

1101 Zm=PHfF(6,3)*XA(3)+PHIF(6,4)* 
1XA(4)^PHIFf6 > 5]l*XA(5)+1.0«XA{6)+DF(6 J 1) w U(1)+DF(6 > 2}«U(2)-F 
aDFtOjSj^UOJ+DFtG.UJ^UtljJ-ET 

I F ( ABS ( ZETA) . LE, EPS ) GO TO 1102 
t F( IT.GE.NIT) GO TO 1003 
DTOFF=CON( IRICH)*TFFC1 
T=TFFC1+DTOFF ' 

CALL STRAN{T, PHI FT, DFT, F2, 02,6,4, 6,4, 1 , 0, ERROR, D I FMAX, ITER, 1 FLAG) 
ZETAN=PHIFT(6,3)*XA(3)-f-PHIFT(6,4 

1 )*XA( 43+PH I FT( 6, 5}*XA( 5 )+1 . 0«XA( 6 )+DFT( 6, 1 )*U( 1 J+DFT( 6, 2)* 
2U(2)+DFT(6,3)*U(3)*DFT(6,4)*U(4}~ET 
DZ ET A=Z ETAN ~Z ETA , 

SLOPE^DZETA/DTOFF 

TFFC1=TFFC1«{ZETA/SL0PEJ 

T=TFFC1 

CALL STRANf T, PH I F, DF, F2, G2, 6, 4, 6, 4, 1 , 0, ERROR, D I FMAX, ITER, IFLAG) 
IT= JT+1 
GO TO 1101 

1102 DO 13 J=1, 6 
TEMP=0. 00 - ' 

DO 14 K-1 ,6 

DO 81 MENT=1,5 

8i phif(m£nt,6)=o.o 

PHI F(6,6)=1 .0 

14 TEMP=TEMP+PHI F{J,K)*XA(K) 

DO 15 K=1 , 4 

15 T£MP=TEMP+DF(J,K)*U(K) 

FUN1 ( J )=TEMP 

13 CONTINUE 

TFFC2=T0FFB 

IT-0 

T=T0FF8 

CALL STRAN(T, PH I F,DF, F2, G2,6,4,6, 4, 1 ,0, ERROR, D I FMAX, ITER, IFLAG) 

1103 ZETA=PH]F(6,3)*XB(3)+PHIF(6,4)tt 
1XB(4)+PIUF{6,5)*}<B(5}+1.0*XB(6)+DF(6,1)*U(1)-5DF(6,2)*U(2}+ 
2DF(.6,3)*U[3)+DF(6,4)*U{4)-ET 

IF(ABS(ZETA).LE.EPS) GO TO 1104 
IF( IT,GE.NIT) GO TO 1003 
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DTOFF=CON{ 1 R I CH )*TFFC2 
T=TFFC2+DT0FF 

CALL STRANt T, PHI FT, DFT, F2,G2,6,4, 6,4, 1,0, ERROR, D1 FMAX, ITER, t FLAG) 
2 £TAN= PH I FT ( 6 , 3 ) *XB ( 3 ) + PH I F T ( 6 , 4 

1)' if XB(4)+PHIFT(6,51*XB(5)+1.0*XB(6)+DFT(6,1>*U(1)*DFT(6,2)* 
2U(2)WT(6,3)->rut3)+DFT(6,4)*u{4)-ET 
DZETA=2ETAN-ZETA 
SLOPE=DZETA/DTOFF 
T F FC2=T FFC2- { ZET A/S LOPE ) 

T=TFFC2 

CALL STRAN{T, PHI F,DF, F2,G2,6,4,6,4, 1,0,ERROR,DI FMAX, ITER, I FLAG) 
IT= I T+1 
GO TO 1103 
1104 DO 62 J=1,6 
T£MP=Q.OO 
DO 63 K=1 ,6 
DO 82 MENT=1,5 
82 PHI F(M£NT,6)=0.0 
PHI F(6,<5 )=1 .0 

63 TEHP=TEMP+PH I F{ J, K )*XB( K) 

DO 64 K=1,4 

64 TEMP=TEMP+DF(J,K)*HJ(K) 

FUN2( J )=TEMP 

62 CONTINUE 

DO 2011 J=1 ,6 

2011 RVAL{ J, IRICH)={FUN1(J)-FUN2(JJ )/(2.0*DX( I ) ) 

2070 CONTINUE 

DO 65 J-1,6 

RVAL(J,5MRVAL£J,2)-0.25*RVALfJ. 1)1/0.75 
RVAL£«j> 6)=£ RVA{.£ J,3)-0.25*RVAE( J,2) }/G«75 

65 RVAL( Jj7)=( RVAL(J f 4)~0.25^RVAL{ J>3 ) )/0. ?5 
DO 66 J=1 ,6 

RVAL(J,8)=(RVAL(J > 6)-0.0625*RVAL(J,i>))/0,9375 

66 RVALf J, 9)"£ RVAL{ J, 7 )-0. 0625*RVAL( J , 6 ) }/0. 9375 
DO 67 J— 7 , 6 

67 TEMP2£J, I )=£RVAL£J,9)-0. Ol5625 4f RVAL£J,8) )/0. 984375 
2000 CONTINUE 

’DO 17 j=1 ,6 
DO 17 J=1,6 
TEMP=0.00 
DO 79 K=1 ,6 

19 TEMP^TEMP-JPNINC !,K)*TEMP2{K,J) 

PSI £ l ,:J )=TEMP 
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17 CONTINUE 

WRIT E{6,92) 

p2 FORMAT! 1 HO, 37X, ' MATRIX PS I— ') 

DO 30 1=1,6 ‘ 

30 WR ITE( 6,93) { PS I ( I, J)>J=1,6) 

93 F0RMAT(6F15. 10) 

GO TO 1004 

1003 WRITE! 6, 100) I 

100 FORMAT! //> ' CONVERGENCE NOT OBTAINED FOR X{J),I= *,14} 

1004 CONTINUE 

CALL E l GRF( PS 1 , 6, 6, 2, W,Z, 6, WK, TER) 

WRITEf 6, 94) 

94 FORMAT! 1 HQ, 37X, ’ THE EIGENVALUES ARE— f ) 

WRITE! 6, 103 J (W( -6) 

103 FORMAT! Ill .6F15.10) . 

WRITE! 6, 104) WK[1), IER 

104 FORMAT! 1 HO, 37X, * CONVERGENCE TOLERANCE=% FI 5,10, r IER=’,I4) 
905 CONTINUE 

STOP 
END 


SUBROUT I NE STRANfTAU, PH I , THETA, A, B, NA, NB, MA, MB, MODE, NTERMS, 
1 TOL, 0 1 FMAK, i TER, I FUG ) 

REAL PHI I NA, NA), THETA{ NA,NB),A( NA,NA) ,B(NA, NB), 

1 WORK1 ! 1 0, 1 0 ) , WORH2 C 1 0 , 1 0 ) , DUMMY! 1 , 1 ) , FACT , CON 1 , FAC2, C0N2 
C1=1 
TS=TAU 
DO 4 l-1,MA 
DO 2 J=T,MA 
V/ORKi [ I , J)=0 
W0RR2( r,j)=o 
PHI ( I ,J)=0 

2 CONTINUE . 

W0RKin,|J=1 
W0RK2( I, I }=TAU 
PIIK I, t)=1 
4 CONTINUE 
DO 6 1=7, MA 
DO 6 J=1,MB 
6 THETA{ I, J)=0 
0 1 FMAX=1 . E8 
NDO=50 
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r;-r;a=r» 


I F{NTERMS. GT.O) NDO=NTERMS 
FAC1=1 
I FLAG=0 

DO 1000 1=1, NDO 

CALL MXMU L ( DUMMY, WOR K1 , A, IS, 1 , 1 , 1 0, 1 0, NA, NA, MA, HA, MA, MA, 1 ) 

FAC1=FAC1*l 

C0N1=1./FAC1 

J Ff NTERMS. EQ.Q. AND. I .GE.4) CALL SERRORf PH I ,W0RK1 , CONI , NA, NA, 

1 10,10,MA, MA, Dl FHAX) 

CALL MXAD.D£ Dl/MHY, PHI ,WORK1,C1,CON1 ,1,1, NA, NA,10, 10,MA, MA, 1 } 
IF(M0DE.EQ.2) GO TO 500 
FAC2=FAC1*t 1*1} 

C0N2=TAD/FAC2 ' 

CALL MXADD{ DUMMY, W0RK2, WGRK1 , Cl , C0N2, 1, 1 ,10,10,10, 10,MA,HA, 1 ) 
500 CONTINUE 

ITER= I . ' . 

1 F( NTERMS. GT.O; OR. I. LT.4) GO TO 1000 
1 F ( D I FMAX. LE. TOL) GO TO 1100 
1000 CONTINUE 
1100 CONTINUE 

I F( MODE. EQ.1 ) CALL MXMUL( THETA, V/0RK2,B, Cl , NA, NB, 1 0, 10, NA, NB, 

1 MA, MA, MA, MB, 2 } 

I F( ITER. EQ.NDO. AND. NTERMS. EQ.Q) 1FLAG=1 

RETURN 

END 


SUBROUTINE SERROR{AMK, BMX.CCC, IA,JA, IB,JB, IDO, JDO,DI FMAX) 

D I HENS 1 ON ANX( I A, JA ) , BMX( t B, J B } 

Dl FMAX=1 . E-30 
DO 100 1=1, IDO 
DO 50 J=1,JD0 

1 F( AMX( 1 , J ) ,EQ. 0 , 0 } GO TO 50 
GijANGE= ABS(BMX( 1,J)*CCC/AMXU,J}) 

I F( CHANGE, GT.D I FMAX) D1FMAX=CHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 
C 

c 

SUBROUTINE MXADD(RMX,AMX,BHX,ACC,BGC, IR,JR, 1A,JA, IB,J8, 1DO,JDO, 
1 MODE) 
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DIMENSION AMX( fA,JA),BMX{ !B,JB),RMX( IR,JP.) 

I F{ IA. LT. IDO.OR.JA.LT, JDO) GO TO 999 
IFt IB.LT. JDO.OR.JB.LT.JDQJ GO TO 99 9 
GO TO (10,1 00); MODE 
10 CONTINUE 

DO 50 (=1, iOO 
DO 50 J=1,JD0 

50 AMX( I r J )=AMX( I , J )*AGC+8MX( J , J )*BCC 
GO TO 300 
100 CONTINUE 

IF( JR. LT. IDO.OR*JR, LT.JDO) GO TO 999 
DO 200 1=1, IDO 
DO 200 J=1 , JDO 

200 RMX( 1,J)=AMX£ J,J)*ACC+BMX(1,J)*BCC * 

300 RETURN 
999 CONTINUE 
RETURN 
END 
C 
C 

SUBROUTINE MXMUL(RMX,AMX,BHX,CCC, JR, JR, IA,JA, IB, JB, IDA, JDA,IDB,JDB 
1 , MODE ) 

DIMENSION AMX(IA,JA),BMX( 1B,JB),RMX( m,JR),TEMP(20) 

I Ff I DA^IDB«JDA*JDB.GT. 1A"1B«JA*JB) GO TO 999 
I F( JDB.GT. JDA) GO TO 999 • 

GO TO (10,210), MODE 
10 DO 100 1=1, IDA 
DO 20 L=1, JDA 
20 TEMPI L)=AMX(I,L) 

DO 80 J=1, JOB 
SUM=0 

DO J|0 K=1 ,JDA 

HO SUM=SUH*TEMF(K}*BMX{K,J) ' ‘ 

AMX( I , J)=SUM W CCC 
80 CONTINUE 
100 CONTINUE 
GO TO 600 

210 CONTINUE . 

I Ft IR. LT. IDA.OR. JR.LT. JDB) GO TO 999 
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340 SUM=SUM+AMX{ I , K } *BMX ( K, J ) 
RHXi !,J)=SUM*CCC 
380 CONTINUE 
400 CONTINUE 
600 RETURN 
999 CONTINUE 
RETURN 
END 
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C 

C * STABILITY ANALYSIS OF A BUCK CONVERTER « 

D 

C THIS PROGRAM COMPUTES THE EIGENVALUES OF A BUGK 
C REGULATOR SYSTEM WITH INPUT FILTER AND 

C FEEDFORWARD. 

C 

DIMENSION F1{7,7LF2(7,7),G1(7,4),G2(7,4),U{4),PHJF(7,7), 

1PHIN( 7,7),DF(7,4},DN(7,4),X(7),ZT(7),AA(7,7),BB(7),TEMPl(7,4}, 

2 PH I FT[7,7),DFT(7,4),XN1(7),ZNT(7) 

DIMENSION PSI{7,7),XA(7) J FUNH7J J FUN2l7) > XB{7),C0NtH),RVAL(7,9) 
DIMENSION TEHP2(7,7),DX£7) 

DIMENSION WKAREA{90},WK{90) 

REAL LI , L, Kl , K2,K3 
COMPLEX W{7),Z(7,7),ZN 

DATA RL'!,L1,C1,L,RO,RN,RC > C/O.20E0 > ?.425E-3,220.E-6 > 23D.E-6,0.2EO # 
10.65EO,0.067EG,300.E-6/ 

DATA RL,R11,R12,R13,R14,C2,R4,CP1/10.0EQ,33.3E3,16.7E3,2.E5, 

147. E3, 0.01E-8,40.7E3, 5600. E-12/ 

DATA ER,ET,Vl,EO/6.7EO,7.EO,25.EO,20.EO/ 

DATA VITON,£Q,ED/0.88E-3,0.7EO,0.2EO/ 

DATA NIT, EPS/150,5. £-6/ 

DATA EPS1/0.5E-6/ 

DATA ERROR/O. 5E-6/ 

DATA C0N/0.2 E-1 , 0 . IE- 1 , 0 . 05E-1 , 0 . 025E-1 / 

DATA GF,RF1, RF2, TOLL/2. 7E-5,5.1E3, 220. OE0,O. 11/ 

C 

Kl=( R12/{ R14*( R1 1+R12) ) )+( l ./R13 )-( RN/R4 } 

K1=K1- [ RN*RO*{ RL+RC ) )/{ R4*RL*RC) 

l<2={ RN/R4 ) - ( 1 ./Rl 3 ) *• ( R1 2/ ( R 1 4*{ Rll^-RI 2 J ) J 

K3=RF2/{RF1+RF2) 

T0N=V I TON/VI 
TOFFB=TON*({VI~EO)/EO) 

DO 1 1=1,7 
DO 2 J=1,4 
G1( l,J)=0. 

G2( I,J}=0. 

DF( I,J)=0. 

DN(I,J)=0. 

2 CONTINUE 
DO 3 K=T, 7 
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FI ( I , K)=0. 

F2{ l,K)=0. 

AA( J,K)=0. 

PHI F( 1 *K]“0. 

PHINC 1*10=0.' 

CONTINUE 
CONTINUE 
DO 5 1=1,7 
ZT{ I )=0. 

BB( I )=0, 

FUN1 ( I }=0. 

FUN2[ I )=0. 

XA[ I )=0. 

XB(l}=0. 

XN7{ i }=0. 

ZN1(I)=0. 

continue 

FT(1,1)*(-1i«RL1.J/U 

(f— -j J ^ n i rtn il 4 

n 

FI 
FI 
FI 
FI 

F1£3,4)=-1./C1 

Fl(li,21=1.0D/L 

F1(4,4}=M.a0*RO)/LH“1.00KRC*RJ-5/(l*(RL+RC)) 

FI (4, 5 ) =-1 . 00*RIV { L* ( RL+RC ) ) 
F1(5,4)=RL/(C#(RG+RL) } 

FI £5,5 )=-1 , 00/ ( C*( RC+RLH 

FI { 6,4 )=£ ■*!, *RC)/{ C2»R13*( RC+RL) } 

F1(S,5)=Ru/<C2*R13*(RC+RL)) 

FT{6,6)=-1.00/£C2«R13) 

F1{7,U“-1.*K3/C1 
FI £7,2}=(-1,00»RN)/£CP1*R4) 

FI ( 7 , 3 )=K3*( C1+C F)/£ Cl *CF*( RF1+RF2) ) 

FI £ 7, 4 }={ K3/C1 )-( ( RL*RC*K1 )/( CP1*{ RL+RC) ] ) 

FI £ 7, 5 ) = ( RL*K2 ) / { CP 1 *.( RL+RC) ) 

FI { 7, 6 )=1 . 00/CCP1*R13 ) 

61 (1,1 1=1. 00/L1 
G1 £ 4, 3 )=-1 , 00/L 
61 ( 7, 2 )=1 . 00/ £ CPI *R1 4 ) 

G1£7,3)=RN/(CP1*R4) 


\ I t c. j I » uu/ u I 

(2,1 )=1 ,00/C1 

(2,3 J=-1 »/(C1*(RF1+RF2JJ 

£2,4)=-1 ,Q0/C1 

(3,1 J=l./C1 

(3,3)=-1-*(C1-l-CF)/(C1*CF*£RF1+RF2)} 
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U(1)=VI 
U(2)=ER 
U( 3 )=EQ 
U(4)=ED 
DO 4 1-1,7 
DO 4 J=l,7 

4 F2f I ,J )=F1 ( I , J ) 

F2{2,4)=0. 

F2t 3,4)-0. 

F2{ 4,2 )—Q. 

F2f7,2J=0. 

F2(7,4 )={-!. *RL»RC»K1)/£CP1*CRL+RC) ) 

G2{1„1)=G1{1,1) 

G2(7,2)=G1(7,2) 

G2(4,4)=-1.00/L 

G2(7,4)=RN/(CP1*R4) 

C 

T=TON 

CALL STRAN ( T, PH I N, DN, FI , G1 , 7 , 4,7,4, 1 , 0, ERROR, D I FMAX, I TER, 1 FLAG ) 
T=TOFFB 

CALL STRAN(T, PHJF,DF, F2,G2, ?,4,7,4, 1,0, ERROR, D1 FMAX, ITER, I FLAG) 

*#**«■»« 

» CALCULATION OF THE APPROXIMATE STEADY STATE * 


DO 5000 1=1,7 
DO 5000 0=1,7 
AA£ I,J)=0. 

DO 5001 K=1,7 

5001 AA{ 1 , J}=AA( l,J)+PHIN{ I , K)*PHI F( K, J ) 

AA( l,J)=-1.*AA( 1,J) 

5000 CONTINUE 

DO 5002 1=1,7 

5002 AA( J , I )=1 ,+AA( 1,1) 

DO 5003 1=1,7 

DO 5003 J=1,4 
TEMPI { I ,J )=0. 

DO 5004 K=l,7 

5004 TEMPI (1,J)=TEMP1{ l,J)*PHIN( I , K)*DF{ K, J ) 

5003 TEMPI ( l,J)=TEMP1( J,J)+DN( l,J) 

DO 5005 1=1,7 

BB( I )= 0 . 
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00 5006 J=1 ,4 

5006 BB( l)=BB( I J+TEMP1 £ l,J)*U(J) 

5005 CONTINUE 

MM=1 

NN=6 

1 AA=7 

I DGT=0 

CALL LEQT2F ( AA, MM, NN, IAA,BB, I DGT, V/KAREA, 1 ER) 

DO 5007 1=1,6 

5007 >i( i )=BB( I ) 

X{7)=ET-PHIF(7,1 )^Cl)“PHIF(7,2)*xra)"PniF(7,3)*X(3)-PHlFf7,4)«X( 
l4)-PHIF(7,5)*X(5)-DF(7,1)*un)-.DF(7,2}*U{2)-DF{7,3)*U(3)-DF{7,4}* 
2U{4)-PHIF(7,6)*X(6) 

WR I T E( 6,5008) 

5008 FO>. Ar(1HT,37X, 1 APPROXIMATE STEADY STATE VALUES,X= * ) 
WRITE(6,5009)(X{ J ), 1=1,7) 

5009 FORMAT! 1H ,7F15.10) 

WR1TE(6,5010)T0FFB 

5010 FORMAT! 1 HO, 37X, 1 APPROXIMATE VALUE OF OFF TIME=* ,F15.1Q) ‘ 

* CALCULATION OF THE EXACT STEADY STATE VALUES * 

IT1=0 
900 DO 26 1=1,6 
ZT( I )=0, 00 
DO 27 K— 1 , 7 

27 ZT! I )=ZT( I J+PH1 F( l,K)*X(K) 

DO 28 J=1 ,4 

28 ZT! I)=ZT{ ! )+DF( l,J)*U(J) 

26 CONTINUE 

IT1=IT1+1 

SMAT=X! 7 ) -PH I N (7, 1 )*ZT{ 1 )-PH J N( 7, 2)*ZT( 2 ) -PH 1 N(7, 3 )*ZT( 3 )- PH I N( 7,4 
1)*ZT{4)-PH!N(7,5)^ZT(5)-PmN{7,7)*ET-DN(7,1)*U(1)-DN{7,2)*U{2)- 
2DN{7,3)*U!3)-DN(7,4)*U(4)-PHIN{7,6)«ZT!6) 

I F(ABS( SHAT) . LE. EPS1 ) GO TO 70 
IF( IT1.GT.60) GO TO 70 
DTOFF=0.01*TOFFB 
T=TOFFB+DTOFF 

CALL STRAN(T, PHI FT, DFT, F2,G2,7,4,7, 4, 1,0, ERROR, Df FMAX, ITER, J FLAG) 
DO 29 1=1,7 
DO 29 J=1 ,7 
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DO 56 K=l,7 

56 AA{ l,J)=AA{ I,J]+PHIN{ l,K)*PHIFT(K,J) 

AA( I, J)=~1 . OD*AA( I, J) 

29 CONTINUE 

DO 31 1=1,7 

31 AA(I,I)=1.00+AA(I,I) 

DO 32 1=1,7 

DO 32 J=1,4 
TEMPI ( J , J )=C,00 
DO 34 K=1,7 

34 TEMPI ( I , J J=TEMP1 { I, j)+PHIN( I ,K)*DFT(K,J} 

32 TEMPI { I , J }=TEHP1 ( 1 , J )+DN( l , J ) 

DO 35 1=1,7 

BB( I )=0.00 
DO 36 J=l,4 

36 BB( I )=BR[ l )+TEMP1( I 

35 CONTINUE 

MM=1 : ' 

NN=6 

IAA=7 

IDGT=0 

CALL LEQT2F{AA,MM,NN, IAA,BB, I DGT, WKAREA, 1 ER ) 

DO 37 1=1,6 

37 XN1( I }=BB£ I ) 

XN1 (7 J=ET**PH I FT(7, 1 )*XNl ( 1 )-*PHI FT( 7, 2)*XN1 (2)-PHIFT(7,3 )*XN1 (3)- 
lPHIFr(7,4)*XN1{4)-PHlFT|7,5)*XNl(5J-DFT(7,l)*U(l)-DFT(7,2J^U(2)- 
2DFT(7,3)«U[3}-DFT{7,4)*U£4)-PHIFT{7,6)-»XN1{6) 

DO 38 1=1,6 
2N1 ( I )=0, 00 
DO 39 K=1 ,7 

39 ZN1 { I )=ZNT{ I )*PH I FT( I , K)*XN1 ( K) 

DO 40 J=1,4 

40 ZNT( I)=2N1{ I )+DFT{ l,J)*U{J) ? 

38 CONTINUE 

SMATN=XNl(7)-PH!N{7,1)*ZNim-PHm{7,2)*ZNl(2)-PHIN£7,3)*ZN1(3)~ 

1 PH I N(7>4 }*ZN1 (4 )*-PH!N( 7,5) W ZN1 (5)-PHJ N(7,7) W ET-DN{7,1 )*U( 1 )-DN{7,2 
2)*U(2)-DN[7,3)*U(3)-DN(7,4)*U(4)-PH!N£7,6)*ZN1(6) 

DSMAT=SMATN-SMAT 

TOFFB=TOFFB+( ( -1 . OO w SMAT ) / ( DSMAT/DTOFF ) ) 

C 

T=TOFFB 

CALL STRANfT, PHI F, DF, F2,G2,7,4, 7,4, 1,0, ERROR, D I FMAX, ITER, t FLAG) 
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DO 112 1=1,7 
DO 42 J=1 , 7 
AA( l,J)=0.00 
DO 43 K=1,7 

43 AA{ I , J )=AA{ I , J )+PH I N{ 1 , K)*PH I F{ K, J ) 

AA(M)=-1.0Q*AA(l>J) 

42 CONTINUE ■ 

DO 55 1=1,7 
55 AA£ I , I )=AA{ 1i 
DO 44 1=1,7 
DO 44 J=1,4 
TEMPI { ) , J )=0,00 
DO 46 K=1,7 

46 TEMPI ( l , J )=TEMP1 ( I , J )+PHIN( J ,K)*DF( K, J ) 

44 TEMPI ( 1,J)=T£MP1( I, J)+DN{ l,J) 

DO 47 1=1,7 

BB( I)=0.00 • 

DO 48 J=1,4 

48 B3( I )=BB{ I )-KFEMP1{ 

47 CONTINUE 
MM=1 

NN=6 to 

I AA=7 *S* 

IDGT=0 cp> 

CALL LEQT2 F £ AA, MMji NN , 1AA, BB, IDGT, WKAREA, IER) 

DO 49 1=1,6 

49 X( I )=BB( 1 ) 

X(7)=ET-PHIF(7,1}*X(1 )-PHIF(7,2]*X(2)-PmF(7,3)*X{3)“PHtF(7,4)* 
1Xf4)“PlllF(7,5}«Xt5J-DF(7,1)«Un)-DF(7,2)*U(2)-DF(7,3)*U(3)*- 
2DF( 7, 4)*U(4)-PHI F(7,6)*X(6) 

GO TO 900 
70 CONTINUE 

WRITE(6,50) • 

50 FORMAT ( 1 HO, 37X, * EXACT STEADY STATE VALUES, X=’ ) 

VmiTE(6,l99)(X(M,l=1,7) 

199 F0RMAT(1H ,7F15.10) 

WR|TE(6, 198) TOFFB 

198 FORMAT ( 1 HO, 37X, * EXACT OFF TIME= ’,F15.10) 

WRITE{ 6,51 ) 1T1 

51 FORMAT [ 1 HO, 37X , 1 NO. OF ITERATIONS REQD.= ’,14) 

C 

C * CALCULATION OF THE MATRIX PS1 AND ITS EIGENVALUES * 
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DO 2000 I -I, 7 
DO 2010 IRICH=1,4 
DO 61 12=1,7 
XA( IZ)=X( 12) 

61 X8( IZ)=X{ IZ) • 

PX( | — C0H'( 1 RICH)*ABS(X( I ) } 

IF(ABS{X{ I ) ), LE.TOLL) DX( I J=G0N( 1RIGH) 

XA( I )=X( 1 )+DX{ I ) 

XBt I )=X( I ]^-DX{ i) 

TFFC1-T0FFB 
I T— 0 
T=TOFFB 

CALL STRAN( T, PH I F, DF, F2, G2, 7, 4,7, 4, 1 , 0, ERROR, D I FMAX, ITER, I FLSG ) 

1101 ZETA=0.0*XA(1 )+0.0*XA(2)+PHlF(7,3)«XA(3)+PHlFC7,4)* 

1XA(4)+PH I F( 7,5 )*XA(5)+1 . 0*XA( 7 }+0. O w U£ 1 )+DF{7, 2)*U(2)+ 

2 D F{ 7, 4 ) *U [ 4 ) -ET+PH I F( 7, 6 ) W XA ( 6 ) 

I F( ABS(ZETA) , LE, EPS) GO TO 1102 
I F{ IT.GE.NIT) GO TO 1003 
DTOFF=CON{ IRtCH)*TFFC1 
T=TFFC1+DT0FF 

CALL STRAN( T, PH I FT, DFT, F2, 02,7, 4, 7, 4, 1,0, ERROR, D l FMAX, ITER, I FLAG) 
ZETAN=0 . 0*XA(1 }+0 . 0*XA( 2 )+PH I FT( 7, 3 )*XA( 3 )+PH I FT( 7, 4 
1)*XA(4)+PHIFT(7,5)*XA{5)+1.0*XA{7)+0.0*U(1)*DFT(7,2)* 
2U(2)+DFT{7,4) a U{4)-ET+PHIFT(7,6)-*XA£6> 

DZETA==ZETAN-ZETA 

SLOPE=DZETA/DTOFF 

T Ff Cl =TF FC1 - { ZETA/SLO PE ) 

T=TFFC1 

CALL STRAN ( T, PH I F, DF, F2, G2, 7, 4, 7, 4, 1 , 0, ERROR, D 1 FMAX, ITER, I FLAG ) 

it=it+i 

GO TO 1101 

1102 DO 13 J=1 ,7 
TEMP"0.00 
DO 14 K=r,7 

DO 81 MENT=1 ,6 
31 PHI F(MENT,7)=0.0 
PH J F(7,7)=1 .0 
PHI F(7,1 )=0.0 
PHIFC7,2)=0.0 

14 TEMP=TEMPTPH t F( J, K)*XA( !() 

DO 15 1(5=1, 4 

DF( 7, 1)5=0.0 

15 TEMP=TEMP+DF{J,K)*U{K) 
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FUN7(J)=TEHP 
13 CONTINUE 

TFFC2=T0FFB 
I T=G 
T-TOFFB 

CALL STRAN (T, PH ! F, DF, F2, G2, 7, 4, 7, 4,1, 0, ERROR, D ! FMAX, ITER, J FLAG) 
•3103.ZETA=0.0*XI3PT J -0..0*XB(2)+PHIF(7,3)*XB{3) + PH|. , F(7,4)* 

1XBf4)+PHJF( f,5)*XB(5)-f-1.0*XB{7}+0.0*U(1}+DF(7,2}*UC2)+ 

2DF{7, 4)*U(4)-ET+P!II F(7, 6) Jf XB{6) 

IF(ABS(2ETA).LE.EPS) GO TO 11 04 
1 Ft IT.GE.N IT) GO TO 1003 
DTOFF=CON( I R I CH )*TFFC2 
T=TFFC2+DT0FF 

CALL STRAN ( T, PH I FT, OFT, F2, G2, 7, 4, 7, 4. 1,0, ERROR, D j FMAX, ITER, 1 FLAG) 
ZETAN=0 , 0«XB (7 )+0 . O w XB ( 2 ) +PH i FT( 7, 3 )*XB ( 3 )+PH I FT [ 7, 4 
• 1 )*XB£4)*PH! FTf7,5) # XB{5)+1 ,0*XB(7)+0. 0^0(1 )+DFT(7,2)* 
2U{2)*DFT(7,4)*U{4)“ET+PIHFTC7,6}*XB(6) 

DZETA=2ETAM-2ETA 

SL0PE=D2ETA/DT0FF 

TFFC2=TFFC2-{Z£TA/SL0PE) 

T=TFFC2 ' : 

CALL STRAN ( T, PH I F, DF, F2, G2, 7 , 4, 7, 4, 1 , 0, ERROR , D I FMAX, J TER, I FLAG ) 
IT=IT+1 
GO TO 1103 
1104 DO 62 J=1 ,7 
TEHP=0.00 
DO 63 K— 1,7 
DO 82 HENT=1,6 
82 PHIF(MENT,7)=0.0 
PHI F(7, 7)=1 .0 
PHI F(7, 1 }=0. 0 
PHIF(7,2)=0.0 

63 TEMP=TEMP+PHIF(J,KJ*XB{1<) 

DO 64 !(=1,4 
DF(7,1)=0.0 

64 TEMP=TEMPTDF(J,K)*U(K} 

FUN2(J>TEMP 

t/i CONTINUE 

DO 2011 J=1 , 7 

Loll RVAL(J, |RICH)=(FUN1{J)-FUN2(J))/(2.0 M DX{ 1)) 

2010 CONTINUE 

DO 65 J=1,7 

RVAL( J,5)=( RVAL( J, 2 1-0 . 25*RVAL( J, 1 ) )/0 > 75 
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RVALf J,6)=[RVAL( J,3)-0. 25*RVAL ( J , 2 ) ) /0 . 75 

65 RVAL{J,7}~{RVAL(J,4)-0.25*RVALfJ,3))/0,75 
DO 66 J=1 , 7 

RVAL{J,6)=(RVAL{J,6}-0.G625*RVAL{J,5))/0.9375 

66 RVAL{d;9)=(RVAL(J,7)~0.0625*RVAL(J,6)}/Q.9375 

DO 67 J~1 ,7 ' 

67 TEMP2{J, 1 )=!RVAL(J,9)-0.015625*RVAL(J,8) )/0. 984375 
2000 CONTINUE 

DO 17 1=1,7 
DO 17 J=1 ,7 
TEMP=0.00 
DO 19 K=1 ,7 

19 TEMP=TEMP+PH I N{ I , K)*TEMPH( K, J ) 

PS! { 1 ,J )“TEHP 
17 CONTINUE 

WRITE! 6, 92) 

92 FORMAT ( 1 HD, 37X, * MATRIX PSi= 1 ) 

DO 30 1=1,7 ' 

30 WR!TE{6,93)(PSl( l,J),J=1,7) 

93 FORMAT (7F1 5. 10) 

GO TO 1004 

1003 WRITE! 6, 100) I 

100 FORMAT!//, * CONVERGENCE NOT OBTAINED FOR X{f),l= *,14) 

1204 CONTINUE 

CALL ElGRFtPS1,7,7,2,W,2,7,WI<, IER) 

WRITE! 6,94) . 

94 FORMAT(1KO,37X, * THE EIGENVALUES ARE— * ) 

WRITE(6,103)(W{ I ), 1=1,7) 

103 FORMAT (1 H ,7F15,10) 

WRITE! 6,104) WK(1 ), i ER - ■ 

104 FORMAT! 1 HO, 37X, ’ CONVERGENCE TOLERANCE^ , F15* 10, * IER= *,I4) 
905 CONTINUE 

STOP 
END 

SUBR0UT1 NE STRAN( TAU, PH I, THETA, A, B, NA, NB,MA, MB, MODE, NTERMS, 

1 TOL, D I FMAX, I TER, i FLAG ) 

REAL PHI ( NA, NA),THETA(NA,'NB), A{NA, NA},B£ NA,NB), 

1 W0RK1 (10,10), W0RK2 (10,10), DUMMY! 1,1), FACT , CON 1 , FAC2, C0N2 
C1=1 
TS=TAU 
DO 4 1-1, MA 




r- 




jtwwpw] 


l 


ORIGINAL PASS Cfc 
OF POOR QUALITY 


oo 


V' 


1 • 




DO 2 J=1,HA 
WORK1((,J)=0 
W0RI<2( 1 , J )=0 
I’ll f £ 1 , J )=0 
2 CONTINUE 
W0RK1{ l,l)=1 
W0RK2( I , I )=TAU 
PH t C I , I )-1 
4 CONTINUE 
DO 6 I =1 , MA 
DO 6 J-1 ,MB 
6 ^HETAt l,J )=0 
0 1 FMAX=1 . E8 
ND0=50 

1F£NTERMS.GT.0J ndo=nterms 
FAC1-1 
i : IFLAG=0 

■,! DO TOGO 1=1, NDO 

CALL MXMUL( DUMMY, W0RK1 , A, TS, 1 , 1 , TO, TO, NA, NA, HA, MA,MA,MA, T ) 

FACT=FACT*I 

CON T=T./ FACT 

TF(NTERMS. EQ.O.AND. I .GE.4) CALL SERRORf PHI, WORK!, CONI, NA,NA, 

1 TO, 10,MA,HA, DIFMAX) 

CALL MXADD( DUMMY, Pill , WORK! , Cl , CONI , T , T , NA, NA, TO, 10, MA, MA, 1 ) 
IF(MuDE.EQ.2) GO TO 500 
FAC2=FACT*( 1+1) 

C0N2=TAU/FAC2 

CALL MXADDf DUMMY, W0RK2, WORK! , CT ,CON2, T , 1, TO, TO, TO, T 0, MA, MA, T ) 
500 CONTINUE 
ITER- I 

I F(NTERMS,GT.£>.OR, I .LT.4J GO TO 1000 
: IF(D( FMAX.LE.TOL) GO TO TTOO 
1000 CONTINUE 
1100 CONTINUE 

IF(MODE.EQ.I) CALL MXMUL£THETA,W0RK2,B,C1,NA,NB, 10, 10,NA,NB, 

1 MA, MA, MA, MB , 2 ) 

J F{ I TER. EQ. NDO, AND.NTERMS, EQ.OJ I FLAG=1 
RETURN 
END 

SUBROUTINE SERROR[ AMX, BMX, CCC, IA,JA, i DO, JDO, D I FMAX) 
DIMENSION AMX( IA,JA),BMX( IB,JB) 


to 

tn 
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01 FHAX=1 -E-30 
DO TOD 1=1, J DO 
DO 50 J=1,JD0 

J F( AMX( f , J) . EQ.0.0 ) GO TO 50 
CHANGE^ ABS( BMX( I , J )*CCC/AMX( I , J ) ) 

I F ( CHANG E . GT . D I FMAX } D J FMAX=CHANGE 
50 CONTINUE 
100 CONTINUE 
RETURN 
END 


SUBROUTINE MXADD(RMX,AMX,BMX J! AGG,BCC, IR,JR, 1A,JA, 18,06, IDO,JDO, 
1 MODE) 

DIMENSION AMX( lA,JA),8MXf I B, JB),RMX{ IR,JR) 

IF( IA.LT. IOO.OR.JA.LT.JQG) GO TO 999 
JF( IB.LT. IDO.OR.JB.LT.JDO) GO TO 999 
GO TO (10, TOO), MODE 
10 CONTINUE 

DO 50 1=1, IDO 
DO 50 J=1,JD0 

50 AMX( i,J)=AMX( 1,J)*ACC+BMX( l,J)*BCC 
GO TO 300 
100 CONTINUE 

I F{ IR. L.T. I DO. OR. JR. LT. JDO) GO TO 999 
DO 200 1=1, IDO 
DO 200 J=1, JDO 

200 RMX( f,J)=AMX( l, J ) # ACC+BMX( l,J)*BCC 
300 RETURN 
999 CONTINUE 
RETURN 
END 


SUBROUTINE MXMULf RMX, AMX, BMX, CCC, IR,JR, IA,JA, IB,JB, IDA,JDA, IDB,JDB 
1 ,MODE) 

DIMENSION AWX{ IA,JA),BMX( IB,JB),RMX(IR,JR),TEMP(20) 

IF( 1DA*IDB«JDA*JDB.GT. I A* I B*JA*JB) GO TO 999 
IF(JDB.GT.JDA) GO TO 999 
GO TO (10, 210), MODE 
10 DO 100 1=1, IDA 
DO 20 L=1 , JDA 
20 TEMPf L)=AMX( I , L) 
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DO 80 J— 1 , JDB 
SU[«t=0 

DO 40 K=1,JDA 

^0 SUM=SUM+TEHP(K)«BHX(K,JJ 
AMX( l,JJ=SUM*CCC 

ao continue: 

100 CONTINUE 
GO TO 800 
210 CONTINUE 

I F( IN. LT. IDA. OR.JR.LT, JDB) GO TO 999 
DO i|00 1=1, IDA 
DO 380 J=1 , JDB 
SUM=0 

DO 3NO K=1, JDA 

3<i0 SUM=SUM+AMX( I,K)*BMX(K,J) 

RMX{ l,J]=SUM*CCC 
380 CONTINUE 
IlOQ CONTINUE 
600 RETURN 
999 CONTINUE 
RETURN 

END to 

U1 
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